Transcriptome/Degradome-wide discovery of microRNAs and transcript targets in two Paulownia australis genotypes

Integrated analysis of small RNAs, transcriptome and degradome sequencing reveal the drought stress network in Agropyron mongolicum Keng

Agropyron mongolicum (A. mongolicum) is an excellent gramineous forage with extreme drought tolerance, which lives in arid and semiarid desert areas. However, the mechanism that underlies the response of microRNAs (miRNAs) and their targets in A. mongolicum to drought stress is not well understood. In this study, we analyzed the transcriptome, small RNAome (specifically the miRNAome) and degradome to generate a comprehensive resource that focused on identifying key regulatory miRNA-target circuits under drought stress. The most extended transcript in each collection is known as the UniGene, and a total of 41,792 UniGenes and 1,104 miRNAs were identified, and 99 differentially expressed miRNAs negatively regulated 1,474 differentially expressed target genes. Among them, eight miRNAs were unique to A. mongolicum, and there were 36 target genes. A weighted gene co-expression network analysis identified five hub genes. The miRNAs of five hub genes were screened with an integration analysis of the degradome and sRNAs, such as osa-miR444a-3p.2-MADS47, bdi-miR408-5p_1ss19TA-CCX1, tae-miR9774_L-2R-1_1ss11GT-carC, ata-miR169a-3p-PAO2, and bdi-miR528-p3_2ss15TG20CA-HOX24. The functional annotations revealed that they were involved in mediating the brassinosteroid signal pathway, transporting and exchanging sodium and potassium ions and regulating the oxidation-reduction process, hydrolase activity, plant response to water deprivation, abscisic acid (ABA) and the ABA-activated signaling pathway to regulate drought stress. Five hub genes were discovered, which could play central roles in the regulation of drought-responsive genes. These results show that the combined analysis of miRNA, the transcriptome and degradation group provides a useful platform to investigate the molecular mechanism of drought resistance in A. mongolicum and could provide new insights into the genetic engineering of Poaceae crops in the future.

Analysis of small RNAs revealed differential expressions during pollen and embryo sac development in autotetraploid rice

Background

Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart.

Results

A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during pollen and embryo sac development in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of DEM associated with embryo sac and pollen development revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development, while only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development in autotetraploid rice. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression patterns during pollen and embryo sac development. Of these DEM, osa-miR1436_L + 3_1ss5CT and osa-miR167h-3p were associated with the female meiosis, while osa-miR159a.1 and osa-MIR159a-p5 were related with the male meiosis. 21 nt-phasiRNAs were detected during both pollen and embryo sac development, while 24 nt-phasiRNAs were found only in pollen development, which displayed down-regulation in autotetraploid compared to diploid rice and their spatial-temporal expression patterns were similar to osa-miR2275d. 24 nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development.

Conclusion

The above results not only provide the small RNAs changes during four landmark stages of pollen and embryo sac development in autotetraploid rice but also have identified specifically expressed miRNAs, especially meiosis-related miRNAs, pollen-specific-24 nt-phasiRNAs and TEs-siRNAs in autotetraploid rice. Together, these findings provide a foundation for understanding the effect of polyploidy on small RNAs expression patterns during pollen and embryo sac development that may lead to different abnormalities in autotetraploid rice.

Electronic supplementary material

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

Identification of novel miRNAs from drought tolerant rice variety Nagina 22

MicroRNAs regulate a spectrum of developmental and biochemical processes in plants and animals. Thus, knowledge of the entire miRNome is essential to understand the complete regulatory schema of any organism. The current study attempts to unravel yet undiscovered miRNA genes in rice. Analysis of small RNA libraries from various tissues of drought-tolerant ‘aus’ rice variety Nagina 22 (N22) identified 71 novel miRNAs. These were validated based on precursor hairpin structure, small RNA mapping pattern, ‘star’ sequence, conservation and identification of targets based on degradome data. While some novel miRNAs were conserved in other monocots and dicots, most appear to be lineage-specific. They were segregated into two different classes based on the closeness to the classical miRNA definition. Interestingly, evidence of a miRNA-like cleavage was found even for miRNAs that lie beyond the classical definition. Several novel miRNAs displayed tissue-enriched and/or drought responsive expression. Generation and analysis of the degradome data from N22 along with publicly available degradome identified several high confidence targets implicated in regulation of fundamental processes such as flowering and stress response. Thus, discovery of these novel miRNAs considerably expands the dimension of the miRNA-mediated regulation in rice.

Transcriptome and Degradome of microRNAs and Their Targets in Response to Drought Stress in the Plants of a Diploid and Its Autotetraploid Paulownia australis

MicroRNAs (miRNAs) are small, non-coding RNAs that play vital roles in plant growth, development, and stress response. Increasing numbers of studies aimed at discovering miRNAs and analyzing their functions in plants are being reported. In this study, we investigated the effect of drought stress on the expression of miRNAs and their targets in plants of a diploid and derived autotetraploid Paulownia australis. Four small RNA (sRNA) libraries and four degradome libraries were constructed from diploid and autotetraploid P. australis plants treated with either 75% or 25% relative soil water content. A total of 33 conserved and 104 novel miRNAs (processing precision value > 0.1) were identified, and 125 target genes were identified for 36 of the miRNAs by using the degradome sequencing. Among the identified miRNAs, 54 and 68 were differentially expressed in diploid and autotetraploid plants under drought stress (25% relative soil water content), respectively. The expressions of miRNAs and target genes were also validated by quantitative real-time PCR. The results showed that the relative expression trends of the randomly selected miRNAs were similar to the trends predicted by Illumina sequencing. And the correlations between miRNAs and their target genes were also analyzed. Furthermore, the functional analysis showed that most of these miRNAs and target genes were associated with plant development and environmental stress response. This study provided molecular evidence for the possible involvement of certain miRNAs in the drought response and/or tolerance in P. australis, and certain level of differential expression between diploid and autotetraploid plants.

Identification of microRNAs and their targets in Paulownia fortunei plants free from phytoplasma pathogen after methyl methane sulfonate treatment

MicroRNAs (miRNAs) play major roles in plant responses to various biotic and abiotic stresses by regulating gene expression at the transcriptional and post-transcriptional levels. Paulownia witches' broom (PaWB) disease caused by phytoplasmas reduces Paulownia production worldwide. In this study, we investigated the miRNA-mediated plant response to PaWB phytoplasma by Illumina sequencing and degradome analysis of Paulownia fortunei small RNAs (sRNAs). The sRNA and degradome libraries were constructed from healthy and diseased P. fortunei plants and the plants free from phytoplasma pathogen after 60 mg L(-1) methyl methane sulfonate treatment. A total of 96 P. fortunei-conserved miRNAs and 83 putative novel miRNAs were identified. Among them, 37 miRNAs (17 conserved, 20 novel) were found to be differentially expressed in response to PaWB phytoplasma infection. In addition, 114 target genes for 18 of the conserved miRNA families and 33 target genes for 15 of the novel miRNAs in P. fortunei were detected. The expression patterns of 14 of the PaWB phytoplasma-responsive miRNAs and 12 target genes were determined by quantitative real-time polymerase chain reaction (qPCR) experiments. A functional analysis of the miRNA targets indicated that these targeted genes may regulate transcription, stress response, nitrogen metabolism, and various other activities. Our results will help identify the potential roles of miRNAs involved in protecting P. fortunei from diseases.

Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulator Sedum alfredii uncovers a complex regulatory network and provides insights into cadmium phytoremediation

The hyperaccumulating ecotype of Sedum alfredii Hance is a cadmium (Cd)/zinc/lead co-hyperaccumulating species of Crassulaceae. It is a promising phytoremediation candidate accumulating substantial heavy metal ions without obvious signs of poisoning. However, few studies have focused on the regulatory roles of miRNAs and their targets in the hyperaccumulating ecotype of S. alfredii. Here, we combined analyses of the transcriptomics, sRNAs and the degradome to generate a comprehensive resource focused on identifying key regulatory miRNA-target circuits under Cd stress. A total of 87 721 unigenes and 356 miRNAs were identified by deep sequencing, and 79 miRNAs were differentially expressed under Cd stress. Furthermore, 754 target genes of 194 miRNAs were validated by degradome sequencing. A gene ontology (GO) enrichment analysis of differential miRNA targets revealed that auxin, redox-related secondary metabolism and metal transport pathways responded to Cd stress. An integrated analysis uncovered 39 pairs of miRNA targets that displayed negatively correlated expression profiles. Ten miRNA-target pairs also exhibited negative correlations according to a real-time quantitative PCR analysis. Moreover, a coexpression regulatory network was constructed based on profiles of differentially expressed genes. Two hub genes, ARF4 (auxin response factor 4) and AAP3 (amino acid permease 3), which might play central roles in the regulation of Cd-responsive genes, were uncovered. These results suggest that comprehensive analyses of the transcriptomics, sRNAs and the degradome provided a useful platform for investigating Cd hyperaccumulation in S. alfredii, and may provide new insights into the genetic engineering of phytoremediation.

Comparative Analysis and Identification of miRNAs and Their Target Genes Responsive to Salt Stress in Diploid and Tetraploid Paulownia fortunei Seedlings

Salt stress is a global environmental problem that affects plant growth and development. Paulownia fortunei is an adaptable and fast-growing deciduous tree native to China that is environmentally and economically important. MicroRNAs (miRNAs) play important regulatory roles in growth, development, and stress responses in plants. MiRNAs that respond to biotic stresses have been identified; however, how miRNAs in P. fortunei respond to salt stress has not yet been reported. To identify salt-stress-responsive miRNAs and predict their target genes, four small RNA and four degradome libraries were constructed from NaCl-treated and NaCl-free leaves of P. fortunei seedlings. The results indicated that salt stress had different physiological effects on diploid and tetraploid P. fortunei. We detected 53 conserved miRNAs belonging to 17 miRNA families and 134 novel miRNAs in P. fortunei. Comparing their expression levels in diploid and tetraploid P. fortunei, we found 10 conserved and 10 novel miRNAs that were significantly differentially expressed under salt treatment, among them eight were identified as miRNAs probably associated with higher salt tolerance in tetraploid P. fortunei than in diploid P. fortunei. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to predict the functions of the target genes of the conserved and novel miRNAs. The expressions of 10 differentially expressed miRNAs were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on P. fortunei miRNAs and their target genes under salt stress. The results provided information at the physiological and molecular levels for further research into the response mechanisms of P. fortunei to salt stress.

Transcriptome, microRNA, and degradome analyses of the gene expression of Paulownia with phytoplamsa

BACKGROUND

Paulownia witches' broom (PaWB) is a fatal disease of Paulownia caused by a phytoplasma. In previous studies, we found that plants with PaWB symptoms would revert to a healthy morphology after methyl methane sulfonate (MMS) treatment. To completely understand the gene expression profiles of the Paulownia-phytoplasma interaction, three high-throughput sequencing technologies were used to investigate changes of gene expression and microRNAs (miRNAs) in healthy Paulownia tomentosa plantlets, PaWB-infected plantlets, and PaWB-infected plantlets treated with 60 mg · L(-1) MMS.

METHODS

Transcriptome, miRNAs and degradome sequencing were performed to explore the global gene expression profiles in the process of Paulownia tomentosa with phytoplasma infection.

RESULTS

A total of 98,714 all-unigenes, 62 conserved miRNAs, and 35 novel miRNAs were obtained, among which 902 differentially expressed genes (DEGs) and 24 miRNAs were found to be associated with PaWB disease. Subsequently, the target genes of these miRNAs were predicted by degradome sequencing. Interestingly, we found that 19 target genes of these differentially expressed miRNAs were among the 902 DEGs. The targets of pau-miR156g, pau-miR403, and pau-miR166c were significantly up-regulated in the P. tomentosa plantlets infected with phytoplasma. Interaction of miRNA -target genes mediated gene expression related to PaWB were identified.

CONCLUSIONS

The results elucidated the possible roles of the regulation of genes and miRNAs in the Paulownia-phytoplasma interaction, which will enrich our understanding of the mechanisms of PaWB disease in this plant.

Plant-Pathogen Interaction-Related MicroRNAs and Their Targets Provide Indicators of Phytoplasma Infection in Paulownia tomentosa × Paulownia fortunei

Paulownia witches’ broom (PaWB) caused by a phytoplasma, has caused extensive losses in the yields of paulownia timber and resulted in significant economic losses. However, the molecular mechanisms in Paulownia that underlie the phytoplasma stress are poorly characterized. In this study, we use an Illumina platform to sequence four small RNA libraries and four degradome sequencing libraries derived from healthy, PaWB-infected, and PaWB-infected 15 mg·L⁻¹ and 30 mg·L⁻¹ methyl methane sulfonate (MMS)-treated plants. In total, 125 conserved and 118 novel microRNAs (miRNAs) were identified and 33 miRNAs responsive to PaWB disease were discovered. Furthermore, 166 target genes for 18 PaWB disease-related miRNAs were obtained, and found to be involved in plant-pathogen interaction and plant hormone signal transduction metabolic pathways. Eleven miRNAs and target genes responsive to PaWB disease were examined by a quantitative real-time PCR approach. Our findings will contribute to studies on miRNAs and their targets in Paulownia, and provide new insights to further understand plant-phytoplasma interactions.

Discovery of microRNAs and transcript targets related to witches' broom disease in Paulownia fortunei by high-throughput sequencing and degradome approach

Paulownia witches' broom (PaWB) caused by the phytoplasma is a devastating disease of Paulownia trees. It has caused heavy yield losses to Paulownia production worldwide. However, knowledge of the transcriptional and post-transcriptional regulation of gene expression by microRNAs (miRNAs), especially miRNAs responsive to PaWB disease stress, is still rudimentary. In this study, to identify miRNAs and their transcript targets that are responsive to PaWB disease stress, six sequencing libraries were constructed from healthy (PF), PaWB-infected (PFI), and PaWB-infected, 20 mg L(-1) methyl methane sulfonate-treated (PFI20) P. fortunei seedlings. As a result, 95 conserved miRNAs belonging to 18 miRNA families, as well as 122 potential novel miRNAs, were identified. Most of them were found to be a response to PaWB disease-induced stress, and the expression levels of these miRNAs were validated by quantitative real-time PCR analysis. The study simultaneously identified 109 target genes from the P. fortunei for 14 conserved miRNA families and 24 novel miRNAs by degradome sequencing. Furthermore, the functions of the miRNA targets were annotated based on Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis. The results presented here provide the groundwork for further analysis of miRNAs and target genes responsive to the PaWB disease stress, and could be also useful for addressing new questions to better understand the mechanisms of plant infection by phytoplasma in the future.