Characterization of the small RNA component of leaves and fruits from four different cucurbit species

Schwann cells-derived exosomal miR-21 participates in high glucose regulation of neurite outgrowth

As a common complication of diabetes, the pathogenesis of diabetic peripheral neuropathy (DPN) is closely related to high glucose but has not been clarified. Exosomes can mediate crosstalk between Schwann cells (SC) and neurons in the peripheral nerve. Herein, we found that miR-21 in serum exosomes from DPN rats was decreased. SC proliferation was inhibited, cell apoptosis was increased, and the expression of miR-21 in cells and exosomes was downregulated when cultured in high glucose. Increasing miR-21 expression reversed these changes, while knockdown of miR-21 led to the opposite results. When co-cultured with exosomes derived from SC exposed to high glucose, neurite outgrowth was inhibited. On the contrary, neurite outgrowth was accelerated when incubated with exosomes rich in miR-21. We further demonstrated that the SC-derived exosomal miR-21 participates in neurite outgrowth probably through the AKT signaling pathway. Thus, SC-derived exosomal miR-21 contributes to high glucose regulation of neurite outgrowth.

•

The miR-21 was decreased in serum exosomes and sciatic nerve of DPN rats

•

High glucose inhibited SC viability and downregulated the expression of miR-21

•

Exosomes derived from SC cultured in high glucose inhibited the neurite outgrowth

•

SC-derived exosomes rich in miR-21 accelerated the neurite outgrowth of neuron

Combined Stress Conditions in Melon Induce Non-additive Effects in the Core miRNA Regulatory Network

Climate change has been associated with a higher incidence of combined adverse environmental conditions that can promote a significant decrease in crop productivity. However, knowledge on how a combination of stresses might affect plant development is still scarce. MicroRNAs (miRNAs) have been proposed as potential targets for improving crop productivity. Here, we have combined deep-sequencing, computational characterization of responsive miRNAs and validation of their regulatory role in a comprehensive analysis of response of melon to several combinations of four stresses (cold, salinity, short day, and infection with a fungus). Twenty-two miRNA families responding to double and/or triple stresses were identified. The regulatory role of the differentially expressed miRNAs was validated by quantitative measurements of the expression of the corresponding target genes. A high proportion (ca. 60%) of these families (mainly highly conserved miRNAs targeting transcription factors) showed a non-additive response to multiple stresses in comparison with that observed under each one of the stresses individually. Among those miRNAs showing non-additive response to stress combinations, most interactions were negative, suggesting the existence of functional convergence in the miRNA-mediated response to combined stresses. Taken together, our results provide compelling pieces of evidence that the response to combined stresses cannot be easily predicted from the study individual stresses.

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

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

BACKGROUND

MicroRNA-mediated gene regulatory networks play a significant role in plant growth and development and environmental stress responses.

RESULTS

We identified 79 microRNAs (miRNAs) and multiple miRNA variants (isomiRs) belonging to 26 miRNA families in the primary root growth zone of maize seedlings grown at one of three water potentials: well-watered (- 0.02 MPa), mild water deficit stress (- 0.3 MPa), and severe water deficit stress (- 1.6 MPa). The abundances of 3 miRNAs (mild stress) and 34 miRNAs representing 17 families (severe stress) were significantly different in water-deficit stressed relative to well-watered controls (FDR < 0.05 and validated by stem loop RT-qPCR). Degradome sequencing revealed 213 miRNA-regulated transcripts and trancriptome profiling revealed that the abundance of 77 (miRNA-regulated) were regulated by water-defecit stress. miR399e,i,j-3p was strongly regulated by water-defcit stress implicating the possibility of nutrient deficiency during stress.

CONCLUSIONS

We have identified a number of maize miRNAs that respond to specific water deficits applied to the primary root growth zone. We have also identified transcripts that are targets for miRNA regulation in the root growth zone under water-deficit stress. The miR399e,i,j-3p that is known to regulate phosphate uptake in response to nutrient deficiencies responds to water-deficit stress, however, at the seedling stage the seed provides adequate nutrients for root growth thus miR399e,i,j-3p may play a separate role in water-deficit responses. A water-deficit regulated maize transcript, similar to known miR399 target mimics, was identified and we hypothesized that it is another regulatory player, moderating the role of miR399e,i,j-3p, in primary root growth zone water deficit responses.

Identification and Temporal Expression Analysis of Conserved and Novel MicroRNAs in the Leaves of Winter Wheat Grown in the Field

Cold acclimation and vegetative/reproductive transition are two important evolutionary adaptive mechanisms for winter wheat surviving the freezing temperature in winter and successful seeds setting in the next year. MicroRNA (miRNA) is a class of regulatory small RNAs (sRNAs), which plays critical roles in the growth and development of plants. However, the regulation mechanism of miRNAs during cold acclimation and vegetative/reproductive transition of winter wheat is not much understood. In this study, four sRNA libraries from leaves of winter wheat grown in the field at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage were analyzed to identify known and novel miRNAs and to understand their potential roles in the growth and development of winter wheat. We examined miRNA expression using a high-throughput sequencing technique. A total of 373 known, 55 novel, and 27 putative novel miRNAs were identified. Ninety-one miRNAs were found to be differentially expressed at the four stages. Among them, the expression of six known and eight novel miRNAs was significantly suppressed at the winter dormancy stage, whereas the expression levels of seven known and eight novel miRNAs were induced at this stage; three known miRNAs and three novel miRNAs were significantly induced at the spring green-up stage; six known miRNAs were induced at the spring green-up stage and reached the highest expression level at the jointing stage; and 20 known miRNAs and 10 novel miRNAs were significantly induced at the jointing stage. Expression of a number of representative differentially expressed miRNAs was verified using quantitative real-time polymerase chain reaction (qRT-PCR). Potential target genes for known and novel miRNAs were predicted. Moreover, six novel target genes for four Pooideae species-specific miRNAs and two novel miRNAs were verified using the RNA ligase-mediated 5'-rapid amplification of cDNA ends (RLM-5'RACE) technique. These results indicate that miRNAs are key non-coding regulatory factors modulating the growth and development of wheat. Our study provides valuable information for in-depth understanding of the regulatory mechanism of miRNAs in cold acclimation and vegetative/reproductive transition of winter wheat grown in the field.

Novel and Conserved miRNAs Among Brazilian Pine and Other Gymnosperms

The knowledge about plant miRNAs has increased exponentially, with thousands of miRNAs been reported in different plant taxa using high throughput sequencing technologies and bioinformatic tools. Nevertheless, several groups of plants remain unexplored, and the gap of knowledge about conifer miRNAs is considerable. There is no sequence or functional information available on miRNAs in Araucariaceae. This group is represented in Brazil by only one species, Araucaria angustifolia, an endangered species known as Brazilian pine. In the present study, Brazilian pine has its transcriptome explored with respect to small RNAs, representing the first description in a member of the Araucariaceae family. The screening for conserved miRNAs in Brazilian pine revealed 115 sequences of 30 miRNA families. A total of 106 precursors sequences were predicted. Forty one comprised conserved miRNAs from 16 families, whereas 65 were annotated as novel miRNAs. The comparison of Brazilian pine precursors with sRNA libraries of other five conifer species indicates that 9 out 65 novel miRNAs are conserved among gymnosperms, while 56 seems to be specific for Brazilian pine or restricted to Araucariaceae family. Analysis comparing novel Brazilian pine miRNAs precursors and Araucaria cunninghamii RNA-seq data identified seven orthologs between both species. Mature miRNA identified by bioinformatics predictions were validated using stem-loop RT-qPCR assays. The expression pattern of conserved and novel miRNAs was analyzed in five different tissues of 3-month-old Araucaria seedlings. The present study provides insights about the nature and composition of miRNAs in an Araucariaceae species, with valuable information on miRNAs diversity and conservation in this taxon.

Inferring the regulatory network of the miRNA-mediated response to biotic and abiotic stress in melon

BACKGROUND

MiRNAs have emerged as key regulators of stress response in plants, suggesting their potential as candidates for knock-in/out to improve stress tolerance in agricultural crops. Although diverse assays have been performed, systematic and detailed studies of miRNA expression and function during exposure to multiple environments in crops are limited.

RESULTS

Here, we present such pioneering analysis in melon plants in response to seven biotic and abiotic stress conditions. Deep-sequencing and computational approaches have identified twenty-four known miRNAs whose expression was significantly altered under at least one stress condition, observing that down-regulation was preponderant. Additionally, miRNA function was characterized by high scale degradome assays and quantitative RNA measurements over the intended target mRNAs, providing mechanistic insight. Clustering analysis provided evidence that eight miRNAs showed a broad response range under the stress conditions analyzed, whereas another eight miRNAs displayed a narrow response range. Transcription factors were predominantly targeted by stress-responsive miRNAs in melon. Furthermore, our results show that the miRNAs that are down-regulated upon stress predominantly have as targets genes that are known to participate in the stress response by the plant, whereas the miRNAs that are up-regulated control genes linked to development.

CONCLUSION

Altogether, this high-resolution analysis of miRNA-target interactions, combining experimental and computational work, Illustrates the close interplay between miRNAs and the response to diverse environmental conditions, in melon.

Genotype- and tissue-specific miRNA profiles and their targets in three alfalfa (Medicago sativa L) genotypes

BACKGROUND

Alfalfa (Medicago sativa L.) is a forage legume with significant agricultural value worldwide. MicroRNAs (miRNAs) are key components of post-transcriptional gene regulation and essentially regulate many aspects of plant growth and development. Although miRNAs were reported in alfalfa, their expression profiles in different tissues and the discovery of novel miRNAs as well as their targets have not been described in this plant species.

RESULTS

To identify tissue-specific miRNA profiles in whole plants, shoots and roots of three different alfalfa genotypes (Altet-4, NECS-141and NF08ALF06) were used. Small RNA libraries were generated and sequenced using a high-throughput sequencing platform. Analysis of these libraries enabled identification of100 miRNA families; 21 of them belong to the highly conserved families while the remaining 79 families are conserved at the minimum between M. sativa and the model legume and close relative, M. truncatula. The profiles of the six abundantly expressed miRNA families (miR156, miR159, miR166, miR319, miR396 and miR398) were relatively similar between the whole plants, roots and shoots of these three alfalfa genotypes. In contrast, robust differences between shoots and roots for miR160 and miR408 levels were evident, and their expression was more abundant in the shoots. Additionally, 17 novel miRNAs were identified and the relative abundance of some of these differed between tissue types. Further, the generation and analysis of degradome libraries from the three alfalfa genotypes enabled confirmation of 69 genes as targets for 31 miRNA families in alfalfa.

CONCLUSIONS

The miRNA profiles revealed both similarities and differences in the expression profiles between tissues within a genotype as well as between the genotypes. Among the highly conserved miRNA families, miR166 was the most abundantly expressed in almost all tissues from the three genotypes. The identification of conserved and novel miRNAs as well as their targets in different tissues of multiple genotypes increased our understanding of miRNA-mediated gene regulation in alfalfa and could provide valuable insights for practical research and plant improvement applications in alfalfa and related legume species.

In silico identification and computational characterization of endogenous small interfering RNAs from diverse grapevine tissues and stages

Small interfering RNAs (siRNAs) are effectors of regulatory pathways underlying plant development, metabolism, and stress- and nutrient-signaling regulatory networks. The endogenous siRNAs are generally not conserved between plants; consequently, it is necessary and important to identify and characterize siRNAs from various plants. To address the nature and functions of siRNAs, and understand the biological roles of the huge siRNA population in grapevine (Vitis vinifera L.). The high-throughput sequencing technology was used to identify a large set of putative endogenous siRNAs from six grapevine tissues/organs. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to classify the target genes of siRNA. In total, 520,519 candidate siRNAs were identified and their expression profiles exhibited typical temporal characters during grapevine development. In addition, we identified two grapevine trans-acting siRNA (TAS) gene homologs (VvTAS3 and VvTAS4) and the derived trans-acting siRNAs (tasiRNAs) that could target grapevine auxin response factor (ARF) and myeloblastosis (MYB) genes. Furthermore, the GO and KEGG analysis of target genes showed that most of them covered a broad range of functional categories, especially involving in disease-resistance process. The large-scale and completely genome-wide level identification and characterization of grapevine endogenous siRNAs from the diverse tissues by high throughput technology revealed the nature and functions of siRNAs in grapevine.

Genome-wide identification and comprehensive analysis of microRNAs and phased small interfering RNAs in watermelon

BACKGROUND

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs involved in the post-transcriptional gene regulation and play a critical role in plant growth, development and stress responses. Watermelon (Citrullus lanatus L.) is one of the important agricultural crops worldwide. However, the watermelon miRNAs and phasiRNAs and their functions are not well explored.

RESULTS

Here we carried out computational and experimental analysis of miRNAs and phased small interfering RNAs (phasiRNAs) in watermelon by analyzing 14 small RNA profiles from roots, leaves, androecium, petals, and fruits, and one published small RNA profile of mixed tissues. To identify the targets of miRNAs and phasiRNAs, we generated a degradome profile for watermelon leaf which is analyzed using the SeqTar algorithm. We identified 97 conserved pre-miRNAs, of which 58 have not been reported previously and 348 conserved mature miRNAs without precursors. We also found 9 novel pre-miRNAs encoding 18 mature miRNAs. One hundred and one 21 nucleotide (nt) PHAS loci, and two hundred and forty one 24 nt PHAS loci were also identified. We identified 127 conserved targets of the conserved miRNAs and TAS3-derived tasiRNAs by analyzing a degradome profile of watermelon leaf.

CONCLUSIONS

The presented results provide a comprehensive view of small regulatory RNAs and their targets in watermelon.

Comparative Analysis of Fruit Ripening-Related miRNAs and Their Targets in Blueberry Using Small RNA and Degradome Sequencing

MicroRNAs (miRNAs) play vital roles in the regulation of fruit development and ripening. Blueberry is an important small berry fruit crop with economical and nutritional value. However, nothing is known about the miRNAs and their targets involved in blueberry fruit ripening. In this study, using high-throughput sequencing of small RNAs, 84 known miRNAs belonging to 28 families and 16 novel miRNAs were identified in white fruit (WF) and blue fruit (BF) libraries, which represent fruit ripening onset and in progress, respectively. Among them, 41 miRNAs were shown to be differentially expressed during fruit maturation, and 16 miRNAs representing 16 families were further chosen to validate the sRNA sequencing data by stem-loop qRT-PCR. Meanwhile, 178 targets were identified for 41 known and 7 novel miRNAs in WF and BF libraries using degradome sequencing, and targets of miR160 were validated using RLM-RACE (RNA Ligase-Mediated (RLM)-Rapid Amplification of cDNA Ends) approach. Moreover, the expression patterns of 6 miRNAs and their targets were examined during fruit development and ripening. Finally, integrative analysis of miRNAs and their targets revealed a complex miRNA-mRNA regulatory network involving a wide variety of biological processes. The findings will facilitate future investigations of the miRNA-mediated mechanisms that regulate fruit development and ripening in blueberry.

Small RNA profiles from Panax notoginseng roots differing in sizes reveal correlation between miR156 abundances and root biomass levels

Plant genomes encode several classes of small regulatory RNAs (sRNAs) that play critical roles in both development and stress responses. Panax notoginseng (Burk.) F.H. Chen (P. notoginseng) is an important traditional Chinese herbal medicinal plant species for its haemostatic effects. Therefore, the root yield of P. notoginseng is a major economically important trait since the roots of P. notoginseng are the parts used to produce medicine. To identify sRNAs that are critical for the root biomass of P. notoginseng, we performed a comprehensive study of miRNA transcriptomes from P. notoginseng roots of different biomasses. We identified 675 conserved miRNAs, of which 180 pre-miRNAs are also identified, and three TAS3 loci in P. notoginseng. By using degradome sequencing, we identified 79 conserved miRNA:target or tasiRNA:target interactions, of which eight were further confirmed with the RLM 5'-RACE experiments. More importantly, our results revealed that a member of miR156 family and one of its SPL target genes have inverse expression levels, which is tightly correlated with greater root biomass contents. These results not only contributes to overall understanding of post-transcriptional gene regulation in roots of P. notoginseng but also could serve as markers for breeding P. notoginseng with greater root yield.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

Background

microRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.

Results

In this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.

Conclusions

Our results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Electronic supplementary material

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

MicroRNA dynamics in a wild and cultivated species of Convolvulaceae exposed to drought stress

Agricultural productivity is severely hampered by drought in many parts of the globe. It is well-known that wild plant species can tolerate drought better when compared with their closely related cultivated plant species. Better drought adaptation of wild species over cultivated ones is accounted for their ability to differentially regulate gene expression. miRNAs, known to regulate gene expression at the post-transcriptional level, are admitted to play an important role in plant adaptation to stresses. This study aims at evaluating miRNA dynamics in a drought-tolerant wild Ipomoea campanulata L. and drought-sensitive cultivated Jacquemontia pentantha (Jacq.) of the family Convolvulaceae under ex situ drought. Sequencing profiles revealed that 34 conserved miRNA families were analogous between the two species. Drought altered expression levels of several of these miRNAs in both the species. Drought-tolerant I. campanulata showed upregulation of miR398, miR168, miR858, miR162 and miR408, while miR394 and miR171 were downregulated. Drought-sensitive J. pentantha showed upregulation of miR394, miR156, miR160, miR164, miR167, miR172, miR319, miR395, miR396, miR403 and downregulation of miR157. Basal miRNA levels and their drought mediated regulation were very different between the two species. Differential drought sensitivities of these two plant species can be attributed to these innate variations in miRNA levels and their expression.

Intron Lariat RNA Inhibits MicroRNA Biogenesis by Sequestering the Dicing Complex in Arabidopsis

Lariat RNAs formed as by-products of splicing are quickly degraded by the RNA debranching enzyme 1 (DBR1), leading to their turnover. Null dbr1 mutants in both animals and plants are embryo lethal, but the mechanism underlying the lethality remains unclear. Here we characterized a weak mutant allele of DBR1 in Arabidopsis, dbr1-2, and showed that a global increase in lariat RNAs was unexpectedly accompanied by a genome-wide reduction in miRNA accumulation. The dbr1-2 mutation had no effects on expression of miRNA biogenesis genes or primary miRNAs (pri-miRNAs), but the association of pri-miRNAs with the DCL1/HYL1 dicing complex was impaired. Lariat RNAs were associated with the DCL1/HYL1 dicing complex in vivo and competitively inhibited the binding of HYL1 with pri-miRNA. Consistent with the impacts of lariat RNAs on miRNA biogenesis, over-expression of lariat RNAs reduced miRNA accumulation. Lariat RNAs localized in nuclear bodies, and partially co-localize with HYL1, and both DCL1 and HYL1 were mis-localized in dbr1-2. Together with our findings that nearly four hundred lariat RNAs exist in wild type plants and that these lariat RNAs also associate with the DCL1/HYL1 dicing complex in vivo, we thus propose that lariat RNAs, as decoys, inhibit miRNA processing, suggesting a hitherto unknown layer of regulation in miRNA biogenesis.

It is known that lariat RNAs formed during pre-mRNA splicing are debranched by DBR1 (RNA debranching enzyme 1). Loss of function of DBR1 causes embryo lethality in both animals and plants. In animals, some debranched lariat RNAs could be further processed into mirtron miRNAs, a class of nonconventional miRNAs that bypass the microprocessor for their biogenesis. However, no mirtron has been functionally validated in plants, and how the accumulation of lariat RNA in dbr1 results in embryo lethality remains unclear. Here, we show that DBR1 is necessary for the regulation of genome-wide miRNA biogenesis in plants. By investigating the correlation between lariat RNA accumulation and miRNA processing, we showed that the DBR1-mediated lariat RNA debranching process provides a safeguard role for the binding of the dicing complex with miRNA precursors. As both the DBR1-mediated lariat RNA debranching process and miRNA biogenesis are common features in higher eukaryotes, the finding that lariat RNAs sequester the dicing complex in plants may have a broad implications for the non-coding RNA field.

Integrated mRNA and miRNA transcriptome reveal a cross-talk between developing response and hormone signaling for the seed kernels of Siberian apricot

Recently, our transcriptomic analysis has identified some functional genes responsible for oil biosynthesis in developing SASK, yet miRNA-mediated regulation for SASK development and oil accumulation is poorly understood. Here, 3 representative periods of 10, 30 and 60 DAF were selected for sRNA sequencing based on the dynamic patterns of growth tendency and oil content of developing SASK. By miRNA transcriptomic analysis, we characterized 296 known and 44 novel miRNAs in developing SASK, among which 36 known and 6 novel miRNAs respond specifically to developing SASK. Importantly, we performed an integrated analysis of mRNA and miRNA transcriptome as well as qRT-PCR detection to identify some key miRNAs and their targets (miR156-SPL, miR160-ARF18, miR164-NAC1, miR171h-SCL6, miR172-AP2, miR395-AUX22B, miR530-P2C37, miR393h-TIR1/AFB2 and psi-miRn5-SnRK2A) potentially involved in developing response and hormone signaling of SASK. Our results provide new insights into the important regulatory function of cross-talk between development response and hormone signaling for SASK oil accumulation.

Genome-wide identification of microRNAs in pomegranate (Punica granatum L.) by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs), a class of small non-coding endogenous RNAs that regulate gene expression post-transcriptionally, play multiple key roles in plant growth and development and in biotic and abiotic stress response. Knowledge and roles of miRNAs in pomegranate fruit development have not been explored.

RESULTS

Pomegranate, which accumulates a large amount of anthocyanins in skin and arils, is valuable to human health, mainly because of its antioxidant properties. In this study, we developed a small RNA library from pooled RNA samples from young seedlings to mature fruits and identified both conserved and pomegranate-specific miRNA from 29,948,480 high-quality reads. For the pool of 15- to 30-nt small RNAs, ~50 % were 24 nt. The miR157 family was the most abundant, followed by miR156, miR166, and miR168, with variants within each family. The base bias at the first position from the 5' end had a strong preference for U for most 18- to 26-nt sRNAs but a preference for A for 18-nt sRNAs. In addition, for all 24-nt sRNAs, the nucleotide U was preferred (97 %) in the first position. Stem-loop RT-qPCR was used to validate the expression of the predominant miRNAs and novel miRNAs in leaves, male and female flowers, and multiple fruit developmental stages; miR156, miR156a, miR159a, miR159b, and miR319b were upregulated during the later stages of fruit development. Higher expression of miR156 in later fruit developmental may positively regulate anthocyanin biosynthesis by reducing SPL transcription factor. Novel miRNAs showed variation in expression among different tissues. These novel miRNAs targeted different transcription factors and hormone related regulators. Gene ontology and KEGG pathway analyses revealed predominant metabolic processes and catalytic activities, important for fruit development. In addition, KEGG pathway analyses revealed the involvement of miRNAs in ascorbate and linolenic acid, starch and sucrose metabolism; RNA transport; plant hormone signaling pathways; and circadian clock.

CONCLUSION

Our first and preliminary report of miRNAs will provide information on the synthesis of biochemical compounds of pomegranate for future research. The functions of the targets of the novel miRNAs need further investigation.

Transfer RNA Derived Small RNAs Targeting Defense Responsive Genes Are Induced during Phytophthora capsici Infection in Black Pepper (Piper nigrum L.)

Small RNAs derived from transfer RNAs were recently assigned as potential gene regulatory candidates for various stress responses in eukaryotes. In this study, we report on the cloning and identification of tRNA derived small RNAs from black pepper plants in response to the infection of the quick wilt pathogen, Phytophthora capsici. 5'tRFs cloned from black pepper were validated as highly expressed during P. capsici infection. A high-throughput systematic analysis of the small RNAome (sRNAome) revealed the predominance of 5'tRFs in the infected leaf and root. The abundance of 5'tRFs in the sRNAome and the defense responsive genes as their potential targets indicated their regulatory role during stress response in black pepper. The 5'Ala(CGC) tRF mediated cleavage was experimentally mapped at the tRF binding sites on the mRNA targets of Non-expresser of pathogenesis related protein (NPR1), which was down-regulated during pathogen infection. Comparative sRNAome further demonstrated sequence conservation of 5'Ala tRFs across the angiosperm plant groups, and many important genes in the defense response were identified in silico as their potential targets. Our findings uncovered the diversity, differential expression and stress responsive functional role of tRNA-derived small RNAs during Phytophthora infection in black pepper.

Characterization of drought- and heat-responsive microRNAs in switchgrass

Recent investigations revealed that microRNAs (miRNAs) play crucial roles in plant acclimation to stress conditions. Switchgrass, one of the important biofuel crop species can withstand hot and dry climates but the molecular basis of stress tolerance is relatively unknown. To identify miRNAs that are important for tolerating drought or heat, small RNAs were profiled in leaves of adult plants exposed to drought or heat. Sequence analysis enabled the identification of 29 conserved and 62 novel miRNA families. Notably, the abundances of several conserved and novel miRNAs were dramatically altered following drought or heat. Using at least one fold (log2) change as cut off, we observed that 13 conserved miRNA families were differentially regulated by both stresses, and, five and four families were specifically regulated by drought and heat, respectively. Similarly, using a more stringent cut off of two fold (log2) regulation, we found 5 and 16 novel miRNA families were upregulated but 6 and 7 families were downregulated under drought and heat, respectively. The stress-altered expression of a subset of miRNAs and their targets was confirmed using quantitative PCR. Overall, the switchgrass plants exposed to drought or heat revealed similarities as well as differences with respect to miRNA regulation, which could be important for enduring different stress conditions.

Transcriptome profiling of differentially expressed genes in floral buds and flowers of male sterile and fertile lines in watermelon

Background

Male sterility is an important mechanism for the production of hybrid seeds in watermelon. Although fruit development has been studied extensively in watermelon, there are no reports on gene expression in floral organs. In this study, RNA-sequencing (RNA-seq) was performed in two near-isogenic watermelon lines (genic male sterile [GMS] line, DAH3615-MS and male fertile line, DAH3615) to identify the differentially expressed genes (DEGs) related to male sterility.

Results

DEG analysis showed that 1259 genes were significantly associated with male sterility at a FDR P-value of < 0.01. Most of these genes were only expressed in the male fertile line. In addition, 11 functional clusters were identified using DAVID functional classification analysis. Of detected genes in RNA-seq analysis, 19 were successfully validated by qRT-PCR.

Conclusions

In this study, we carried out a comprehensive floral transcriptome sequence comparison of a male fertile line and its near-isogenic male sterile line in watermelon. This analysis revealed essential genes responsible for stamen development, including pollen development and pollen tube elongation, and allowed their functional classification. These results provided new information on global mechanisms related to male sterility in watermelon.

Electronic supplementary material

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

Identification of miRNAs with potential roles in regulation of anther development and male-sterility in 7B-1 male-sterile tomato mutant

BACKGROUND

The 7B-1 tomato line (Solanum lycopersicum cv. Rutgers) is a photoperiod-sensitive male-sterile mutant, with potential application in hybrid seed production. Small RNAs (sRNAs) in tomato have been mainly characterized in fruit development and ripening, but none have been studied with respect to flower development and regulation of male-sterility. Using sRNA sequencing, we identified miRNAs that are potentially involved in anther development and regulation of male-sterility in 7B-1 mutant.

RESULTS

Two sRNA libraries from 7B-1 and wild type (WT) anthers were sequenced and thirty two families of known miRNAs and 23 new miRNAs were identified in both libraries. MiR390, miR166, miR159 were up-regulated and miR530, miR167, miR164, miR396, miR168, miR393, miR8006 and two new miRNAs, miR#W and miR#M were down-regulated in 7B-1 anthers. Ta-siRNAs were not differentially expressed and likely not associated with 7B-1 male-sterility. miRNA targets with potential roles in anther development were validated using 5'-RACE. QPCR analysis showed differential expression of miRNA/target pairs of interest in anthers and stem of 7B-1, suggesting that they may regulate different biological processes in these tissues. Expression level of most miRNA/target pairs showed negative correlation, except for few. In situ hybridization showed predominant expression of miR159, GAMYBL1, PMEI and cystatin in tapetum, tetrads and microspores.

CONCLUSION

Overall, we identified miRNAs with potential roles in anther development and regulation of male-sterility in 7B-1. A number of new miRNAs were also identified from tomato for the first time. Our data could be used as a benchmark for future studies of the molecular mechanisms of male-sterility in other crops.

High Throughput Sequencing of Small RNAs in the Two Cucurbita Germplasm with Different Sodium Accumulation Patterns Identifies Novel MicroRNAs Involved in Salt Stress Response

MicroRNAs (miRNAs), a class of small non-coding RNAs, recognize their mRNA targets based on perfect sequence complementarity. MiRNAs lead to broader changes in gene expression after plants are exposed to stress. High-throughput sequencing is an effective method to identify and profile small RNA populations in non-model plants under salt stresses, significantly improving our knowledge regarding miRNA functions in salt tolerance. Cucurbits are sensitive to soil salinity, and the Cucurbita genus is used as the rootstock of other cucurbits to enhance salt tolerance. Several cucurbit crops have been used for miRNA sequencing but salt stress-related miRNAs in cucurbit species have not been reported. In this study, we subjected two Cucurbita germplasm, namely, N12 (Cucurbita. maxima Duch.) and N15 (Cucurbita. moschata Duch.), with different sodium accumulation patterns, to Illumina sequencing to determine small RNA populations in root tissues after 4 h of salt treatment and control. A total of 21,548,326 and 19,394,108 reads were generated from the control and salt-treated N12 root tissues, respectively. By contrast, 19,108,240 and 20,546,052 reads were obtained from the control and salt-treated N15 root tissues, respectively. Fifty-eight conserved miRNA families and 33 novel miRNAs were identified in the two Cucurbita germplasm. Seven miRNAs (six conserved miRNAs and one novel miRNAs) were up-regulated in salt-treated N12 and N15 samples. Most target genes of differentially expressed novel miRNAs were transcription factors and salt stress-responsive proteins, including dehydration-induced protein, cation/H⁺ antiporter 18, and CBL-interacting serine/threonine-protein kinase. The differential expression of miRNAs between the two Cucurbita germplasm under salt stress conditions and their target genes demonstrated that novel miRNAs play an important role in the response of the two Cucurbita germplasm to salt stress. The present study initially explored small RNAs in the response of pumpkin to salt stress, and provided valuable information on novel miRNAs and their target genes in Cucurbita.

Grafting-responsive miRNAs in cucumber and pumpkin seedlings identified by high-throughput sequencing at whole genome level

Grafting is an important agricultural technique widely used for improving growth, yields and tolerance of crops to abiotic and biotic stresses. As one type of endogenous, non-coding small RNAs, microRNAs (miRNAs) regulate development and responsiveness to biotic and abiotic stresses by negatively mediating expression of target genes at the post-transcriptional level. However, there have been few detailed studies to evaluate the role of miRNAs in mediation of grafting-induced physiological processes in plants. Cucumis sativus and Cucurbita moschata are important vegetables worldwide. We constructed eight small RNA libraries from leaves and roots of seedlings that were grafted in the following four ways: (1) hetero-grafting, using cucumber as scion and pumpkin as rootstock; (2) hetero-grafting, with pumpkin as scion and cucumber as rootstock; (3) auto-grafting of cucumbers and (4) auto-grafting of pumpkins. High-throughput sequencing was employed, and more than 120 million raw reads were obtained. We annotated 112 known miRNAs belonging to 40 miRNA families and identified 48 new miRNAs in the eight libraries, and the targets of these known and novel miRNAs were predicted by bioinformatics. Grafting led to changes in expression of most miRNAs and their predicted target genes, suggesting that miRNAs may play significant roles in mediating physiological processes of grafted seedlings by regulating the expression of target genes. The potential role of the grafting-responsive miRNAs in seedling growth and long-distance transport of miRNA was discussed. These results are useful for functional characterization of miRNAs in mediation of grafting-dependent physiological processes.

MiR-182 is up-regulated and targeting Cebpa in hepatocellular carcinoma

MicroRNAs (miRNAs) are endogenous small non-coding RNAs that repress their targets at post transcriptional level. Existing studies have shown that miRNAs are important regulatory genes in hepatocellular carcinoma (HCC), as either tumor suppressors or oncogenes. MiR-122 is normally downregulated in HCC and regarded as a tumor suppressor. Recently miR-122 has been reported to be regulated by CEBPA, which is then involved in a novel pathway to influence proliferation of tumor cells. However it is unknown whether CEBPA is regulated by miRNAs in HCC. In this study, we find that miR-182 is upregulated in HCC model rat, and represses CEBPA in both rat and human. This further improves the current CEBPA/miR-122 pathway that controls the proliferation of tumor cells. These results suggest that miR-182 is a potential oncogene in HCC and could be used as a diagnostic marker and drug target of HCC.

Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-20) contains supplementary material, which is available to authorized users.

Peculiar evolutionary history of miR390-guided TAS3-like genes in land plants

PCR-based approach was used as a phylogenetic profiling tool to probe genomic DNA samples from representatives of evolutionary distant moss taxa, namely, classes Bryopsida, Tetraphidopsida, Polytrichopsida, Andreaeopsida, and Sphagnopsida. We found relatives of all Physcomitrella patens miR390 and TAS3-like loci in these plant taxa excluding Sphagnopsida. Importantly, cloning and sequencing of Marchantia polymorpha genomic DNA showed miR390 and TAS3-like sequences which were also found among genomic reads of M. polymorpha at NCBI database. Our data suggest that the ancient plant miR390-dependent TAS molecular machinery firstly evolved to target AP2-like mRNAs in Marchantiophyta and only then both ARF- and AP2-specific mRNAs in mosses. The presented analysis shows that moss TAS3 families may undergone losses of tasiAP2 sites during evolution toward ferns and seed plants. These data confirm that miR390-guided genes coding for ARF- and AP2-specific ta-siRNAs have been gradually changed during land plant evolution.

The hot pepper (Capsicum annuum) microRNA transcriptome reveals novel and conserved targets: a foundation for understanding MicroRNA functional roles in hot pepper

MicroRNAs (miRNAs) are a class of non-coding RNAs approximately 21 nt in length which play important roles in regulating gene expression in plants. Although many miRNA studies have focused on a few model plants, miRNAs and their target genes remain largely unknown in hot pepper (Capsicum annuum), one of the most important crops cultivated worldwide. Here, we employed high-throughput sequencing technology to identify miRNAs in pepper extensively from 10 different libraries, including leaf, stem, root, flower, and six developmental stage fruits. Based on a bioinformatics pipeline, we successfully identified 29 and 35 families of conserved and novel miRNAs, respectively. Northern blot analysis was used to validate further the expression of representative miRNAs and to analyze their tissue-specific or developmental stage-specific expression patterns. Moreover, we computationally predicted miRNA targets, many of which were experimentally confirmed using 5′ rapid amplification of cDNA ends analysis. One of the validated novel targets of miR-396 was a domain rearranged methyltransferase, the major de novo methylation enzyme, involved in RNA-directed DNA methylation in plants. This work provides the first reliable draft of the pepper miRNA transcriptome. It offers an expanded picture of pepper miRNAs in relation to other plants, providing a basis for understanding the functional roles of miRNAs in pepper.