Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana

Plant Immunity: From Signaling to Epigenetic Control of Defense

Pathogen recognition by plants results in the activation of signaling pathways that induce defense reactions. There is growing evidence indicating that epigenetic mechanisms directly participate in plant immune memory. Here, we discuss current knowledge of diverse epigenomic processes and elements, such as noncoding RNAs, DNA and RNA methylation, histone post-translational modifications, and chromatin remodeling, that have been associated with the regulation of immune responses in plants. Furthermore, we discuss the currently limited evidence of transgenerational inheritance of pathogen-induced defense priming, together with its potentials, challenges, and limitations for crop improvement and biotechnological applications.

Regulation mechanism of long non-coding RNA in plant response to stress

Long non-coding RNA (lncRNA) is a non-coding RNA greater than 200 nucleotides in length. LncRNAs can regulate gene expression at transcription and post-transcription, epigenetic level, and plays an important role in a wide range of biological processes such as genomic imprinting, chromatin remodeling, transcriptional activation, transcriptional interference and cell cycle. It becomes the current hot topics in the study of molecular biology and genetics. Emerging evidence proposed that lncRNAs play important roles in response to both abiotic and biotic stress. In this review, we discuss the role of lncRNAs in drought resistance, salt resistance, disease resistance, and immunity of plants, providing strong evidence for exploring the important role of lncRNAs in plant resistance, in order to explore new ideas and new targets for prevention and control.

Analysis of long-non-coding RNAs associated with ethylene in tomato

Long-Non-Coding RNAs (LncRNAs) are a class of non-coding endogenous RNAs contributing to numerous biological processes. LeERF1 is a tomato ethylene response factor (ERF) near the end of the ethylene signal transduction pathway. To identify lncRNAs in tomato and elucidate their roles in ethylene signaling, deep sequencing was deployed in over-expression and repression LeERF1 transgenic and control tomato fruits. A total of 397 lncRNAs were identified, including 169 tomato lncRNAs that had not previously been identified. Among these, 12 were differentially expressed between the transgenic and control tomato fruits. Numerous lncRNA target genes were identified including many associated with ethylene signaling including auxin response factors and auxin-induced proteins, F-box proteins, ERFs and MADS-box proteins. In addition, two lncRNAs were found to be the precursor of three miRNAs and four lncRNAs could be targeted by five miRNAs. We propose a regulatory model highlighting the relationships between lncRNAs and their targets involved in ethylene signal transduction which establishes a foundation for addressing the role of LncRNAs in ethylene response.

Genome-wide analysis of long non-coding RNAs affecting roots development at an early stage in the rice response to cadmium stress

Background

Long non-coding RNAs (lncRNAs) have been found to play a vital role in several gene regulatory networks involved in the various biological processes in plants related to stress response. However, systematic analyses of lncRNAs expressed in rice Cadmium (Cd) stress are seldom studied. Thus, we presented the characterization and expression of lncRNAs in rice root development at an early stage in response to Cd stress.

Results

The lncRNA deep sequencing revealed differentially expressed lncRNAs among Cd stress and normal condition. In the Cd stress group, 69 lncRNAs were up-regulated and 75 lncRNAs were down-regulated. Furthermore, 386 matched lncRNA-mRNA pairs were detected for 120 differentially expressed lncRNAs and 362 differentially expressed genes in cis, and target gene-related pathway analyses exhibited significant variations in cysteine and methionine metabolism pathway-related genes. For the genes in trans, overall, 28,276 interaction relationships for 144 lncRNAs and differentially expressed protein-coding genes were detected. The pathway analyses found that secondary metabolites, such as phenylpropanoids and phenylalanine, and photosynthesis pathway-related genes were significantly altered by Cd stress. All of these results indicate that lncRNAs may regulate genes of cysteine-rich peptide metabolism in cis, as well as secondary metabolites and photosynthesis in trans, to activate various physiological and biochemical reactions to respond to excessive Cd.

Conclusion

The present study could provide a valuable resource for lncRNA studies in response to Cd treatment in rice. It also expands our knowledge about lncRNA biological function and contributes to the annotation of the rice genome.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4807-6) contains supplementary material, which is available to authorized users.

Re-analysis of long non-coding RNAs and prediction of circRNAs reveal their novel roles in susceptible tomato following TYLCV infection

BACKGROUND

Long Noncoding-RNAs (LncRNAs) are known to be involved in some biological processes, but their roles in plant-virus interactions remain largely unexplored. While circular RNAs (circRNAs) have been studied in animals, there has yet to be extensive research on them in a plant system, especially in tomato-tomato yellow leaf curl virus (TYLCV) interaction.

RESULTS

In this study, RNA transcripts from the susceptible tomato line JS-CT-9210 either infected with TYLCV or untreated, were sequenced in a pair-end strand-specific manner using ribo-zero rRNA removal library method. A total of 2056 lncRNAs including 1767 long intergenic non-coding RNA (lincRNAs) and 289 long non-coding natural antisense transcripts (lncNATs) were obtained. The expression patterns in lncRNAs were similar in susceptible tomato plants between control check (CK) and TYLCV infected samples. Our analysis suggested that lncRNAs likely played a role in a variety of functions, including plant hormone signaling, protein processing in the endoplasmic reticulum, RNA transport, ribosome function, photosynthesis, glulathione metabolism, and plant-pathogen interactions. Using virus-induced gene silencing (VIGS) analysis, we found that reduced expression of the lncRNA S-slylnc0957 resulted in enhanced resistance to TYLCV in susceptible tomato plants. Moreover, we identified 184 circRNAs candidates using the CircRNA Identifier (CIRI) software, of which 32 circRNAs were specifically expressed in untreated samples and 83 circRNAs in TYLCV samples. Approximately 62% of these circRNAs were derived from exons. We validated the circRNAs by both PCR and Sanger sequencing using divergent primers, and found that most of circRNAs were derived from the exons of protein coding genes. The silencing of these circRNAs parent genes resulted in decreased TYLCV virus accumulation.

CONCLUSION

In this study, we identified novel lncRNAs and circRNAs using bioinformatic approaches and showed that these RNAs function as negative regulators of TYLCV infection. Moreover, the expression patterns of lncRNAs in susceptible tomato plants were different from that of resistant tomato plants, while exonic circRNAs expression positively associated with their respective protein coding genes. This work provides a foundation for elaborating the novel roles of lncRNAs and circRNAs in susceptible tomatoes following TYLCV infection.

Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton

Long noncoding RNAs (lncRNAs) have several known functions in plant development, but their possible roles in responding to plant disease remain largely unresolved. In this study, we described a comprehensive disease-responding lncRNA profiles in defence against a cotton fungal disease Verticillium dahliae. We further revealed the conserved and specific characters of disease-responding process between two cotton species. Conservatively for two cotton species, we found the expression dominance of induced lncRNAs in the Dt subgenome, indicating a biased induction pattern in the co-existing subgenomes of allotetraploid cotton. Comparative analysis of lncRNA expression and their proposed functions in resistant Gossypium barbadense cv. '7124' versus susceptible Gossypium hirsutum cv. 'YZ1' revealed their distinct disease response mechanisms. Species-specific (LS) lncRNAs containing more SNPs displayed a fiercer inducing level postinfection than the species-conserved (core) lncRNAs. Gene Ontology enrichment of LS lncRNAs and core lncRNAs indicates distinct roles in the process of biotic stimulus. Further functional analysis showed that two core lncRNAs, GhlncNAT-ANX2- and GhlncNAT-RLP7-silenced seedlings, displayed an enhanced resistance towards V. dahliae and Botrytis cinerea, possibly associated with the increased expression of LOX1 and LOX2. This study represents the first characterization of lncRNAs involved in resistance to fungal disease and provides new clues to elucidate cotton disease response mechanism.

Genome-wide profiling of long non-coding RNAs from tomato and a comparison with mRNAs associated with the regulation of fruit ripening

BACKGROUND

Long non-coding RNAs (lncRNAs) are involved in multiple biological processes in both mammals and plants. There is growing evidence that they are associated with development; but their expression and regulation during fruit ripening in the model plant tomato (Solanum lycopersicum) has yet to be described.

RESULTS

Following integration of 134 RNA-seq data sets, we identified 79,322 putative lncRNAs, consisting of 70,635 lincRNAs, 8085 antisense non-coding RNAs (ancRNAs) and 602 sense lncRNAs (slncRNAs). lncRNAs had specific features that are distinct from mRNAs, including tissue-specificity, and shorter and fewer exons. Notably, more than 5000 of the novel lincRNAs were found to be expressed across the mature green (MG), breaker (BR) and breaker plus 7 days (BR + 7) developmental stages. The differently expressed lincRNAs had different DNA methylation profiles from the mRNAs.

CONCLUSIONS

Integrating transcriptome datasets and genome-wide screening enabled the identification of a comprehensive set of tomato lncRNAs. Here, we found that the lncRNAs DNA methylation profiles were different from those of mRNAs. This will help future investigation of lncRNA function, especially for the dissection of the molecular mechanisms involved in the regulation of fruit development.

Long non-coding RNAs of switchgrass (Panicum virgatum L.) in multiple dehydration stresses

BACKGROUND

Long non-coding RNAs (lncRNAs) play important roles in plant growth and stress responses. Studies of lncRNAs in non-model plants are quite limited, especially those investigating multiple dehydration stresses. In this study, we identified novel lncRNAs and analyzed their functions in dehydration stress memory in switchgrass, an excellent biofuel feedstock and soil-conserving plant in the Gramineae family.

RESULTS

We analyzed genome-wide transcriptional profiles of leaves of 5-week-old switchgrass plantlets grown via tissue culture after primary and secondary dehydration stresses (D1 and D2) and identified 16,551 novel lncRNAs, including 4554 annotated lncRNAs (targeting 3574 genes), and 11,997 unknown lncRNAs. Gene ontology and pathway enrichment analysis of annotated genes showed that the differentially expressed lncRNAs were related to abscisic acid (ABA) and ethylene (ETH) biosynthesis and signal transduction, and to starch and sucrose metabolism. The upregulated lncRNAs and genes were related to ABA synthesis and its signal transduction, and to trehalose synthesis. Meanwhile, lncRNAs and genes related to ETH biosynthesis and signal transduction were suppressed. LncRNAs and genes involved in ABA metabolism were verified using quantitative real-time PCR, and the endogenous ABA content was determined via high performance liquid chromatography mass spectrometry (HPLC-MS). These results showed that ABA accumulated significantly during dehydration stress, especially in D2. Furthermore, we identified 307 dehydration stress memory lncRNAs, and the ratios of different memory types in switchgrass were similar to those in Arabidopsis and maize.

CONCLUSIONS

The molecular responses of switchgrass lncRNAs to multiple dehydration stresses were researched systematically, revealing novel information about their transcriptional regulatory behavior. This study provides new insights into the response mechanism to dehydration stress in plants. The lncRNAs and pathways identified in this study provide valuable information for genetic modification of switchgrass and other crops.

A novel antisense long noncoding RNA, TWISTED LEAF, maintains leaf blade flattening by regulating its associated sense R2R3-MYB gene in rice

Natural antisense long noncoding RNAs (lncRNAs) are widespread in many organisms. However, their biological functions remain largely unknown, particularly in plants. We report the identification and characterization of an endogenous lncRNA, TWISTED LEAF (TL), which is transcribed from the opposite strand of the R2R3 MYB transcription factor gene locus, OsMYB60, in rice (Oryza sativa). TL and OsMYB60 were found to be coexpressed in many different tissues, and the expression level of TL was higher than that of OsMYB60. Downregulation of TL by RNA interference (RNAi) and overexpression of OsMYB60 resulted in twisted leaf blades in transgenic rice. The expression level of OsMYB60 was significantly increased in TL-RNAi transgenic plants. This suggests that TL may play a cis-regulatory role on OsMYB60 in leaf morphological development. We also determined that the antisense transcription suppressed the sense gene expression by mediating chromatin modifications. We further discovered that a C2H2 transcription factor, OsZFP7, is an OsMYB60 binding partner and involved in leaf development. Taken together, these findings reveal that the cis-natural antisense lncRNA plays a critical role in maintaining leaf blade flattening in rice. Our study uncovers a regulatory mechanism of lncRNA in plant leaf development.

Landscape of long non-coding RNAs in Trichophyton mentagrophytes-induced rabbit dermatophytosis lesional skin and normal skin

Emerging evidences suggest that long non-coding RNAs (lncRNAs) play important role in disease development. However, the role of rabbit lncRNAs in the pathogenesis of dermatophytosis remains elusive. The present study aimed to study and characterize lncRNA transcriptome in 8 T. mentagrophytes-induced female rabbit dermatophytosis lesional (TM) and 4 normal saline-infected (NS) skin biopsies using RNAseq. We identified 5883 lncRNAs in 12 strand-specific RNA-seq libraries and found 64 differentially expressed lncRNAs (q < 0.05) in TM relative to NS. As in other mammalian counterparts, rabbit lncRNAs were distributed in all chromosomes except the Y chromosome and were generally smaller in size and fewer in exon numbers compared to protein coding genes. Next, co-expression analysis revealed that 107 pairs between 32 DE lncRNAs and 96 protein coding genes showed a highly correlated expression (|r| > 0.8). Moreover, miRPara analysis of the lncRNAs revealed 173 lncRNAs with precursor sequences for 9561 probable novel miRNAs. Finally, q-PCR results validated the RNA-seq results with eight randomly selected lncRNAs. To the best of our knowledge, this is the first report on rabbit lncRNAs, and our results highlighted the potential role of lncRNAs in the pathogenesis of dermatophytosis.

Long Non-Coding RNAs Responsive to Salt and Boron Stress in the Hyper-Arid Lluteño Maize from Atacama Desert

Long non-coding RNAs (lncRNAs) have been defined as transcripts longer than 200 nucleotides, which lack significant protein coding potential and possess critical roles in diverse cellular processes. Long non-coding RNAs have recently been functionally characterized in plant stress-response mechanisms. In the present study, we perform a comprehensive identification of lncRNAs in response to combined stress induced by salinity and excess of boron in the Lluteño maize, a tolerant maize landrace from Atacama Desert, Chile. We use deep RNA sequencing to identify a set of 48,345 different lncRNAs, of which 28,012 (58.1%) are conserved with other maize (B73, Mo17 or Palomero), with the remaining 41.9% belonging to potentially Lluteño exclusive lncRNA transcripts. According to B73 maize reference genome sequence, most Lluteño lncRNAs correspond to intergenic transcripts. Interestingly, Lluteño lncRNAs presents an unusual overall higher expression compared to protein coding genes under exposure to stressed conditions. In total, we identified 1710 putatively responsive to the combined stressed conditions of salt and boron exposure. We also identified a set of 848 stress responsive potential trans natural antisense transcripts (trans-NAT) lncRNAs, which seems to be regulating genes associated with regulation of transcription, response to stress, response to abiotic stimulus and participating of the nicotianamine metabolic process. Reverse transcription-quantitative PCR (RT-qPCR) experiments were performed in a subset of lncRNAs, validating their existence and expression patterns. Our results suggest that a diverse set of maize lncRNAs from leaves and roots is responsive to combined salt and boron stress, being the first effort to identify lncRNAs from a maize landrace adapted to extreme conditions such as the Atacama Desert. The information generated is a starting point to understand the genomic adaptabilities suffered by this maize to surpass this extremely stressed environment.

Roles of Non-Coding RNA in Sugarcane-Microbe Interaction

Studies have highlighted the importance of non-coding RNA regulation in plant-microbe interaction. However, the roles of sugarcane microRNAs (miRNAs) in the regulation of disease responses have not been investigated. Firstly, we screened the sRNA transcriptome of sugarcane infected with Acidovorax avenae. Conserved and novel miRNAs were identified. Additionally, small interfering RNAs (siRNAs) were aligned to differentially expressed sequences from the sugarcane transcriptome. Interestingly, many siRNAs aligned to a transcript encoding a copper-transporter gene whose expression was induced in the presence of A. avenae, while the siRNAs were repressed in the presence of A. avenae. Moreover, a long intergenic non-coding RNA was identified as a potential target or decoy of miR408. To extend the bioinformatics analysis, we carried out independent inoculations and the expression patterns of six miRNAs were validated by quantitative reverse transcription-PCR (qRT-PCR). Among these miRNAs, miR408-a copper-microRNA-was downregulated. The cleavage of a putative miR408 target, a laccase, was confirmed by a modified 5'RACE (rapid amplification of cDNA ends) assay. MiR408 was also downregulated in samples infected with other pathogens, but it was upregulated in the presence of a beneficial diazotrophic bacteria. Our results suggest that regulation by miR408 is important in sugarcane sensing whether microorganisms are either pathogenic or beneficial, triggering specific miRNA-mediated regulatory mechanisms accordingly.

Long noncoding RNAs that respond to Fusarium oxysporum infection in 'Cavendish' banana (Musa acuminata)

Long noncoding RNAs (lncRNAs) are a class of genes that influence a variety of biological functions through acting as signal, decoy, guide, and scaffold molecules. In banana (Musa spp.), an important economic fruit crop, particularly in Southeast Asia, the wilt disease caused by Fusarium oxysporum f. sp. cubense (Foc), especially strain Foc TR4, is disastrous. In banana, how the biogenesis of these lncRNAs is regulated in response to pathogen infection is still largely unknown. In this study, strand-specific paired-end RNA sequencing of banana samples was performed on susceptible and resistant cultivars inoculated with Foc, with three biological replicates and at two different times after infection. Overall, 5,294 lncRNAs were predicted with high confidence through strict filtration, including long intergenic ncRNA (lincRNA) and antisense lncRNA. Differentially expressed (DE) lncRNAs were identified in response to Foc infection in the inoculated versus the mock-inoculated banana of the susceptible 'BX' and resistant 'NK' cultivars. Through KEGG, GO, and the expression levels of the DE lncRNAs, some DE lncRNAs were predicted to be involved in plant-pathogen interactions and phytohormone signal transduction. In this study, this catalog of lncRNAs and their properties will facilitate further experimental studies and functional classifications of these genes.

Seed weight differences between wild and domesticated soybeans are associated with specific changes in gene expression

KEY MESSAGE

Our study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication. In addition, potential contributions of lncRNAs to seed weight were also investigated. Soybeans have a long history of domestication in China, and there are several significant phenotypic differences between cultivated and wild soybeans, for example, seeds of cultivars are generally larger and heavier than those from wild accessions. We analyzed seed transcriptomes from thirteen soybean samples, including six landraces and seven wild accessions using strand-specific RNA sequencing. Differentially expressed genes related to seed weight were identified, and some of their homologs were associated with seed development in Arabidopsis. We also identified 1251 long intergenic noncoding RNAs (lincRNAs), 243 intronic RNAs and 81 antisense lncRNAs de novo from these soybean transcriptomes. We then profiled the expression patterns of lncRNAs in cultivated and wild soybean seeds, and found that transcript levels of a number of lncRNAs were sample-specific. Moreover, gene transcript and lincRNA co-expression network analysis showed that some soybean lincRNAs might have functional roles as they were hubs of co-expression modules. In conclusion, this study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication, and will provide a useful future resource for molecular breeding of soybeans.

Whole transcriptome sequencing of Pseudomonas syringae pv. actinidiae-infected kiwifruit plants reveals species-specific interaction between long non-coding RNA and coding genes

An outbreak of kiwifruit bacterial canker disease caused by Pseudomonas syringae pv. actinidiae (Psa) beginning in 2008 caused disaster to the kiwifruit industry. However the mechanisms of interaction between kiwifruit and Psa are unknown. Long noncoding RNAs (lncRNAs) are known to regulate many biological processes, but comprehensive repertoires of kiwifruit lncRNAs and their effects on the interaction between kiwifruit and Psa are unknown. Here, based on in-depth transcriptomic analysis of four kiwifruit materials at three stages of infection with Psa, we identified 14,845 transcripts from 12,280 loci as putative lncRNAs. Hierarchical clustering analysis of differentially-expressed transcripts reveals that both protein-coding and lncRNA transcripts are expressed species-specifically. Comparing differentially-expressed transcripts from different species, variations in pattern-triggered immunity (PTI) were the main causes of species-specific responses to infection by Psa. Using weighted gene co-expression network analysis, we identified species-specific expressed key lncRNAs which were closely related to plant immune response and signal transduction. Our results illustrate that different kiwifruit species employ multiple different plant immunity layers to fight against Psa infection, which causes distinct responses. We also discovered that lncRNAs might affect kiwifruit responses to Psa infection, indicating that both protein-coding regions and noncoding regions can affect kiwifruit response to Psa infection.

Dehydration stress extends mRNA 3' untranslated regions with noncoding RNA functions in Arabidopsis

The 3′ untranslated regions (3′ UTRs) of mRNAs play important roles in the regulation of mRNA localization, translation, and stability. Alternative cleavage and polyadenylation (APA) generates mRNAs with different 3′ UTRs, but the involvement of this process in stress response has not yet been clarified. Here, we report that a subset of stress-related genes exhibits 3′ UTR extensions of their mRNAs during dehydration stress. These extended 3′ UTRs have characteristics of long noncoding RNAs and likely do not interact with miRNAs. Functional studies using T-DNA insertion mutants reveal that they can act as antisense transcripts to repress expression levels of sense genes from the opposite strand or can activate the transcription or lead to read-through transcription of their downstream genes. Further analysis suggests that transcripts with 3′ UTR extensions have weaker poly(A) signals than those without 3′ UTR extensions. Finally, we show that their biogenesis is partially dependent on a trans-acting factor FPA. Taken together, we report that dehydration stress could induce transcript 3′ UTR extensions and elucidate a novel function for these stress-induced 3′ UTR extensions as long noncoding RNAs in the regulation of their neighboring genes.

Emerging roles of long non-coding RNAs in plant response to biotic and abiotic stresses

Spectacular progress in high-throughput transcriptome sequencing and expression profiling using next-generation sequencing technologies have recently revolutionized molecular biology and allowed massive advances in identifying the genomic regions and molecular mechanisms underlying transcriptional regulation of growth, development, and stress response. Through recent research, non-coding RNAs, in particular long non-coding RNAs, have emerged as key regulators of transcription in eukaryotes. Long non-coding RNAs are vastly heterogeneous groups of RNAs that execute a broad range of essential roles in various biological processes at the epigenetic, transcriptional, and post-transcriptional levels. They modulate transcription through diverse mechanisms. Recently, numerous lncRNAs have been identified to be associated with defense responses to biotic and abiotic stresses. These have been suggested to perform indispensable roles in plant immunity and adaptation to environmental conditions. However, only a few lncRNAs have been functionally characterized in plants. In this paper, we summarize the present knowledge of lncRNAs, review the recent advances in understanding regulatory functions of lncRNAs, and highlight the emerging roles of lncRNAs in regulating immune responses in plants.

Transposable elements (TEs) contribute to stress-related long intergenic noncoding RNAs in plants

Noncoding RNAs have been extensively described in plant and animal transcriptomes by using high-throughput sequencing technology. Of these noncoding RNAs, a growing number of long intergenic noncoding RNAs (lincRNAs) have been described in multicellular organisms, however the origins and functions of many lincRNAs remain to be explored. In many eukaryotic genomes, transposable elements (TEs) are widely distributed and often account for large fractions of plant and animal genomes yet the contribution of TEs to lincRNAs is largely unknown. By using strand-specific RNA-sequencing, we profiled the expression patterns of lincRNAs in Arabidopsis, rice and maize, and identified 47 611 and 398 TE-associated lincRNAs (TE-lincRNAs), respectively. TE-lincRNAs were more often derived from retrotransposons than DNA transposons and as retrotransposon copy number in both rice and maize genomes so did TE-lincRNAs. We validated the expression of these TE-lincRNAs by strand-specific RT-PCR and also demonstrated tissue-specific transcription and stress-induced TE-lincRNAs either after salt, abscisic acid (ABA) or cold treatments. For Arabidopsis TE-lincRNA11195, mutants had reduced sensitivity to ABA as demonstrated by longer roots and higher shoot biomass when compared to wild-type. Finally, by altering the chromatin state in the Arabidopsis chromatin remodelling mutant ddm1, unique lincRNAs including TE-lincRNAs were generated from the preceding untranscribed regions and interestingly inherited in a wild-type background in subsequent generations. Our findings not only demonstrate that TE-associated lincRNAs play important roles in plant abiotic stress responses but lincRNAs and TE-lincRNAs might act as an adaptive reservoir in eukaryotes.

Present Scenario of Long Non-Coding RNAs in Plants

Small non-coding RNAs have been extensively studied in plants over the last decade. In contrast, genome-wide identification of plant long non-coding RNAs (lncRNAs) has recently gained momentum. LncRNAs are now being recognized as important players in gene regulation, and their potent regulatory roles are being studied comprehensively in eukaryotes. LncRNAs were first reported in humans in 1992. Since then, research in animals, particularly in humans, has rapidly progressed, and a vast amount of data has been generated, collected, and organized using computational approaches. Additionally, numerous studies have been conducted to understand the roles of these long RNA species in several diseases. However, the status of lncRNA investigation in plants lags behind that in animals (especially humans). Efforts are being made in this direction using computational tools and high-throughput sequencing technologies, such as the lncRNA microarray technique, RNA-sequencing (RNA-seq), RNA capture sequencing, (RNA CaptureSeq), etc. Given the current scenario, significant amounts of data have been produced regarding plant lncRNAs, and this amount is likely to increase in the subsequent years. In this review we have documented brief information about lncRNAs and their status of research in plants, along with the plant-specific resources/databases for information retrieval on lncRNAs.

Comparative transcriptome analysis between resistant and susceptible tomato allows the identification of lncRNA16397 conferring resistance to Phytophthora infestans by co-expressing glutaredoxin

The rapid development of omics sequencing technology has facilitated the identification of thousands of long non-coding (lnc)RNAs in plant species, but the role of lncRNAs in plant-pathogen interactions remains largely unexplored. We used comparative transcriptome analysis of Phytophthora infestans-resistant and -susceptible tomatoes to identify differentially expressed genes (DEGs) and lncRNAs (DELs), and examine lncRNA-mRNA networks. A total of 1037 DEGs and 688 DELs were identified between P. infestans-resistant and -susceptible tomatoes. The co-localization networks, including 128 DEGs and 127 DELs, were performed. We found that lncRNA16397 acted as an antisense transcript of SlGRX22 to regulate its expression, and also induced SlGRX21 expression when lncRNA16397 was overexpressed. In addition, disease symptoms and reactive oxygen species (ROS) accumulation in tomatoes overexpressing lncRNA16397 and SpGRX were fewer and lower than those in wild-type after P. infestans infection. This result suggests that tomato lncRNA16397 induces SlGRX expression to reduce ROS accumulation and alleviate cell membrane injury, resulting in enhanced resistance to P. infestans. Our results provide insight into lncRNAs involved in the response of tomato to P. infestans infection, demonstrate that the lncRNA16397-GRXs network is an important component of the P. infestans network in tomato, and provide candidates for breeding to enhance biotic stress-resistance in tomato.

Massive Effect on LncRNAs in Human Monocytes During Fungal and Bacterial Infections and in Response to Vitamins A and D

Mycoses induced by C.albicans or A.fumigatus can cause important host damage either by deficient or exaggerated immune response. Regulation of chemokine and cytokine signaling plays a crucial role for an adequate inflammation, which can be modulated by vitamins A and D. Non-coding RNAs (ncRNAs) as transcription factors or cis-acting antisense RNAs are known to be involved in gene regulation. However, the processes during fungal infections and treatment with vitamins in terms of therapeutic impact are unknown. We show that in monocytes both vitamins regulate ncRNAs involved in amino acid metabolism and immune system processes using comprehensive RNA-Seq analyses. Compared to protein-coding genes, fungi and bacteria induced an expression change in relatively few ncRNAs, but with massive fold changes of up to 4000. We defined the landscape of long-ncRNAs (lncRNAs) in response to pathogens and observed variation in the isoforms composition for several lncRNA following infection and vitamin treatment. Most of the involved antisense RNAs are regulated and positively correlated with their sense protein-coding genes. We investigated lncRNAs with stimulus specific immunomodulatory activity as potential marker genes: LINC00595, SBF2-AS1 (A.fumigatus) and RP11-588G21.2, RP11-394l13.1 (C.albicans) might be detectable in the early phase of infection and serve as therapeutic targets in the future.

The Role of Long Noncoding RNAs in Plant Stress Tolerance

Plants must adapt to multiple biotic and abiotic stresses ; thus, sensing and responding to environmental signals is imperative for their survival. Moreover, understanding these responses is imperative for efforts to improve plant yield and consistency. Regulation of transcript levels is a key aspect of the plant response to environmental signals. Long noncoding RNAs (lncRNAs) have gained widespread attention in recent years with the advance of high-throughput sequencing technologies. As important biological regulators, lncRNAs have been implicated in a wide range of developmental processes and diseases in animals. However, knowledge of the role that lncRNAs play in plant stress tolerance remains limited. Here, we review recent studies on the identification, characteristics, classification, and biological functions of lncRNAs in response to various stresses, including bacterial pathogens, excess light, drought, salinity, hypoxia, extreme temperatures, and nitrogen/phosphate deficiency. We also discuss possible directions for future research.

CMV2b-AGO Interaction Is Required for the Suppression of RDR-Dependent Antiviral Silencing in Arabidopsis

Using a transient plant system, it was previously found that the suppression of Cucumber mosaic virus (CMV) 2b protein relies on its double-strand (ds) RNA binding capacity, but it is independent of its interaction with ARGONAUTE (AGO) proteins. Thus, the biological meaning of the 2b-AGO interaction in the context of virus infection remains elusive. In this study, we created infectious clones of CMV mutants that expressed the 2b functional domains of dsRNA or AGO binding and tested the effect of these CMV mutants on viral pathogenicity. We found that the mutant CMV2b_(1–76) expressing the 2b dsRNA-binding domain exhibited the same virulence as wild-type CMV in infection with either wild-type Arabidopsis or rdr1/6 plants with RDR1- and RDR6-deficient mutations. However, remarkably reduced viral RNA levels and increased virus (v)siRNAs were detected in CMV2b_(1–76)-infected Arabidopsis in comparison to CMV infection, which demonstrated that the 2b_(1–76) deleted AGO-binding domain failed to suppress the RDR1/RDR6-dependent degradation of viral RNAs. The mutant CMV2b_(8–111) expressing mutant 2b, in which the N-terminal 7 amino acid (aa) was deleted, exhibited slightly reduced virulence, but not viral RNA levels, in both wild-type and rdr1/6 plants, which indicated that 2b retained the AGO-binding activity acquired the counter-RDRs degradation of viral RNAs. The deletion of the N-terminal 7 aa of 2b affected virulence due to the reduced affinity for long dsRNA. The mutant CMV2b_(18–111) expressing mutant 2b lacked the N-terminal 17 aa but retained its AGO-binding activity greatly reduced virulence and viral RNA level. Together with the instability of both 2b_(18–111)-EGFP and RFP-AGO4 proteins when co-expressed in Nicotiana benthamiana leaves, our data demonstrates that the effect of 2b-AGO interaction on counter-RDRs antiviral defense required the presence of 2b dsRNA-binding activity. Taken together, our findings demonstrate that the dsRNA-binding activity of the 2b was essential for virulence, whereas the 2b-AGO interaction was necessary for interference with RDR1/6-dependent antiviral silencing in Arabidopsis.

Genome-wide identification of potato long intergenic noncoding RNAs responsive to Pectobacterium carotovorum subspecies brasiliense infection

BACKGROUND

Long noncoding RNAs (lncRNAs) represent a class of RNA molecules that are implicated in regulation of gene expression in both mammals and plants. While much progress has been made in determining the biological functions of lncRNAs in mammals, the functional roles of lncRNAs in plants are still poorly understood. Specifically, the roles of long intergenic nocoding RNAs (lincRNAs) in plant defence responses are yet to be fully explored.

RESULTS

In this study, we used strand-specific RNA sequencing to identify 1113 lincRNAs in potato (Solanum tuberosum) from stem tissues. The lincRNAs are expressed from all 12 potato chromosomes and generally smaller in size compared to protein-coding genes. Like in other plants, most potato lincRNAs possess single exons. A time-course RNA-seq analysis between a tolerant and a susceptible potato cultivar showed that 559 lincRNAs are responsive to Pectobacterium carotovorum subsp. brasiliense challenge compared to mock-inoculated controls. Moreover, coexpression analysis revealed that 17 of these lincRNAs are highly associated with 12 potato defence-related genes.

CONCLUSIONS

Together, these results suggest that lincRNAs have potential functional roles in potato defence responses. Furthermore, this work provides the first library of potato lincRNAs and a set of novel lincRNAs implicated in potato defences against P. carotovorum subsp. brasiliense, a member of the soft rot Enterobacteriaceae phytopathogens.

Genome Wide Identification and Functional Prediction of Long Non-Coding RNAs Responsive to Sclerotinia sclerotiorum Infection in Brassica napus

Sclerotinia stem rot caused by Sclerotinia sclerotiorum affects canola production worldwide. Emerging evidence suggests that long non-coding RNAs (lncRNAs) play important roles in the regulation of gene expression in plants, in response to both abiotic and biotic stress. So far, identification of lncRNAs has been limited to a few model plant species, and their roles in mediating responses to biotic stresses are yet to be characterized in Brassica napus. The present study reports the identification of novel lncRNAs responsive to S. sclerotiorum infection in B. napus at two time points after infection (24 hpi and 48 hpi) using a stranded RNA-Sequencing technique and a detection pipeline for lncRNAs. Of the total 3,181 lncRNA candidates, 2,821 lncRNAs were intergenic, 111 were natural antisense transcripts, 76 possessed exonic overlap with the reference coding transcripts while the remaining 173 represented novel lnc- isoforms. Forty one lncRNAs were identified as the precursors for microRNAs (miRNAs) including miR156, miR169 and miR394, with significant roles in mediating plant responses to fungal phytopathogens. A total of 931 differentially expressed lncRNAs were identified in response to S. sclerotiorum infection and the expression of 12 such lncRNAs was further validated using qRT-PCR. B. napus antisense lncRNA, TCONS_00000966, having 90% overlap with a plant defensin gene, showed significant induction at both infection stages, suggesting its involvement in the transcriptional regulation of defense responsive genes under S. sclerotiorum infection. Additionally, nine lncRNAs showed overlap with cis-regulatory regions of differentially expressed genes of B. napus. Quantitative RT-PCR verification of a set of S. sclerotiorum responsive sense/antisense transcript pairs revealed contrasting expression patterns, supporting the hypothesis that steric clashes of transcriptional machinery may lead to inactivation of sense promoter. Our findings highlight the potential contributions of lncRNAs in regulating expression of plant genes that respond to biotic stress.

Differential transcriptomic responses to Fusarium graminearum infection in two barley quantitative trait loci associated with Fusarium head blight resistance

BACKGROUND

Fusarium graminearum causes Fusarium head blight (FHB), a major disease problem worldwide. Resistance to FHB is controlled by quantitative trait loci (QTL) of which two are located on barley chromosomes 2H bin8 and 6H bin7. The mechanisms of resistance mediated by FHB QTL are poorly defined.

RESULTS

Near-isogenic lines (NILs) carrying Chevron-derived resistant alleles for the two QTL were developed and exhibited FHB resistance in field trials. To understand the molecular responses associated with resistance, transcriptomes of the NILs and recurrent parents (M69 and Lacey) were investigated with RNA sequencing (RNA-Seq) after F. graminearum or mock inoculation. A total of 2083 FHB-responsive transcripts were detected and provide a gene expression atlas for the barley-F. graminearum interaction. Comparative analysis of the 2Hb8 resistant (R) NIL and M69 revealed that the 2Hb8 R NIL exhibited an elevated defense response in the absence of fungal infection and responded quicker than M69 upon fungal infection. The 6Hb7 R NIL displayed a more rapid induction of a set of defense genes than Lacey during the early stage of fungal infection. Overlap of differentially accumulated genes were identified between the two R NILs, suggesting that certain responses may represent basal resistance to F. graminearum and/or general biotic stress response and were expressed by both resistant genotypes. Long noncoding RNAs (lncRNAs) have emerged as potential key regulators of transcription. A total of 12,366 lncRNAs were identified, of which 604 were FHB responsive.

CONCLUSIONS

The current transcriptomic analysis revealed differential responses conferred by two QTL during F. graminearum infection and identified genes and lncRNAs that were associated with FHB resistance.

Coexpression network analysis of the genes regulated by two types of resistance responses to powdery mildew in wheat

Powdery mildew disease caused by Blumeria graminis f. sp. tritici (Bgt) inflicts severe economic losses in wheat crops. A systematic understanding of the molecular mechanisms involved in wheat resistance to Bgt is essential for effectively controlling the disease. Here, using the diploid wheat Triticum urartu as a host, the genes regulated by immune (IM) and hypersensitive reaction (HR) resistance responses to Bgt were investigated through transcriptome sequencing. Four gene coexpression networks (GCNs) were developed using transcriptomic data generated for 20 T. urartu accessions showing IM, HR or susceptible responses. The powdery mildew resistance regulated (PMRR) genes whose expression was significantly correlated with Bgt resistance were identified, and they tended to be hubs and enriched in six major modules. A wide occurrence of negative regulation of PMRR genes was observed. Three new candidate immune receptor genes (TRIUR3_13045, TRIUR3_01037 and TRIUR3_06195) positively associated with Bgt resistance were discovered. Finally, the involvement of TRIUR3_01037 in Bgt resistance was tentatively verified through cosegregation analysis in a F2 population and functional expression assay in Bgt susceptible leaf cells. This research provides insights into the global network properties of PMRR genes. Potential molecular differences between IM and HR resistance responses to Bgt are discussed.

Regulation of miR394 in Response to Fusarium oxysporum f. sp. cepae (FOC) Infection in Garlic (Allium sativum L)

MicroRNAs (miRNAs) are a class of post-transcriptional regulators that negatively regulate gene expression through target mRNA cleavage or translational inhibition and play important roles in plant development and stress response. In the present study, six conserved miRNAs from garlic (Allium sativum L.) were analyzed to identify differentially expressed miRNAs in response to Fusarium oxysporum f. sp. cepae (FOC) infection. Stem-loop RT-PCR revealed that miR394 is significantly induced in garlic seedlings post-treatment with FOC for 72 h. The induction of miR394 expression during FOC infection was restricted to the basal stem plate tissue, the primary site of infection. Garlic miR394 was also upregulated by exogenous application of jasmonic acid. Two putative targets of miR394 encoding F-box domain and cytochrome P450 (CYP450) family proteins were predicted and verified using 5' RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends) assay. Quantitative RT-PCR showed that the transcript levels of the predicted targets were significantly reduced in garlic plants exposed to FOC. When garlic cultivars with variable sensitivity to FOC were exposed to the pathogen, an upregulation of miR394 and down regulation of the targets were observed in both varieties. However, the expression pattern was delayed in the resistant genotypes. These results suggest that miR394 functions in negative modulation of FOC resistance and the difference in timing and levels of expression in variable genotypes could be examined as markers for selection of FOC resistant garlic cultivars.

lncRNAs in Stress Response

All living organisms sense and respond to harmful changes in their intracellular and extracellular environment through complex signaling pathways that lead to changes in gene expression and cellular function in order to maintain homeostasis. Long non-coding RNAs (lncRNAs), a large and heterogeneous group of functional RNAs, play important roles in cellular response to stressful conditions. lncRNAs constitute a significant fraction of the genes differentially expressed in response to diverse stressful stimuli and, once induced, contribute to the regulation of downstream cellular processes, including feedback regulation of key stress response proteins. While many lncRNAs seem to be induced in response to a specific stress, there is significant overlap between lncRNAs induced in response to different stressful stimuli. In addition to stress-induced RNAs, several constitutively expressed lncRNAs also exert a strong regulatory impact on the stress response. Although our understanding of the contribution of lncRNAs to the cellular stress response is still highly rudimentary, the existing data point to the presence of a complex network of lncRNAs, miRNAs, and proteins in regulation of the cellular response to stress.

Genome-wide analysis of tomato long non-coding RNAs and identification as endogenous target mimic for microRNA in response to TYLCV infection

Recently, a large number of long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes in animals and plants. However, how lncRNAs function during plant DNA virus infection is largely unknown. We performed strand-specific paired-end RNA sequencing of tomato samples infected with Tomato yellow leaf curl virus (TYLCV) with three biological replicates. Overall, we predicted 1565 lncRNAs including long intergenic ncRNAs (lincRNAs) and natural antisense transcripts (lncNATs) and definitively identified lnRNAs that are involved in TYLCV infection by virus-induced gene silencing (VIGS). We also verified the functions of a set of lncRNAs that were differentially expressed between 0 and 7 days post inoculation (dpi). More importantly, we found that several lncRNAs acted as competing endogenous target mimics (eTMs) for tomato microRNAs involved in the TYLCV infection. These results provide new insight into lncRNAs involved in the response to TYLCV infection that are important components of the TYLCV network in tomatoes.

RNA sequencing and functional analysis implicate the regulatory role of long non-coding RNAs in tomato fruit ripening

Recently, long non-coding RNAs (lncRNAs) have been shown to play critical regulatory roles in model plants, such as Arabidopsis, rice, and maize. However, the presence of lncRNAs and how they function in fleshy fruit ripening are still largely unknown because fleshy fruit ripening is not present in the above model plants. Tomato is the model system for fruit ripening studies due to its dramatic ripening process. To investigate further the role of lncRNAs in fruit ripening, it is necessary and urgent to discover and identify novel lncRNAs and understand the function of lncRNAs in tomato fruit ripening. Here it is reported that 3679 lncRNAs were discovered from wild-type tomato and ripening mutant fruit. The lncRNAs are transcribed from all tomato chromosomes, 85.1% of which came from intergenic regions. Tomato lncRNAs are shorter and have fewer exons than protein-coding genes, a situation reminiscent of lncRNAs from other model plants. It was also observed that 490 lncRNAs were significantly up-regulated in ripening mutant fruits, and 187 lncRNAs were down-regulated, indicating that lncRNAs could be involved in the regulation of fruit ripening. In line with this, silencing of two novel tomato intergenic lncRNAs, lncRNA1459 and lncRNA1840, resulted in an obvious delay of ripening of wild-type fruit. Overall, the results indicated that lncRNAs might be essential regulators of tomato fruit ripening, which sheds new light on the regulation of fruit ripening.

Identification and characterization of long non-coding RNAs involved in osmotic and salt stress in Medicago truncatula using genome-wide high-throughput sequencing

BACKGROUND

Long non-coding RNAs (lncRNAs) have been shown to play crucially regulatory roles in diverse biological processes involving complex mechanisms. However, information regarding the number, sequences, characteristics and potential functions of lncRNAs in plants is so far overly limited.

RESULTS

Using high-throughput sequencing and bioinformatics analysis, we identified a total of 23,324 putative lncRNAs from control, osmotic stress- and salt stress-treated leaf and root samples of Medicago truncatula, a model legume species. Out of these lncRNAs, 7,863 and 5,561 lncRNAs were identified from osmotic stress-treated leaf and root samples, respectively. While, 7,361 and 7,874 lncRNAs were identified from salt stress-treated leaf and root samples, respectively. To reveal their potential functions, we analyzed Gene Ontology (GO) terms of genes that overlap with or are neighbors of the stress-responsive lncRNAs. Enrichments in GO terms in biological processes such as signal transduction, energy synthesis, molecule metabolism, detoxification, transcription and translation were found.

CONCLUSIONS

LncRNAs are likely involved in regulating plant's responses and adaptation to osmotic and salt stresses in complex regulatory networks with protein-coding genes. These findings are of importance for our understanding of the potential roles of lncRNAs in responses of plants in general and M. truncatula in particular to abiotic stresses.

Long non-coding RNAs and their biological roles in plants

With the development of genomics and bioinformatics, especially the extensive applications of high-throughput sequencing technology, more transcriptional units with little or no protein-coding potential have been discovered. Such RNA molecules are called non-protein-coding RNAs (npcRNAs or ncRNAs). Among them, long npcRNAs or ncRNAs (lnpcRNAs or lncRNAs) represent diverse classes of transcripts longer than 200 nucleotides. In recent years, the lncRNAs have been considered as important regulators in many essential biological processes. In plants, although a large number of lncRNA transcripts have been predicted and identified in few species, our current knowledge of their biological functions is still limited. Here, we have summarized recent studies on their identification, characteristics, classification, bioinformatics, resources, and current exploration of their biological functions in plants.

Exploring the secrets of long noncoding RNAs

High-throughput sequencing has revealed that the majority of RNAs have no capacity to encode protein. Among these non-coding transcripts, recent work has focused on the roles of long noncoding RNAs (lncRNAs) of >200 nucleotides. Although many of their attributes, such as patterns of expression, remain largely unknown, lncRNAs have key functions in transcriptional, post-transcriptional, and epigenetic gene regulation; Also, new work indicates their functions in scaffolding ribonuclear protein complexes. In plants, genome-wide identification of lncRNAs has been conducted in several species, including Zea mays, and recent research showed that lncRNAs regulate flowering time in the photoperiod pathway, and function in nodulation. In this review, we discuss the basic mechanisms by which lncRNAs regulate key cellular processes, using the large body of knowledge on animal and yeast lncRNAs to illustrate the significance of emerging work on lncRNAs in plants.

High-throughput sequencing and De Novo assembly of the Isatis indigotica transcriptome

Background

Isatis indigotica, the source of the traditional Chinese medicine Radix isatidis (Ban-Lan-Gen), is an extremely important economical crop in China. To facilitate biological, biochemical and molecular research on the medicinal chemicals in I. indigotica, here we report the first I. indigotica transcriptome generated by RNA sequencing (RNA-seq).

Results

RNA-seq library was created using RNA extracted from a mixed sample including leaf and root. A total of 33,238 unigenes were assembled from more than 28 million of high quality short reads. The quality of the assembly was experimentally examined by cDNA sequencing of seven randomly selected unigenes. Based on blast search 28,184 unigenes had a hit in at least one of the protein and nucleotide databases used in this study, and 8 unigenes were found to be associated with biosynthesis of indole and its derivatives. According to Gene Ontology classification, 22,365 unigenes were categorized into 48 functional groups. Furthermore, Clusters of Orthologous Group and Swiss-Port annotation were assigned for 7,707 and 18,679 unigenes, respectively. Analysis of repeat motifs identified 6,400 simple sequence repeat markers in 4,509 unigenes.

Conclusion

Our data provide a comprehensive sequence resource for molecular study of I. indigotica. Our results will facilitate studies on the functions of genes involved in the indole alkaloid biosynthesis pathway and on metabolism of nitrogen and indole alkaloids in I. indigotica and its related species.

DNA demethylases target promoter transposable elements to positively regulate stress responsive genes in Arabidopsis

Background

DNA demethylases regulate DNA methylation levels in eukaryotes. Arabidopsis encodes four DNA demethylases, DEMETER (DME), REPRESSOR OF SILENCING 1 (ROS1), DEMETER-LIKE 2 (DML2), and DML3. While DME is involved in maternal specific gene expression during seed development, the biological function of the remaining DNA demethylases remains unclear.

Results

We show that ROS1, DML2, and DML3 play a role in fungal disease resistance in Arabidopsis. A triple DNA demethylase mutant, rdd (ros1 dml2 dml3), shows increased susceptibility to the fungal pathogen Fusarium oxysporum. We identify 348 genes differentially expressed in rdd relative to wild type, and a significant proportion of these genes are downregulated in rdd and have functions in stress response, suggesting that DNA demethylases maintain or positively regulate the expression of stress response genes required for F. oxysporum resistance. The rdd-downregulated stress response genes are enriched for short transposable element sequences in their promoters. Many of these transposable elements and their surrounding sequences show localized DNA methylation changes in rdd, and a general reduction in CHH methylation, suggesting that RNA-directed DNA methylation (RdDM), responsible for CHH methylation, may participate in DNA demethylase-mediated regulation of stress response genes. Many of the rdd-downregulated stress response genes are downregulated in the RdDM mutants nrpd1 and nrpe1, and the RdDM mutants nrpe1 and ago4 show enhanced susceptibility to F. oxysporum infection.

Conclusions

Our results suggest that a primary function of DNA demethylases in plants is to regulate the expression of stress response genes by targeting promoter transposable element sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0458-3) contains supplementary material, which is available to authorized users.

Mutation of the RDR1 gene caused genome-wide changes in gene expression, regional variation in small RNA clusters and localized alteration in DNA methylation in rice

BACKGROUND

Endogenous small (sm) RNAs (primarily si- and miRNAs) are important trans/cis-acting regulators involved in diverse cellular functions. In plants, the RNA-dependent RNA polymerases (RDRs) are essential for smRNA biogenesis. It has been established that RDR2 is involved in the 24 nt siRNA-dependent RNA-directed DNA methylation (RdDM) pathway. Recent studies have suggested that RDR1 is involved in a second RdDM pathway that relies mostly on 21 nt smRNAs and functions to silence a subset of genomic loci that are usually refractory to the normal RdDM pathway in Arabidopsis. Whether and to what extent the homologs of RDR1 may have similar functions in other plants remained unknown.

RESULTS

We characterized a loss-of-function mutant (Osrdr1) of the OsRDR1 gene in rice (Oryza sativa L.) derived from a retrotransposon Tos17 insertion. Microarray analysis identified 1,175 differentially expressed genes (5.2% of all expressed genes in the shoot-tip tissue of rice) between Osrdr1 and WT, of which 896 and 279 genes were up- and down-regulated, respectively, in Osrdr1. smRNA sequencing revealed regional alterations in smRNA clusters across the rice genome. Some of the regions with altered smRNA clusters were associated with changes in DNA methylation. In addition, altered expression of several miRNAs was detected in Osrdr1, and at least some of which were associated with altered expression of predicted miRNA target genes. Despite these changes, no phenotypic difference was identified in Osrdr1 relative to WT under normal condition; however, ephemeral phenotypic fluctuations occurred under some abiotic stress conditions.

CONCLUSIONS

Our results showed that OsRDR1 plays a role in regulating a substantial number of endogenous genes with diverse functions in rice through smRNA-mediated pathways involving DNA methylation, and which participates in abiotic stress response.