Small RNAs and Gene Network in a Durable Disease Resistance Gene--Mediated Defense Responses in Rice

OsGELP77, a QTL for broad-spectrum disease resistance and yield in rice, encodes a GDSL-type lipase

Lipids and lipid metabolites have essential roles in plant-pathogen interactions. GDSL-type lipases are involved in lipid metabolism modulating lipid homeostasis. Some plant GDSLs modulate lipid metabolism altering hormone signal transduction to regulate host-defence immunity. Here, we functionally characterized a rice lipase, OsGELP77, promoting both immunity and yield. OsGELP77 expression was induced by pathogen infection and jasmonic acid (JA) treatment. Overexpression of OsGELP77 enhanced rice resistance to both bacterial and fungal pathogens, while loss-of-function of osgelp77 showed susceptibility. OsGELP77 localizes to endoplasmic reticulum and is a functional lipase hydrolysing universal lipid substrates. Lipidomics analyses demonstrate that OsGELP77 is crucial for lipid metabolism and lipid-derived JA homeostasis. Genetic analyses confirm that OsGELP77-modulated resistance depends on JA signal transduction. Moreover, population genetic analyses indicate that OsGELP77 expression level is positively correlated with rice resistance against pathogens. Three haplotypes were classified based on nucleotide polymorphisms in the OsGELP77 promoter where OsGELP77Hap3 is an elite haplotype. Three OsGELP77 haplotypes are differentially distributed in wild and cultivated rice, while OsGELP77Hap3 has been broadly pyramided for hybrid rice development. Furthermore, quantitative trait locus (QTL) mapping and resistance evaluation of the constructed near-isogenic line validated OsGELP77, a QTL for broad-spectrum disease resistance. In addition, OsGELP77-modulated lipid metabolism promotes JA accumulation facilitating grain yield. Notably, the hub defence regulator OsWRKY45 acts upstream of OsGELP77 by initiating the JA-dependent signalling to trigger immunity. Together, OsGELP77, a QTL contributing to immunity and yield, is a candidate for breeding broad-spectrum resistant and high-yielding rice.

Rice microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 modules regulate defenses against bacteria

Rice (Oryza sativa L.) microRNA156/529-SQUAMOSA PROMOTER BINDING PROTEIN-LIKE7/14/17 (miR156/529-SPL7/14/17) modules have pleiotropic effects on many biological pathways. OsSPL7/14 can interact with DELLA protein SLENDER RICE1 (SLR1) to modulate gibberellin acid (GA) signal transduction against the bacterial pathogen Xanthomonas oryzae pv. oryzae. However, whether the miR156/529-OsSPL7/14/17 modules also regulate resistance against other pathogens is unclear. Notably, OsSPL7/14/17 functioning as transcriptional activators, their target genes, and the corresponding downstream signaling pathways remain largely unexplored. Here, we demonstrate that miR156/529 play negative roles in plant immunity and that miR156/529-regulated OsSPL7/14/17 confer broad-spectrum resistance against 2 devastating bacterial pathogens. Three OsSPL7/14/17 proteins directly bind to the promoters of rice Allene Oxide Synthase 2 (OsAOS2) and NONEXPRESSOR OF PATHOGENESIS-RELATED GENES1 (OsNPR1) and activate their transcription, regulating jasmonic acid (JA) accumulation and the salicylic acid (SA) signaling pathway, respectively. Overexpression of OsAOS2 or OsNPR1 impairs the susceptibility of the osspl7/14/17 triple mutant. Exogenous application of JA enhances resistance of the osspl7/14/17 triple mutant and the miR156 overexpressing plants. In addition, genetic evidence confirms that bacterial pathogen-activated miR156/529 negatively regulate pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) responses, such as pattern recognition receptor Xa3/Xa26-initiated PTI. Our findings demonstrate that bacterial pathogens modulate miR156/529-OsSPL7/14/17 modules to suppress OsAOS2-catalyzed JA accumulation and the OsNPR1-promoted SA signaling pathway, facilitating pathogen infection. The uncovered miR156/529-OsSPL7/14/17-OsAOS2/OsNPR1 regulatory network provides a potential strategy to genetically improve rice disease resistance.

Cadmium accumulation regulated by a rice heavy-metal importer is harmful for host plant and leaf bacteria

INTRODUCTION

Cadmium (Cd), one of the major toxic heavy metals, causes severe deleterious effects on all living organisms from prokaryotes to eukaryotes. Cadmium deposition affects bacterial diversity and bacterial population in soil. Cadmium accumulation in plants is mainly controlled by transporters and the resulting Cd enrichment gives rise to phytotoxicity.

OBJECTIVE

This study aimed to mine transporters that control Cd import or accumulation in rice and uncover the underlying mechanisms that how accumulated Cd poses risks to host plant and leaf bacteria.

METHODS

RNA-seq analysis, histochemical assays, and elemental quantification were carried out to reveal the biological roles of OsABCG43 for Cd import. Pathogen inoculation, IC₅₀ value, and bacterial virulence assays were conducted to disclose the effects of Cd on leaf bacteria.

RESULTS

OsABCG43 is characterized as a Cd importer controlling Cd accumulation in rice. OsABCG43 was induced under Cd stress and specifically expressed in the vasculature of leaves and roots. Overexpression of OsABCG43 caused Cd accumulation which inhibits photosynthesis and development and alters the antioxidant system, resulting in phytotoxicity. Moreover, overexpression of OsABCG43 resulted in retarded plant growth and enhanced rice sensitivity to Cd stress. Numerous differentially expressed genes were identified via RNA-seq analysis between the OsABCG43-overexpressing plants and wild type, which functioned in Cd or reactive oxygen species (ROS) homeostasis. In addition, OsABCG43 transcripts were induced by leaf bacteria Xanthomonas oryzae pv. oryzicola (Xoc) and X. oryzae pv. oryzae (Xoo). The enriched Cd directly impaired the formation of virulence factors for the leaf bacteria, preventing colonization or proliferation of Xoc or Xoo in rice leaves.

CONCLUSION

This work reveals that OsABCG43 is expressed specifically in the vascular and plasma membrane-localized OsABCG43 functions as a Cd importer. OsABCG43-mediated import of Cd is harmful for both rice and the corresponding leaf bacteria.

miR395-regulated sulfate metabolism exploits pathogen sensitivity to sulfate to boost immunity in rice

MicroRNAs (miRNAs) play important roles in plant physiological activities. However, their roles and molecular mechanisms in boosting plant immunity, especially through the modulation of macronutrient metabolism in response to pathogens, are largely unknown. Here, we report that an evolutionarily conserved miRNA, miR395, promotes resistance to Xanthomonas oryzae pv. oryzae (Xoo) and X. oryzae pv. oryzicola (Xoc), two destructive bacterial pathogens, by regulating sulfate accumulation and distribution in rice. Specifically, miR395 targets and suppresses the expression of the ATP sulfurylase gene OsAPS1, which functions in sulfate assimilation, and two sulfate transporter genes, OsSULTR2;1 and OsSULTR2;2, which function in sulfate translocation, to promote sulfate accumulation, resulting in broad-spectrum resistance to bacterial pathogens in miR395-overexpressing plants. Genetic analysis revealed that miR395-triggered resistance is involved in both pathogen-associated molecular pattern-triggered immunity and R gene-mediated resistance. Moreover, we found that accumulated sulfate but not S-metabolites inhibits proliferation of pathogenic bacteria, revealing a sulfate-mediated antibacterial defense mechanism that differs from sulfur-induced resistance. Furthermore, compared with other bacteria, Xoo and Xoc, which lack the sulfate transporter CysZ, are sensitive to high levels of extracellular sulfate. Accordingly, miR395-regulated sulfate accumulation impaired the virulence of Xoo and Xoc by decreasing extracellular polysaccharide production and biofilm formation. Taken together, these results suggest that rice miR395 modulates sulfate metabolism to exploit pathogen sensitivity to sulfate and thereby promotes broad-spectrum resistance.

The rice heavy-metal transporter OsNRAMP1 regulates disease resistance by modulating ROS homoeostasis

Crop diseases threaten food security and sustainable agriculture. Consumption of crops containing nonessential toxic metals leads to health risks for humans. Therefore, cultivation of disease-resistant and toxic metal-safe crops is a double-gain strategy that can contribute to food security. Here, we show that rice heavy-metal transporter OsNRAMP1 plays an important role in plant immunity by modulating metal ion and reactive oxygen species (ROS) homoeostasis. OsNRAMP1 expression was induced after pathogenic bacteria and fungi infections. The osnramp1 mutants had an increased content of H₂ O₂ and activity of superoxide dismutase, but decreased activity of catalase, showing enhanced broad-spectrum resistance against bacterial and fungal pathogens. RNA-seq analysis identified a number of differentially expressed genes that were involved in metal ion and ROS homoeostasis. Altered expression of metal ion-dependent ROS-scavenging enzymes genes and lower accumulation of cations such as Mn together induced compromised metal ion-dependent enzyme-catalysing activity and modulated ROS homoeostasis, which together contributed towards disease resistance in osnramp1 mutants. Furthermore, the osnramp1 mutants contained lower levels of toxic heavy metals Cd and Pb and micronutrients Ni and Mn in leaves and grains. Taken together, a proof of concept was achieved that broad-spectrum disease-resistant and toxic heavy-metal-safe rice was engineered by removal of the OsNRAMP1 gene.

Epigenetic regulation of salinity stress responses in cereals

Cereals are important crops and are exposed to various types of environmental stresses that affect the overall growth and yield. Among the various abiotic stresses, salt stress is a major environmental factor that influences the genetic, physiological, and biochemical responses of cereal crops. Epigenetic regulation which includes DNA methylation, histone modification, and chromatin remodelling plays an important role in salt stress tolerance. Recent studies in rice genomics have highlighted that the epigenetic changes are heritable and therefore can be considered as molecular signatures. An epigenetic mechanism under salinity induces phenotypic responses involving modulations in gene expression. Association between histone modification and altered DNA methylation patterns and differential gene expression has been evidenced for salt sensitivity in rice and other cereal crops. In addition, epigenetics also creates stress memory that helps the plant to better combat future stress exposure. In the present review, we have discussed epigenetic influences in stress tolerance, adaptation, and evolution processes. Understanding the epigenetic regulation of salinity could help for designing salt-tolerant varieties leading to improved crop productivity.

Inducible Enrichment of Osa-miR1432 Confers Rice Bacterial Blight Resistance through Suppressing OsCaML2

MicroRNAs (miRNAs) handle immune response to pathogens by adjusting the function of target genes in plants. However, the experimentally documented miRNA/target modules implicated in the interplay between rice and Xanthomonas oryzae pv. oryzae (Xoo) are still in the early stages. Herein, the expression of osa-miR1432 was induced in resistant genotype IRBB5, but not susceptible genotype IR24, under Xoo strain PXO86 attack. Overexpressed osa-miR1432 heightened rice disease resistance to Xoo, indicated by enhancive enrichment of defense marker genes, raised reactive oxygen species (ROS) levels, repressed bacterial growth and shortened leaf lesion length, whilst the disruptive accumulation of osa-miR1432 accelerated rice susceptibility to Xoo infection. Noticeably, OsCaML2 (LOC_Os03g59770) was experimentally confirmed as a target gene of osa-miR1432, and the overexpressing OsCaML2 transgenic plants exhibited compromised resistance to Xoo infestation. Our results indicate that osa-miR1432 and OsCaML2 were differently responsive to Xoo invasion at the transcriptional level and fine-tune rice resistance to Xoo infection, which may be referable in resistance gene discovery and valuable in the pursuit of improving Xoo resistance in rice breeding.

Identification of miRNAs Contributing to the Broad-Spectrum and Durable Blast Resistance in the Yunnan Local Rice Germplasm

MicroRNAs are 20-24 nucleotide non-coding RNAs and play important roles in plant-environment interactions. In recent years, many microRNAs (miRNAs) have been found to regulate rice immunity against rice blast fungus. However, there are limited studies about miRNAs that directly target resistance (R) genes to regulate rice immunity. In this study, by deep sequencing, small RNA libraries were constructed from four-leaf stage seedlings of the resistant variety Ziyu44 and susceptible variety Jiangnanxiangnuo (JNXN) upon Magnaporthe oryzae infection, we found that much more miRNAs were significantly differentially expressed in Ziyu44 than in JNXN. Among these miRNAs, we focused on miR9664, a newly identified rice miRNA in our sequencing, which was upregulated lightly in Ziyu44 and drastically in JNXN at 24-48 h post-inoculation (hpi). The transgenic plants overexpressing miR9664 (miR9664-oe) displayed reduced defense responses to M. oryzae, while those knocking down miR9664 (miR9664-m) displayed enhanced defense responses to M. oryzae. Most of the detected miR9664 predicted target genes were reduced in the miR9664-oe lines while increased in the miR9664-m lines. The cleavage site of LOC_Os08g07774 was confirmed by RLM-RACE. Meanwhile, after being inoculated with M. oryzae, the genes were expressed differently between Ziyu44 and JNXN. The results suggest that miR9664-mediated R gene turnover contributes to Ziyu44 broad-spectrum resistance to rice blast fungus. Taken together, our research identified a new rice miRNA that directly targets R genes to regulate rice immunity against rice blast fungus, adding significant information to the study of rice-M. oryzae interaction.

Genome-Wide Analysis of Potassium Channel Genes in Rice: Expression of the OsAKT and OsKAT Genes under Salt Stress

Potassium (K+), as a vital element, is involved in regulating important cellular processes such as enzyme activity, cell turgor, and nutrient movement in plant cells, which affects plant growth and production. Potassium channels are involved in the transport and release of potassium in plant cells. In the current study, three OsKAT genes and two OsAKT genes, along with 11 nonredundant putative potassium channel genes in the rice genome, were characterized based on their physiochemical properties, protein structure, evolution, duplication, in silico gene expression, and protein-protein interactions. In addition, the expression patterns of OsAKTs and OsKATs were studied in root and shoot tissues under salt stress using real-time PCR in three rice cultivars. K+ channel genes were found to have diverse functions and structures, and OsKATs showed high genetic divergence from other K+ channel genes. Furthermore, the Ka/Ks ratios of duplicated gene pairs from the K+ channel gene family in rice suggested that these genes underwent purifying selection. Among the studied K+ channel proteins, OsKAT1 and OsAKT1 were identified as proteins with high potential N-glycosylation and phosphorylation sites, and LEU, VAL, SER, PRO, HIS, GLY, LYS, TYR, CYC, and ARG amino acids were predicted as the binding residues in the ligand-binding sites of K+ channel proteins. Regarding the coexpression network and KEGG ontology results, several metabolic pathways, including sugar metabolism, purine metabolism, carbon metabolism, glycerophospholipid metabolism, monoterpenoid biosynthesis, and folate biosynthesis, were recognized in the coexpression network of K+ channel proteins. Based on the available RNA-seq data, the K+ channel genes showed differential expression levels in rice tissues in response to biotic and abiotic stresses. In addition, the real-time PCR results revealed that OsAKTs and OsKATs are induced by salt stress in root and shoot tissues of rice cultivars, and OsKAT1 was identified as a key gene involved in the rice response to salt stress. In the present study, we found that the repression of OsAKTs, OsKAT2, and OsKAT2 in roots was related to salinity tolerance in rice. Our findings provide valuable insights for further structural and functional assays of K+ channel genes in rice.

Contribution of Small RNA Pathway to Interactions of Rice with Pathogens and Insect Pests

Small RNAs (sRNAs) are mainly classified into microRNAs (miRNAs) and small interfering RNAs (siRNAs) according to their origin. miRNAs originate from single-stranded RNA precursors, whereas siRNAs originate from double-stranded RNA precursors that are synthesized by RNA-dependent RNA polymerases. Both of single-stranded and double-stranded RNA precursors are processed into sRNAs by Dicer-like proteins. Then, the sRNAs are loaded into ARGONAUTE proteins, forming RNA-induced silencing complexes (RISCs). The RISCs repress the expression of target genes with sequences complementary to the sRNAs through the cleavage of transcripts, the inhibition of translation or DNA methylation. Here, we summarize the recent progress of sRNA pathway in the interactions of rice with various parasitic organisms, including fungi, viruses, bacteria, as well as insects. Besides, we also discuss the hormone signal in sRNA pathway, and the emerging roles of circular RNAs and long non-coding RNAs in rice immunity. Obviously, small RNA pathway may act as a part of rice innate immunity to coordinate with growth and development.

Understanding epigenomics based on the rice model

The purpose of this paper provides a comprehensive overview of the recent researches on rice epigenomics, including DNA methylation, histone modifications, noncoding RNAs, and three-dimensional genomics. The challenges and perspectives for future research in rice are discussed. Rice as a model plant for epigenomic studies has much progressed current understanding of epigenetics in plants. Recent results on rice epigenome profiling and three-dimensional chromatin structure studies reveal specific features and implication in gene regulation during rice plant development and adaptation to environmental changes. Results on rice chromatin regulator functions shed light on mechanisms of establishment, recognition, and resetting of epigenomic information in plants. Cloning of several rice epialleles associated with important agronomic traits highlights importance of epigenomic variation in rice plant growth, fitness, and yield. In this review, we summarize and analyze recent advances in rice epigenomics and discuss challenges and directions for future research in the field.

Transcriptional landscape of pathogen-responsive lncRNAs in rice unveils the role of ALEX1 in jasmonate pathway and disease resistance

Plant defence is multilayered and is essential for surviving in a changing environment. The discovery of long noncoding RNAs (lncRNAs) has dramatically extended our understanding of post-transcriptional gene regulation in diverse biological processes. However, the expression profile and function of lncRNAs in disease resistance are still largely unknown, especially in monocots. Here, we performed strand-specific RNA sequencing of rice leaves infected by Xanthomonas oryzae pv. Oryzae (Xoo) in different time courses and systematically identified 567 disease-responsive rice lncRNAs. Target analyses of these lncRNAs showed that jasmonate (JA) pathway was significantly enriched. To reveal the interaction between lncRNAs and JA-related genes, we studied the coexpression of them and found 39 JA-related protein-coding genes to be interplayed with 73 lncRNAs, highlighting the potential modulation of lncRNAs in JA pathway. We subsequently identified an lncRNA, ALEX1, whose expression is highly induced by Xoo infection. A T-DNA insertion line constructed using enhancer trap system showed a higher expression of ALEX1 and exerted a significant resistance to rice bacterial blight. Functional study revealed that JA signalling is activated and the endogenous content of JA and JA-Ile is increased. Overexpressing ALEX1 in rice further confirmed the activation of JA pathway and resistance to bacterial blight. Our findings reveal the expression of pathogen-responsive lncRNAs in rice and provide novel insights into the connection between lncRNAs and JA pathway in the regulation of plant disease resistance.

Characteristic Dissection of Xanthomonas oryzae pv. oryzae Responsive MicroRNAs in Rice

MicroRNAs (miRNAs) are crucial player in plant-pathogen interaction. While the evidence has demonstrated that rice miRNAs mediate immune response to pathogens invasion, the roles of miRNAs on Xanthomonas oryzae pv. oryzae (Xoo) attack remain be in place. Herein, we monitored the responsive changes of rice miRNAs at 0, 8, 24 h across Xoo strain PXO86 infection in its compatible rice variety IR24 and incompatible variety IRBB5 by small RNA sequencing, and the genes targeted by miRNAs were also detected via degradome technology. The faithfulness of sequencing data was validated through quantitative real-time stem-loop reverse transcription-polymerase chain reaction assay. Bioinformatic analysis showed that the differentially expressed miRNAs could be divided into three immunity-related clusters, and 80 regulatory units were emerged in infection process, which comprises 29 differentially expressed known miRNAs and 38 cleaved targets. Furthermore, the miRNA presumptive function of separate immunity cluster in rice-Xoo interplay was confirmed through overexpressing osa-miR164a, osa-miR167d and osa-miR159b, and the disruption of regulatory units, osa-miR164a/OsNAC60, osa-miR167d-5p/OsWD40-174 and osa-miR159b/OsMYBGA, OsLRR-RLK2, OsMPK20-4, may reset rice defense response to Xoo infestation in a controllable manner. These findings provide new insights into the complex roles of characteristic miRNAs and their targets in rice-Xoo interactions.

Integration Analysis of Small RNA and Degradome Sequencing Reveals MicroRNAs Responsive to Dickeya zeae in Resistant Rice

Rice foot rot disease caused by the pathogen Dickeya zeae (formerly known as Erwinia chrysanthemi pv. zeae), is a newly emerging damaging bacterial disease in China and the southeast of Asia, resulting in the loss of yield and grain quality. However, the genetic resistance mechanisms mediated by miRNAs to D. zeae are unclear in rice. In the present study, 652 miRNAs including osa-miR396f predicted to be involved in multiple defense responses to D. zeae were identified with RNA sequencing. A total of 79 differentially expressed miRNAs were detected under the criterion of normalized reads ≥10, including 51 known and 28 novel miRNAs. Degradome sequencing identified 799 targets predicted to be cleaved by 168 identified miRNAs. Among them, 29 differentially expressed miRNA and target pairs including miRNA396f-OsGRFs were identified by co-expression analysis. Overexpression of the osa-miR396f precursor in a susceptible rice variety showed enhanced resistance to D. zeae, coupled with significant accumulation of transcripts of osa-miR396f and reduction of its target the Growth-Regulating Factors (OsGRFs). Taken together, these findings suggest that miRNA and targets including miR396f-OsGRFs have a role in resisting the infections by bacteria D. zeae.

Genome-wide analysis of the rice PPR gene family and their expression profiles under different stress treatments

Background

Pentatricopeptide-repeat proteins (PPRs) are characterized by tandem arrays of a degenerate 35-amino-acid (PPR motifs), which can bind RNA strands and participate in post-transcription. PPR proteins family is one of the largest families in land plants and play important roles in organelle RNA metabolism and plant development. However, the functions of PPR genes involved in biotic and abiotic stresses of rice (Oryza sativa L.) remain largely unknown.

Results

In the present study, a comprehensive genome-wide analysis of PPR genes was performed. A total of 491 PPR genes were found in the rice genome, of which 246 PPR genes belong to the P subfamily, and 245 genes belong to the PLS subfamily. Gene structure analysis showed that most PPR genes lack intron. Chromosomal location analysis indicated that PPR genes were widely distributed in all 12 rice chromosomes. Phylogenetic relationship analysis revealed the distinct difference between the P and PLS subfamilies. Many PPR proteins are predicted to target chloroplasts or mitochondria, and a PPR protein (LOC_Os10g34310) was verified to localize in mitochondria. Furthermore, three PPR genes (LOC_Os03g17634,LOC_Os07g40820,LOC_Os04g51350) were verified as corresponding miRNA targets. The expression pattern analysis showed that many PPR genes could be induced under biotic and abiotic stresses. Finally, seven PPR genes were confirmed with their expression patterns under salinity or drought stress.

Conclusions

We found 491 PPR genes in the rice genome, and our genes structure analysis and syntenic analysis indicated that PPR genes might be derived from amplification by retro-transposition. The expression pattern present here suggested that PPR proteins have crucial roles in response to different abiotic stresses in rice. Taken together, our study provides a comprehensive analysis of the PPR gene family and will facilitate further studies on their roles in rice growth and development.

Electronic supplementary material

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

Identification of conserved and novel microRNAs in Porphyridium purpureum via deep sequencing and bioinformatics

Background

Porphyridium purpureum has been utilized in important industrial and pharmaceutical fields. The identification of microRNAs (miRNAs) in this unique species is of great importance: such identification can help fill gaps in the small RNA (sRNA) studies of this organism and help to elucidate essential biological processes and their regulation mechanisms in this special micro alga.

Results

In this study, 254 high-confidence miRNAs (203 conserved miRNAs and 51 novel miRNAs) were identified by sRNA deep sequencing (sRNA-seq) combined with bioinformatics. A total of 235 putative miRNA families were predicted, including 192 conserved families and 43 species-specific families. The conservation and diversity of predicted miRNA families were analysed in different plant species. Both the 100 % northern blot validation rate (VR) of four randomly selected miRNAs and the results of stem-loop quantitative real time RT-PCR (qRT-PCR) assays of 25 randomly selected miRNAs demonstrated that the majority of the miRNAs identified in this study are credible. A total of 14,958 and 2184 genes were predicted to be targeted by the 186 conserved and 41 novel miRNAs. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that some target genes likely provide valuable references for further understanding of vital functions in P. purpureum. In addition, a cytoscape network will provide some clues for research into the complex biological processes that occur in this unique alga.

Conclusions

We first identified a large set of conserved and novel miRNAs in P. purpureum. The characteristic and validation analysis on miRNAs demonstrated authenticity of identification data. Functional annotation of target genes and metabolic pathways they involved in illuminated the direction for further utilization and development this micro alga based on its unique properties.

Electronic supplementary material

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