Construction of small RNA-mediated gene regulatory networks in the roots of rice (Oryza sativa)

A miRNA-Encoded Small Peptide, vvi-miPEP171d1, Regulates Adventitious Root Formation

One of the biggest challenges in clonal propagation of grapevine (Vitis vinifera) is difficulty of rooting. Adventitious root initiation and development are the critical steps in the cutting and layering process of grapevine, but the molecular mechanism of these processes remains unclear. Previous reports have found that microRNA (miRNA)-encoded peptides (miPEPs) can regulate plant root development by increasing the transcription of their corresponding primary miRNA. Here, we report the role of a miPEP in increasing adventitious root formation in grapevine. In this study, we performed a global analysis of miPEPs in grapevine and characterized the function of vvi-miPEP171d1, a functional, small peptide encoded by primary-miR171d. There were three small open reading frames in the 500-bp upstream sequence of pre-miR171d. One of them encoded a small peptide, vvi-miPEP171d1, which could increase the transcription of vvi-MIR171d Exogenous application of vvi-miPEP171d1 to grape tissue culture plantlets promoted adventitious root development by activating the expression of vvi-MIR171d Interestingly, neither exogenous application of the vvi-miPEP171d1 peptide nor overexpression of the vvi-miPEP171d1 coding sequence resulted in phenotypic changes in Arabidopsis (Arabidopsis thaliana). Similarly, application of synthetic ath-miPEP171c, the small peptide encoded by the Arabidopsis ortholog of vvi-MIR171d, inhibited the growth of primary roots and induced the early initiation of lateral and adventitious roots in Arabidopsis, while it had no effect on grape root development. Our findings reveal that miPEP171d1 regulates root development by promoting vvi-MIR171d expression in a species-specific manner, further enriching the theoretical research into miPEPs.

Mapping the Salt Stress-Induced Changes in the Root miRNome in Pokkali Rice

A plant's response to stress conditions is governed by intricately coordinated gene expression. The microRNAs (miRs) have emerged as relatively new players in the genetic network, regulating gene expression at the transcriptional and post-transcriptional level. In this study, we performed comprehensive profiling of miRs in roots of the naturally salt-tolerant Pokkali rice variety to understand their role in regulating plant physiology in the presence of salt. For comparisons, root miR profiles of the salt-sensitive rice variety Pusa Basmati were generated. It was seen that the expression levels of 65 miRs were similar for roots of Pokkali grown in the absence of salt (PKNR) and Pusa Basmati grown in the presence of salt (PBSR). The salt-induced dis-regulations in expression profiles of miRs showed controlled changes in the roots of Pokkali (PKSR) as compared to larger variations seen in the roots of Pusa Basmati. Target analysis of salt-deregulated miRs identified key transcription factors, ion-transporters, and signaling molecules that act to maintain cellular Ca2+ homeostasis and limit ROS production. These miR:mRNA nodes were mapped to the Quantitative trait loci (QTLs) to identify the correlated root traits for understanding their significance in plant physiology. The results obtained indicate that the adaptability of Pokkali to excess salt may be due to the genetic regulation of different cellular components by a variety of miRs.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Identification of browning-related microRNAs and their targets reveals complex miRNA-mediated browning regulatory networks in Luffa cylindrica

As a non-coding and endogenous small RNA, MicroRNA (miRNA) takes a vital regulatory role in plant growth and development. Long-term storage and processing of many fruits and vegetables, including Luffa, are subject to influences from browning, a common post-harvest problem that adversely affects flavor, aroma, and nutritional value. The browning regulatory networks mediated by miRNA, however, remain largely unexplored. For a systematic identification of browning-responsive miRNAs and the targets, we built two RNA libraries from Luffa pulps of near-isogenic line, with resistant and sensitive browning characteristics respectively, and then sequenced them using Solexa high-throughput technology. We consequently identified 179 known miRNAs that represent 17 non-conserved miRNA families and 24 conserved families, as well as 84 potential novel miRNAs, among which 16 miRNAs (eight known and eight novel miRNAs) were found to exhibit significant differential expressions and were thus identified as browning-related miRNAs. We then studied those browning-responsive miRNAs and the corresponding targets with RT-qPCR and finally validated their expression patterns. The results revealed that the expression patterns are specific to plant development stages and the miRNAs are identified with 39 target transcripts, which involve in plant development, defense response, transcriptional regulation, and signal transduction. After characterizing these miRNAs and their targets, we propose a browning regulatory network model of miRNA-mediatation in this paper. The findings of the work are helpful for the understanding of miRNA-mediated regulatory mechanisms of browning in Luffa, and will facilitate genetic improvement of pulp characteristics in Luffa.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

Genome-wide changes in microRNA expression during short and prolonged heat stress and recovery in contrasting rice cultivars

MicroRNAs (miRNAs) are known to regulate expression of genes under stress. We report here the deep sequencing of small RNAs expressed during control, short and prolonged heat stress and recovery. Genome-wide identification of miRNAs in tolerant (Nagina 22) and susceptible (Vandana) rice cultivars was performed in 16 samples representing root and shoot of 13-day-old seedlings. The expression profile of miRNAs was analysed in 36 pairwise combinations to identify the genotype-, treatment- and tissue-dependent expression of miRNAs. Small-RNA sequencing of 16 libraries yielded ~271 million high-quality raw sequences; 162 miRNA families were identified. The highly expressed miRNAs in rice tissues were miR166, miR168, miR1425, miR529, mR162, miR1876, and miR1862. Expression of osa-miR1436, osa-miR5076, osa-miR5161, and osa-miR6253 was observed only in stressed tissue of both genotypes indicating their general role in heat stress response. Expression of osa-miR1439, osa-miR1848, osa-miR2096, osa-miR2106, osa-miR2875, osa-miR3981, osa-miR5079, osa-miR5151, osa-miR5484, osa-miR5792, and osa-miR5812 was observed only in Nagina 22 during high temperature, suggesting a specific role of these miRNAs in heat stress tolerance. This study provides details of the repertoire of miRNAs expressed in root and shoot of heat susceptible and tolerant rice genotypes under heat stress and recovery.

Identification of microRNAs in Response to Drought in Common Wild Rice (Oryza rufipogon Griff.) Shoots and Roots

Background

Drought is the most important factor that limits rice production in drought-prone environments. Plant microRNAs (miRNAs) are involved in biotic and abiotic stress responses. Common wild rice (Oryza rufipogon Griff.) contains abundant drought-resistant genes, which provide an opportunity to explore these excellent resources as contributors to improve rice resistance, productivity, and quality.

Results

In this study, we constructed four small RNA libraries, called CL and CR from PEG6000-free samples and DL and DR from PEG6000-treated samples, where ‘R’ indicates the root tissue and ‘L’ indicates the shoot tissue. A total of 200 miRNAs were identified to be differentially expressed under the drought-treated conditions (16% PEG6000 for 24 h), and the changes in the miRNA expression profile of the shoot were distinct from those of the root. At the miRNA level, 77 known miRNAs, which belong to 23 families, including 40 up-regulated and 37 down-regulated in the shoot, and 85 known miRNAs in 46 families, including 65 up-regulated and 20 down-regulated in the root, were identified as differentially expressed. In addition, we predicted 26 new miRNA candidates from the shoot and 43 from the root that were differentially expressed during the drought stress. The quantitative real-time PCR analysis results were consistent with high-throughput sequencing data. Moreover, 88 miRNAs that were differentially-expressed were predicted to match with 197 targets for drought-stress.

Conclusion

Our results suggest that the miRNAs of O. rufipogon are responsive to drought stress. The differentially expressed miRNAs that are tissue-specific under drought conditions could play different roles in the regulation of the auxin pathway, the flowering pathway, the drought pathway, and lateral root formation. Thus, the present study provides an account of tissue-specific miRNAs that are involved in the drought adaption of O. rufipogon.

A network analysis of miRNA mediated gene regulation of rice: crosstalk among biological processes

To understand the network architecture of miRNA mediated regulations at the genomic and functional levels of rice, we have made an unambiguous annotation of the experimentally verified miRNAs, predicted their targets and the possible biological functions they can affect. Some functions, namely translational and protein modifications and photosynthesis are targeted by higher percentage of miRNA. Using transformation procedures, we constructed a genome scale miRNA-miRNA functional synergistic network (MFSN). The analysis of MFSN modules help to identify miRNAs co-regulating target genes having several interrelated biological processes. Some of these target genes are also co-expressed under particular conditions. For example, the genes co-expressed under drought conditions as well as those targeted by miRNAs present in a MFSN module have interdependent biological processes namely, photosynthesis, cell-wall biogenesis, root development and xylan synthesis. The stress-induced miRNAs and their distributions, and the presence of transcription factors in the target set of MFSN modules were also analyzed.

Identification of microRNAs and Their Target Genes Explores miRNA-Mediated Regulatory Network of Cytoplasmic Male Sterility Occurrence during Anther Development in Radish (Raphanus sativus L.)

MicroRNAs (miRNAs) are a type of endogenous non-coding small RNAs that play critical roles in plant growth and developmental processes. Cytoplasmic male sterility (CMS) is typically a maternally inherited trait and widely used in plant heterosis utilization. However, the miRNA-mediated regulatory network of CMS occurrence during anther development remains largely unknown in radish. In this study, a comparative small RNAome sequencing was conducted in floral buds of CMS line 'WA' and its maintainer line 'WB' by high-throughput sequencing. A total of 162 known miRNAs belonging to 25 conserved and 24 non-conserved miRNA families were isolated and 27 potential novel miRNA families were identified for the first time in floral buds of radish. Of these miRNAs, 28 known and 14 potential novel miRNAs were differentially expressed during anther development. Several target genes for CMS occurrence-related miRNAs encode important transcription factors and functional proteins, which might be involved in multiple biological processes including auxin signaling pathways, signal transduction, miRNA target silencing, floral organ development, and organellar gene expression. Moreover, the expression patterns of several CMS occurrence-related miRNAs and their targets during three stages of anther development were validated by qRT-PCR. In addition, a potential miRNA-mediated regulatory network of CMS occurrence during anther development was firstly proposed in radish. These findings could contribute new insights into complex miRNA-mediated genetic regulatory network of CMS occurrence and advance our understanding of the roles of miRNAs during CMS occurrence and microspore formation in radish and other crops.

Are ta-siRNAs only originated from the cleavage site of miRNA on its target RNAs and phased in 21-nt increments?

Trans-acting siRNAs (ta-siRNAs) are a class of small RNAs playing crucial roles in the regulation of plant gene expression. According to the canonical model, specific miRNA-guided cleavage of a TAS transcript triggers and sets the registry for the subsequent production of ta-siRNAs at 21-nt increments from the cleavage site. However, a previously validated 22-nt ta-siR2140 indicated that ta-siRNAs might be initiated from other phase increments and registers, which resulted in massive ta-siRNAs missing in the canonical model. To test this hypothesis, we employed high-throughput sequencing data to thoroughly identify the miR173-triggered ta-siRNAs from TAS1/TAS2 transcripts. As a result, thousands of phased siRNAs not generated through the canonical pathway were identified and 110 novel siRNA-target interactions were further validated based on degradome sequencing data. Based on these results, we propose that the canonical biogenesis model of ta-siRNAs should be modified in order to recruit the previously unidentified ta-siRNA candidates.

Genome-wide identification of microRNAs associated with taproot development in radish (Raphanus sativus L.)

MicroRNAs (miRNAs) are small, endogenous, non-coding RNAs that play vital regulatory roles in plant growth and development. To identify the miRNAs associated with taproot development at the whole genome level, we sequenced five RNA libraries constructed from radish taproots at different developmental stages and generated a total of 148M clean reads. Using an integrative bioinformatics analysis, 494 known miRNAs belonging to 434 families and 220 putative novel miRNAs were identified. Combining the differential expression analysis and target prediction, we found that 77 miRNAs were potentially associated with taproot development. Target transcripts generated significant GO terms relating to cell proliferation, root development and hormone-mediated signaling. The KEGG analyses revealed that plant hormone signal transduction, zeatin biosynthesis, biosynthesis of secondary metabolites, cell cycle, MAPK signaling and p53 signaling were closely associated with taproot development. These findings will provide valuable information for further functional verification of miRNAs and their targets in radish taproot development.

Transcriptome profiling of root microRNAs reveals novel insights into taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Radish (Raphanus sativus L.) is an economically important root vegetable crop, and the taproot-thickening process is the most critical period for the final productivity and quality formation. MicroRNAs (miRNAs) are a family of non-coding small RNAs that play an important regulatory function in plant growth and development. However, the characterization of miRNAs and their roles in regulating radish taproot growth and thickening remain largely unexplored. A Solexa high-throughput sequencing technology was used to identify key miRNAs involved in taproot thickening in radish.

RESULTS

Three small RNA libraries from 'NAU-YH' taproot collected at pre-cortex splitting stage, cortex splitting stage and expanding stage were constructed. In all, 175 known and 107 potential novel miRNAs were discovered, from which 85 known and 13 novel miRNAs were found to be significantly differentially expressed during taproot thickening. Furthermore, totally 191 target genes were identified for the differentially expressed miRNAs. These target genes were annotated as transcription factors and other functional proteins, which were involved in various biological functions including plant growth and development, metabolism, cell organization and biogenesis, signal sensing and transduction, and plant defense response. RT-qPCR analysis validated miRNA expression patterns for five miRNAs and their corresponding target genes.

CONCLUSIONS

The small RNA populations of radish taproot at different thickening stages were firstly identified by Solexa sequencing. Totally 98 differentially expressed miRNAs identified from three taproot libraries might play important regulatory roles in taproot thickening. Their targets encoding transcription factors and other functional proteins including NF-YA2, ILR1, bHLH74, XTH16, CEL41 and EXPA9 were involved in radish taproot thickening. These results could provide new insights into the regulatory roles of miRNAs during the taproot thickening and facilitate genetic improvement of taproot in radish.

Network types and their application in natural variation studies in plants

We are in the age of data-driven biology. Not even a decade after the invention of high-throughput sequencing technologies, there are methods that accurately monitor DNA polymorphisms, transcription profiles, methylation states, transcription factor binding sites, chromatin compactness, nucleosome positions, dynamic histone marks, and so on. We are starting to generate comparable amounts of protein or metabolite data. A key issue is how are we going to make sense of all this information. Network analysis is the most promising method to integrate, query and display large amounts of data for human interpretation. This review shortly summarizes the basic types of networks, their properties and limitations. In addition, I introduce the application of networks to the study of the molecular mechanisms behind natural phenotypic variation.