Unique expression, processing regulation, and regulatory network of peach (Prunus persica) miRNAs

Harnessing the diversity of a lettuce wild relative to identify anthocyanin-related genes transcriptionally responsive to drought stress

Lettuce is a crop particularly vulnerable to drought. A transcriptomic study in the variety 'Romired' and the wild relative Lactuca homblei was conducted to understand the increase in anthocyanins (only significant in L. homblei) in response to drought previously observed. RNA-seq revealed more differentially expressed genes (DEGs), especially upregulated, in the wild species, in which the most abundant and significant GO terms were involved in regulatory processes (including response to water). Anthocyanin synthesis was triggered in L. homblei in response to drought, with 17 genes activated out of the 36 mapped in the phenylpropanoid-flavonoid pathway compared to 7 in 'Romired'. Nineteen candidate DEGs with the strongest change in expression and correlation with both anthocyanin content and drought were selected and validated by qPCR, all being differentially expressed only in the wild species with the two techniques. Their functions were related to anthocyanins and/or stress response and they harboured 404 and 11 polymorphisms in the wild and cultivated species, respectively. Some wild variants had high or moderate predicted impacts on the respective protein function: a transcription factor that responds to abiotic stresses, a heat shock protein involved in stomatal closure, and a phospholipase participating in anthocyanin accumulation under abiotic stress. These genetic variants could explain the differences in the gene expression patterns between the wild (significantly up/downregulated) and the cultivated (no significant changes) species. The diversity of this crop wild relative for anthocyanin-related genes involved in the response to drought could be exploited to improve lettuce resilience against some adverse climate effects.

Changes in MicroRNAs expression mediate molecular mechanism underlying the effect of MeJA on the biosynthesis of podophyllotoxin in Linum album cells

Podophyllotoxin (PTOX), produced by Linum album, is a monolignol that participates in plant defense strategies. Our previous study established that methyl jasmonate (MeJA) significantly stimulates PTOX production in L. album cells. However; the mechanisms by which MeJA regulates PTOX biosynthesis are uncovered. In the present study, we demonstrated that MeJA induces a time-dependent hydrogen peroxide (H2O2) and salicylic acid (SA) accumulation but reduces nitric oxide (NO) generation in L. album cells. PTOX biosynthetic genes such as PAL, CCR, CAD, and PLR were upregulated in response to MeJA exposure. Furthermore, the results of RT-qPCR revealed a positive correlation between the expression of PTOX biosynthetic genes and MeJA-induced upregulation of four miRNAs such as miR156, miR159, miR172, and miR396 at 12 h. Generally, this study revealed that MeJA mediates PTOX biosynthesis in L. album cells by inducing H2O2 and SA formation, which can probably upregulate the expression level of some miRNAs and biosynthetic genes in a redox balance-dependent manner.

euAP2a, a key gene that regulates flowering time in peach (Prunus persica) by modulating thermo-responsive transcription programming

Frequent spring frost damage threatens temperate fruit production, and breeding of late-flowering cultivars is an effective strategy for preventing such damage. However, this effort is often hampered by the lack of specific genes and markers and a lack of understanding of the mechanisms. We examined a Late-Flowering Peach (LFP) germplasm and found that its floral buds require a longer chilling period to release from their dormancy and a longer warming period to bloom than the control cultivar, two key characteristics associated with flowering time. We discovered that a 983-bp deletion in euAP2a, an APETALA2 (AP2)-related gene with known roles in regulating floral organ identity and flowering time, was primarily responsible for late flowering in LFP. This deletion disrupts an miR172 binding site, resulting in a gain-of-function mutation in euAP2a. Transcriptomic analyses revealed that at different stages of floral development, two chilling-responsive modules and four warm-responsive modules, comprising approximately 600 genes, were sequentially activated, forming a unique transcription programming. Furthermore, we found that euAP2a was transiently downregulated during the activation of these thermal-responsive modules at various stages. However, the loss of such transient, stage-specific downregulation of euAP2a caused by the deletion of miR172 binding sites resulted in the deactivation or delay of these modules in the LFP flower buds, suggesting that euAP2a acts as a transcription repressor to control floral developmental pace in peaches by modulating the thermo-responsive transcription programming. The findings shed light on the mechanisms behind late flowering in deciduous fruit trees, which is instrumental for breeding frost-tolerant cultivars.

Phased secondary small interfering RNAs in Camellia sinensis var. assamica

Phased secondary small interfering RNAs (phasiRNAs) in plants play important roles in regulating genome stability, plant development and stress adaption. Camellia sinensis var. assamica has immense economic, medicinal and cultural significance. However, there are still no studies of phasiRNAs and their putative functions in this valuable plant. We identified 476 and 43 PHAS loci which generated 4290 twenty one nucleotide (nt) and 264 twenty four nt phasiRNAs, respectively. Moreover, the analysis of degradome revealed more than 35000 potential targets for these phasiRNAs. We identified several conserved 21 nt phasiRNA generation pathways in tea plant, including miR390 → TAS3, miR482/miR2118 → NB-LRR, miR393 → F-box, miR828 → MYB/TAS4, and miR7122 → PPR in this study. Furthermore, we found that some transposase and plant mobile domain genes could generate phasiRNAs. Our results show that phasiRNAs target genes in the same family in cis- or trans-manners, and different members of the same gene family may generate the same phasiRNAs. The phasiRNAs, generated by transposase and plant mobile domain genes, and their targets, suggest that phasiRNAs may be involved in the inhibition of transposable elements in tea plant. To summarize, these results provide a comprehensive view of phasiRNAs in Camellia sinensis var. assamica.

De novo transcriptome profiling unveils the regulation of phenylpropanoid biosynthesis in unripe Piper nigrum berries

BACKGROUND

Black pepper (Piper nigrum L.) is rich in bioactive compounds that make it an imperative constituent in traditional medicines. Although the unripe fruits have long been used in different Ayurvedic formulations, the mechanism of gene regulation resulting in the production of the bioactive compounds in black pepper is not much investigated. Exploring the regulatory factors favouring the production of bioactive compounds ultimately help to accumulate the medicinally important content of black pepper. The factors that enhance the biosynthesis of these compounds could be potential candidates for metabolic engineering strategies to obtain a high level production of significant biomolecules.

RESULTS

Being a non-model plant, de novo sequencing technology was used to unravel comprehensive information about the genes and transcription factors that are expressed in mature unripe green berries of P. nigrum from which commercially available black pepper is prepared. In this study, the key gene regulations involved in the synthesis of bioactive principles in black pepper was brought out with a focus on the highly expressed phenylpropanoid pathway genes. Quantitative real-time PCR analysis of critical genes and transcription factors in the different developmental stages from bud to the mature green berries provides important information useful for choosing the developmental stage that would be best for the production of a particular bioactive compound. Comparison with a previous study has also been included to understand the relative position of the results obtained from this study.

CONCLUSIONS

The current study uncovered significant information regarding the gene expression and regulation responsible for the bioactivity of black pepper. The key transcription factors and enzymes analyzed in this study are promising targets for achieving a high level production of significant biomolecules through metabolic engineering.

Genome-Wide Identification of Membrane-Bound Fatty Acid Desaturase Genes in Three Peanut Species and Their Expression in Arachis hypogaea during Drought Stress

As a crop irrigated primarily by rain, the quality and yield of peanuts are significantly limited by drought. To date, many studies have indicated that fatty acid desaturase (FAD) genes enhance plant tolerance to drought stresses. In this study, 16, 15, and 31 FADs were identified in Arachis duranensis, Arachis ipaensis, and Arachis hypogaea, respectively. All the FADs were divided into four subfamilies, which had relatively conserved gene structures, motifs, and domains. The synteny relationships and chromosomal position analysis showed that the FADs in subgenome pairs, A. duranensis-A. hypogaea (AA) and A. ipaensis-A. hypogaea (BB), were homologous, and their physical locations were consistent. The Ka/Ks results indicated that nine FAD genes underwent a purifying selection, and Ah|FAD3.2 experienced positive selection during tetraploid peanut speciation. Various cis-acting elements related to hormone signaling and stress responsiveness in promoters and the predicted miRNA targeting Ah|FADs suggested that these genes play crucial roles in drought tolerance. The expression profiles of Ah|FADs in 22 tissues and drought-tolerant and -sensitive cultivars under drought stress suggested that 4 and 6 FADs were putative genes related to oil accumulation and drought, respectively. These findings will help provide insight into the potential functional roles of the FAD genes, which may aid in dealing with plant drought stress.

Small RNA Differential Expression Analysis Reveals miRNAs Involved in Dormancy Progression in Sweet Cherry Floral Buds

In sweet cherry (Prunus avium), as in other temperate woody perennials, bud dormancy allows for survival in adverse environmental conditions during winter. During this process, environmental signals such as short days and/or low temperatures trigger internal signals that enable buds to become tolerant to the cold. The process involves tracking chilling units up to chilling the requirement fulfillment to resume growth, a transition involving transcriptional regulation, metabolic signaling, and epigenetic-related regulatory events. Massive sequencing of small RNAs was performed to identify miRNAs involved in sweet cherry dormancy by comparing their expression in field (regular seasonal) and controlled non-stop (continuous) chilling conditions. miRNAs highlighted by sequencing were validated using specific stem-loop PCR quantification, confirming expression patterns for known miRNAs such as miR156e, miR166c, miR172d, miR391, miR482c, and miR535b, as well as for newly proposed miRNAs. In silico prediction of the target genes was used to construct miRNA/target gene nodes. In particular, the involvement of the sweet cherry version for the miR156/SQUAMOSA PROMOTER-BINDING-LIKE PROTEIN genes whose expression was opposite in the two conditions suggests their involvement on dormancy regulation in sweet cherry. miRNA levels indicate that the regulation of stress-related genes and hormone synthesis modulates the expression of calcium metabolism and cell development-associated genes. Understanding the regulatory networks involved in sweet cherry dormancy, particularly in the context of miRNA involvement, represents the first step in the development of new agricultural strategies that may help overcome the increasing challenges presented by global climate change.

Genome-Wide Analysis of Calmodulin Binding Transcription Activator (CAMTA) Gene Family in Peach (Prunus persica L. Batsch) and Ectopic Expression of PpCAMTA1 in Arabidopsis camta2,3 Mutant Restore Plant Development

Calmodulin-binding transcription activator (CAMTA) is a transcription factor family containing calmodulin (CaM) binding sites and is involved in plant development. Although CAMTAs in Arabidopsis have been extensively investigated, the functions of CAMTAs remain largely unclear in peaches. In this study, we identified five peach CAMTAs which contained conserved CG-1 box, ANK repeats, CaM binding domain (CaMBD) and IQ motifs. Overexpression in tobacco showed that PpCAMTA1/2/3 were located in the nucleus, while PpCAMTA4 and PpCAMTA5 were located in the plasma membrane. Increased expression levels were observed for PpCAMTA1 and PpCAMTA3 during peach fruit ripening. Expression of PpCAMTA1 was induced by cold treatment and was inhibited by ultraviolet B irradiation (UV-B). Driven by AtCAMTA3 promoter, PpCAMTA1/2/3 were overexpressed in Arabidopsis mutant. Here, we characterized peach PpCAMTA1, representing an ortholog of AtCAMTA3. PpCAMTA1 expression in Arabidopsis complements the developmental deficiencies of the camta2,3 mutant, and restored the plant size to the wild type level. Moreover, overexpressing PpCAMTA1 in camta2,3 mutant inhibited salicylic acid (SA) biosynthesis and expression of SA-related genes, resulting in a susceptibility phenotype to Pst DC3000. Taken together, our results provide new insights for CAMTAs in peach fruit and indicate that PpCAMTA1 is associated with response to stresses during development.

Less is more: natural variation disrupting a miR172 gene at the di locus underlies the recessive double-flower trait in peach (P. persica L. Batsch)

BACKGROUND

With the domestication of ornamental plants, artificial selective pressure favored the propagation of mutations affecting flower shape, and double-flower varieties are now readily available for many species. In peach two distinct loci control the double-flower phenotype: the dominant Di2 locus, regulated by the deletion of the binding site for miR172 in the euAP2 PETALOSA gene Prupe.6G242400, and the recessive di locus, of which the underlying factor is still unknown.

RESULTS

Based on its genomic location a candidate gene approach was used to identify genetic variants in a diverse panel of ornamental peach accessions and uncovered three independent mutations in Prupe.2G237700, the gene encoding the transcript for microRNA miR172d: a ~5.0 Kb LTR transposable element and a ~1.2 Kb insertion both positioned upstream of the sequence encoding the pre-miR172d within the transcribed region of Prupe.2G237700, and a ~9.5 Kb deletion encompassing the whole gene sequence. qRT-PCR analysis confirmed that expression of pre-miR172d was abolished in di/di genotypes homozygous for the three variants.

CONCLUSIONS

Collectively, PETALOSA and the mutations in micro-RNA miR172d identified in this work provide a comprehensive collection of the genetic determinants at the base of the double-flower trait in the peach germplasms.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

Metabolomic and transcriptomic analyses reveal the effects of self- and hetero-grafting on anthocyanin biosynthesis in grapevine

Grafting, which joins a scion from a cultivar with the stem of a rootstock from a grapevine wild relative, is commonly used in viticulture. Grafting has crucial effects on various phenotypes of the cultivar, including its phenology, biotic and abiotic resistance, berry metabolome, and coloration, but the underlying genetics and regulatory mechanisms are largely unexplored. In this study, we investigated the phenotypic, metabolomic, and transcriptomic profiles at three developmental stages (45, 75, and 105 days after flowering) of the Crimson Seedless cultivar (Vitis vinifera) grafted onto four rootstocks (three heterografts, CS/101-14, CS/SO4, and CS/110R and one self-graft, CS/CS) with own-rooted graft-free Crimson Seedless (CS) as the control. All the heterografts had a significant effect on berry reddening as early as ~45 days after flowering. The grafting of rootstocks promoted anthocyanin biosynthesis and accumulation in grape berries. The metabolomic features showed that cyanidin 3-O-glucoside, delphinidin 3-O-glucoside, malvidin 3-O-glucoside, peonidin 3-O-glucoside, and petunidin 3-O-glucoside were the pigments responsible for the purplish-red peel color. Transcriptomic analyses revealed that anthocyanin biosynthesis-related genes, from upstream (phenylalanine ammonia-lyase) to downstream (anthocyanidin 3-O-glucosyltransferase and anthocyanidin synthase), were upregulated with the accumulation of anthocyanins in the heterografted plants. At the same time, all these genes were also highly expressed and more anthocyanin was accumulated in self-grafted CS/CS samples compared with own-rooted graft-free CS samples, suggesting that self-grafting may also have promoted berry reddening in grapevine. Our results reveal global transcriptomic and metabolomic features in berry color regulation under different grafting conditions that may be useful for improving berry quality in viticulture.

Synthesis of flavour-related linalool is regulated by PpbHLH1 and associated with changes in DNA methylation during peach fruit ripening

Linalool is one of the common flavour-related volatiles across the plant kingdom and plays an essential role in determining consumer liking of plant foods. Although great process has been made in identifying terpene synthase (TPS) genes associated with linalool synthesis, much less is known about regulation of this pathway. We initiated study by identifying PpTPS3 encoding protein catalysing enantiomer (S)-(+)-linalool synthesis, which is a major linalool component (˜70%) observed in ripe peach fruit. Overexpression of PpTPS3 led to linalool accumulation, while virus-induced gene silencing of PpTPS3 led to a 66.5% reduction in linalool content in peach fruit. We next identified transcription factor (TF) PpbHLH1 directly binds to E-box (CACATG) in the PpTPS3 promoter and activates its expression based on yeast one-hybrid assay and EMSA analysis. Significantly positive correlation was also observed between PpbHLH1 expression and linalool production across peach cultivars. Peach fruit accumulated more linalool after overexpressing PpbHLH1 in peach fruit and reduced approximately 54.4% linalool production after silencing this TF. DNA methylation analysis showed increased PpTPS3 expression was associated with decreased 5 mC level in its promoter during peach fruit ripening, but no reverse pattern was observed for PpbHLH1. Arabidopsis and tomato fruits transgenic for peach PpbHLH1 synthesize and accumulate higher levels of linalool compared with wild-type controls. Taken together, these results would greatly facilitate efforts to enhance linalool production and thus improve flavour of fruits.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Epigenetics for Crop Improvement in Times of Global Change

Epigenetics has emerged as an important research field for crop improvement under the on-going climatic changes. Heritable epigenetic changes can arise independently of DNA sequence alterations and have been associated with altered gene expression and transmitted phenotypic variation. By modulating plant development and physiological responses to environmental conditions, epigenetic diversity-naturally, genetically, chemically, or environmentally induced-can help optimise crop traits in an era challenged by global climate change. Beyond DNA sequence variation, the epigenetic modifications may contribute to breeding by providing useful markers and allowing the use of epigenome diversity to predict plant performance and increase final crop production. Given the difficulties in transferring the knowledge of the epigenetic mechanisms from model plants to crops, various strategies have emerged. Among those strategies are modelling frameworks dedicated to predicting epigenetically controlled-adaptive traits, the use of epigenetics for in vitro regeneration to accelerate crop breeding, and changes of specific epigenetic marks that modulate gene expression of traits of interest. The key challenge that agriculture faces in the 21st century is to increase crop production by speeding up the breeding of resilient crop species. Therefore, epigenetics provides fundamental molecular information with potential direct applications in crop enhancement, tolerance, and adaptation within the context of climate change.

VvMYB114 mediated by miR828 negatively regulates trichome development of Arabidopsis

Trichome is a specialized structure differentiated during the morphogenesis of plant leaf epidermal cells. In recent years, with the continuous researches on trichome development of Arabidopsis and other plants, more and more genes related to trichome morphogenesis have been discovered, including R2R3-type MYB genes. In this study, we cloned a R2R3-type MYB family gene from grape, VvMYB114, a target gene of vvi-miR828. qRT-PCR showed that VvMYB114 mRNA accumulated during grape fruit ripening, and VvMYB114 protein had transcriptional activation activity. Heterologous overexpression of VvMYB114 in Arabidopsis reduced the number of trichome on leaves and stems. Mutating the miR828-binding site in VvMYB114 without altering amino-acid sequence had no effect on trichome development in Arabidopsis. The results showed a different role of the regulation of miR828 to VvMYB114 in Arabidopsis from in grape, which indicated the functional divergence of miRNA targeting homoeologous genes in different species played an important roles in evolution and useful trait selection.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Genome-Wide Changes of Regulatory Non-Coding RNAs Reveal Pollen Development Initiated at Ecodormancy in Peach

Bud dormancy is under the regulation of complex mechanisms including genetic and epigenetic factors. To study the function of regulatory non-coding RNAs in winter dormancy release, we analyzed the small RNA and long non-coding RNA (lncRNA) expression from peach (Prunus persica) floral buds in endodormancy, ecodormancy and bud break stages. Small RNAs underwent a major shift in expression primarily between dormancy and flowering with specific pairs of microRNAs and their mRNA target genes undergoing coordinated differential expression. From endodormancy to ecodormancy, ppe-miR6285 was significantly upregulated while its target gene, an ASPARAGINE-RICH PROTEIN involved in the regulation of abscisic acid signaling, was downregulated. At ecodormancy, ppe-miR2275, a homolog of meiosis-specific miR2275 across angiosperms, was significantly upregulated, supporting microsporogenesis in anthers at a late stage of dormancy. The expression of 785 lncRNAs, unlike the overall expression pattern in the small RNAs, demonstrated distinctive expression signatures across all dormancy and flowering stages. We predicted that a subset of lncRNAs were targets of microRNAs and found 18 lncRNA/microRNA target pairs with both differentially expressed across time points. The genome-wide differential expression and network analysis of non-coding RNAs and mRNAs from the same tissues provide new candidate loci for dormancy regulation and suggest complex noncoding RNA interactions control transcriptional regulation across these key developmental time points.

R2R3-MYB transcription factors, StmiR858 and sucrose mediate potato flavonol biosynthesis

Flavonols and other phenylpropanoids protect plants from biotic and abiotic stress and are dietarily desirable because of their health-promoting properties. The ability to develop new potatoes (Solanum tuberosum) with optimal types and amounts of phenylpropanoids is limited by lack of knowledge about the regulatory mechanisms. Exogenous sucrose increased flavonols, whereas overexpression of the MYB StAN1 induced sucrolytic gene expression. Heterologous StAN1 protein bound promoter fragments from sucrolytic genes (SUSY1 and INV1). Two additional MYBs and one microRNA were identified that regulated potato flavonols. Overexpression analysis showed MYB12A and C increased amounts of flavonols and other phenylpropanoids. Endogenous flavonol amounts in light-exposed organs were much higher those in the dark. Expression levels of StMYB12A and C were high in flowers but low in tubers. Transient overexpression of miR858 altered potato flavonol metabolism. Endogenous StmiR858 expression was much lower in flowers than leaves and correlated with flavonol amounts in these organs. Collectively, these findings support the hypothesis that sucrose, MYBs, and miRNA control potato phenylpropanoid metabolism in a finely tuned manner that includes a feedback loop between sucrose and StAN1. These findings will aid in the development of potatoes with phenylpropanoid profiles optimized for crop performance and human health.

Genome-wide identification, characterization and expression profile of glutaredoxin gene family in relation to fruit ripening and response to abiotic and biotic stresses in banana (Musa acuminata)

Glutaredoxins (GRXs) are disulfide oxidoreductases that are involved in various biological processes. However, little information on the role of GRXs in the regulation of fruit ripening and the response to stress is available. In this study, we isolated 64 GRX genes from banana genome. Their encoded GRX proteins could be classified into four classes: CC, CGFS, CPYC and GRL types. The distribution and synteny of these GRXs on chromosomes, the gene structures, the promoter sequences, and the possible protein subcellular localizations were characterized. Molecular interaction network analysis suggested that MaGRX might interact with glutathione reductase (GR), sulfiredoxin, peroxiredoxin (Prx), and NADPH-dependent thioredoxin reductase C (NTRC), contributing to the antioxidative defense of banana fruit. MicroRNA prediction showed that MaGRX genes might be targeted by different miRNAs. Transcriptome analysis characterized the expression profiles of different MaGRX genes during banana fruit ripening, and in response to different storage stresses. The results suggested that CC-type, CPYC-type and GRL-type MaGRXs might be more active than CGFS-type MaGRXs during banana fruit ripening and the response to stress. Moreover, MaGRX6/7/9/11/17/23/28 and MaGRL3/16/19 might play important roles in regulating fruit ripening or in response to low and high temperature, or Fusarium proliferatum infection.

The MicroRNA828/MYB12 Module Mediates Bicolor Pattern Development in Asiatic Hybrid Lily (Lilium spp.) Flowers

Some Asiatic hybrid lily cultivars develop bicolor tepals, which consist of anthocyanin-pigmented upper halves and un-pigmented lower halves. MYB12, a subgroup 6 member of R2R3-MYB that positively regulates anthocyanin biosynthesis, is downregulated in the lower halves. However, MYB12 is usually expressed over entire tepal regions in numerous lily cultivars. Why MYB12 of bicolor cultivars exhibits variable expression spatially in a single tepal remains unclear. Since the lily MYB12 mRNA harbored a binding site for microRNA828 (miR828), the involvement of miR828 in variable spatial accumulation of MYB12 transcripts was evaluated. We analyzed the cleavage of MYB12 mRNA, mature miR828 accumulation, and MYB12 transcript-derived siRNA generation (microRNA-seq). In the bicolor tepals, mature miR828 was more highly accumulated in the lower halves than in the upper halves, and miR828-directed cleavage of MYB12 transcripts was observed predominantly in the lower halves. Moreover, the cleavage triggered the production of secondary siRNA from MYB12 transcripts, and the siRNAs were accumulated predominantly in the lower halves. Consequently, miR828 suppressed MYB12 transcript accumulation in the white region, and the miR828/MYB12 module participated in the development of bicolor patterns in lily flowers. The results present the first example of a microRNA mediating flower color patterns. Finally, we discuss the potential of miR828 creating flower color variations through suppressing the activity of subgroup 6 R2R3-MYB positive regulators in other species.

Effects of exogenous abscisic acid on oil content, fatty acid composition, biodiesel properties and lipid components in developing Siberian apricot (Prunus sibirica) seeds

Previous studies in Siberian apricot (Prunus sibirica) seed kernel (SASK) have suggested the involvement of abscisic acid (ABA) signaling pathway in oil accumulation. However, there are few reports on the effects of ABA on the metabolism of fatty acids (FA) in seed development. Here, we first evaluated the response of developing SASK to ABA treatment, with a focus on oil content, FA composition, biodiesel properties, lipid compounds and gene expressions. Compared with control samples, the application of exogenous ABA increased the total oil content by 6.55% in mature SASK. The C18:1 content markedly increased in ABA treatment, and conversely C16:0 decreased. Exogenous ABA also improved the biodiesel properties of SASK oil, making it better suited to the specifications of biodiesel standards. Furthermore, the molecular species of phosphatidylcholine (PC), phosphatidic acid (PA), diacylglycerol (DAG) and triacylglycerol (TAG) were detected using lipidomics analysis. The 18:1/18:1 was the main component in PA, PC and DAG, while the main components of 18:1/18:1/18:2, 18:1/18:1/18:3, 18:2/18:2/18:2 and 18:1/18:1/18:1 in TAG. Most lipid species gradually increased with SASK maturity. In addition, the relative contents of TAG-18:1/18:1/18:2 and TAG-18:1/18:1/18:1 in developing SASK increased with the application of exogenous ABA. We also detected elevated gene expression of key genes involved in ABA chemical pathway, which likely affected FA biosynthesis and accumulation. Our results provide insight into the effects of ABA on the oil accumulation in developing SASK, which has direct applications to improving the quality of SASK-derived biodiesel.

Thermal-responsive genetic and epigenetic regulation of DAM cluster controlling dormancy and chilling requirement in peach floral buds

The Dormancy-associated MADS-box (DAM) gene cluster in peach serves as a key regulatory hub on which the seasonal temperatures act and orchestrate dormancy onset and exit, chilling response and floral bud developmental pace. Yet, how different temperature regimes interact with and regulate the six linked DAM genes remains unclear. Here, we demonstrate that chilling downregulates DAM1 and DAM3-6 in dormant floral buds with distinct patterns and identify DAM4 as the most abundantly expressed one. We reveal multiple epigenetic events, with tri-methyl histone H3 lysine 27 (H3K27me3) induced by chilling specifically in DAM1 and DAM5, a 21-nt sRNA in DAM3 and a ncRNA induced in DAM4. Such induction is inversely correlated with downregulation of their cognate DAMs. We also show that the six DAMs were hypermethylated, associating with the production of 24-nt sRNAs. Hence, the chilling-responsive dynamic of the different epigenetic elements and their interactions likely define distinct expression abundance and downregulation pattern of each DAM. We further show that the expression of the five DAMs remains steadily unchanged or continuously downregulated at the ensuing warm temperature after chilling, and this state of regulation correlates with robust increase of sRNA expression, H3K27me3 and CHH methylation, which is particularly pronounced in DAM4. Such robust increase of repressive epigenetic marks may irreversibly reinforce the chilling-imposed repression of DAMs to ensure flower-developmental programming free from any residual DAM inhibition. Taken together, we reveal novel information about genetic and epigenetic regulation of the DAM cluster in peach, which will be of fundamental significance in understanding of the regulatory mechanisms underlying chilling requirement and dormancy release, and of practical application for improvement of plasticity of flower time and bud break in fruit trees to adapt changing climates.

Transcriptomic and biochemical investigations support the role of rootstock-scion interaction in grapevine berry quality

Background

In viticulture, rootstock genotype plays a critical role to improve scion physiology, berry quality and to adapt grapevine (Vitis vinifera L.) to different environmental conditions. This study aimed at investigating the effect of two different rootstocks (1103 Paulsen - P - and Mgt 101–14 - M) in comparison with not grafted plants - NGC - on transcriptome (RNA-seq and small RNA-seq) and chemical composition of berry skin in Pinot noir, and exploring the influence of rootstock-scion interaction on grape quality. Berry samples, collected at veraison and maturity, were investigated at transcriptional and biochemical levels to depict the impact of rootstock on berry maturation.

Results

RNA- and miRNA-seq analyses highlighted that, at veraison, the transcriptomes of the berry skin are extremely similar, while variations associated with the different rootstocks become evident at maturity, suggesting a greater diversification at transcriptional level towards the end of the ripening process. In the experimental design, resembling standard agronomic growth conditions, the vines grafted on the two different rootstocks do not show a high degree of diversity. In general, the few genes differentially expressed at veraison were linked to photosynthesis, putatively because of a ripening delay in not grafted vines, while at maturity the differentially expressed genes were mainly involved in the synthesis and transport of phenylpropanoids (e.g. flavonoids), cell wall loosening, and stress response. These results were supported by some differences in berry phenolic composition detected between grafted and not grafted plants, in particular in resveratrol derivatives accumulation.

Conclusions

Transcriptomic and biochemical data demonstrate a stronger impact of 1103 Paulsen rootstock than Mgt 101–14 or not grafted plants on ripening processes related to the secondary metabolite accumulations in berry skin tissue. Interestingly, the MYB14 gene, involved in the feedback regulation of resveratrol biosynthesis was up-regulated in 1103 Paulsen thus supporting a putative greater accumulation of stilbenes in mature berries.

BrmiR828 Targets BrPAP1, BrMYB82, and BrTAS4 Involved in the Light Induced Anthocyanin Biosynthetic Pathway in Brassica rapa

Comprehensive research in various plants shows that the metabolic pathway of anthocyanin biosynthesis is affected by environmental factors and regulated by microRNAs through post-transcriptional regulation. In seedlings of Brassica rapa Tsuda, the accumulation of anthocyanin is induced by light. However, the roles of BrmiR828 in the light-induced synthesis of anthocyanin in Brassica rapa remain to be explored. Here, a primary transcript of BrmiR828 was identified to be located on the chromosomes of the A03 sub-genome. Five candidate MYB family genes were predicted as targets of BrmiR828 in the database of Brassica rapa (BRAD, V1.1) by using psRNATarget. The transcript abundance of mature BrmiR828 was reduced in seedlings of Brassica rapa Tsuda under blue light irradiation comparing with dark treatment. However, Real-time PCR showed the transcript level of the five candidate targets, Bra004162, Bra022602, Bra001917, Bra029113, and Bra039763 was up-regulated when the seedlings exposed to blue or UV-A light. Trans-acting siRNA gene 4 (BrTAS4) was also identified to have a higher transcript level under blue and UV-A light irradiation than that in dark treatment. RNA ligase mediated 5′amplification of cDNA ends (RLM-5′ RACE) showed that BrmiR828 can splice the mRNA of Bra039763, Bra022602, and BrTAS4 on binding sites. Phylogenetic analysis of candidate BrMYBs targets along with MYBs from Arabidopsis thaliana showed that Bra039763, Bra004162, Bra001917, Bra029113, and Bra022602 are classified to the same group with AtMYB75, AtMYB114, AtMYB90, AtMYB113, and AtMYB82 which are involved in the anthocyanin biosynthetic pathway. As a result, light-induced down-regulation of BrmiR828 can target BrTAS4, BrPAP1 (Bra039763), MYB82 (Bra022602) to negatively regulate their transcript levels leading to the accumulation of MYB transcription factors that positively regulate anthocyanin biosynthesis in light-exposed seedlings of Brassica rapa.

Mapping a double flower phenotype-associated gene DcAP2L in Dianthus chinensis

The double flower is a highly important breeding trait that affects the ornamental value in many flowering plants. To get a better understanding of the genetic mechanism of double flower formation in Dianthus chinensis, we have constructed a high-density genetic map using 140 F2 progenies derived from a cross between a single flower genotype and a double flower genotype. The linkage map was constructed using double-digest restriction site-associated DNA sequencing (ddRAD-seq) with 2353 single nucleotide polymorphisms (SNPs). Quantitative trait locus (QTL) mapping analysis was conducted for 12 horticultural traits, and major QTLs were identified for nine of the 12 traits. Among them, two major QTLs accounted for 20.7% and 78.1% of the total petal number variation, respectively. Bulked segregant RNA-seq (BSR-seq) was performed to search accurately for candidate genes associated with the double flower trait. Integrative analysis of QTL mapping and BSR-seq analysis using the reference genome of Dianthus caryophyllus suggested that an SNP mutation in the miR172 cleavage site of the A-class flower organ identity gene APETALA2 (DcAP2L) is responsible for double flower formation in Dianthus through regulating the expression of DcAG genes.

Isolation of microglia-derived extracellular vesicles: towards miRNA signatures and neuroprotection

The functional preservation of the central nervous system (CNS) is based on the neuronal plasticity and survival. In this context, the neuroinflammatory state plays a key role and involves the microglial cells, the CNS-resident macrophages. In order to better understand the microglial contribution to the neuroprotection, microglia-derived extracellular vesicles (EVs) were isolated and molecularly characterized to be then studied in neurite outgrowth assays. The EVs, mainly composed of exosomes and microparticles, are an important cell-to-cell communication process as they exhibit different types of mediators (proteins, lipids, nucleic acids) to recipient cells. The medicinal leech CNS was initially used as an interesting model of microglia/neuron crosstalk due to their easy collection for primary cultures. After the microglia-derived EV isolation following successive methods, we developed their large-scale and non-targeted proteomic analysis to (i) detect as many EV protein markers as possible, (ii) better understand the biologically active proteins in EVs and (iii) evaluate the resulting protein signatures in EV-activated neurons. The EV functional properties were also evaluated in neurite outgrowth assays on rat primary neurons and the RNAseq analysis of the microglia-derived EVs was performed to propose the most representative miRNAs in microglia-derived EVs. This strategy allowed validating the EV isolation, identify major biological pathways in EVs and corroborate the regenerative process in EV-activated neurons. In parallel, six different miRNAs were originally identified in microglia-derived EVs including 3 which were only known in plants until now. The analysis of the neuronal proteins under the microglial EV activation suggested possible miRNA-dependent regulation mechanisms. Taken together, this combination of methodologies showed the leech microglial EVs as neuroprotective cargos across species and contributed to propose original EV-associated miRNAs whose functions will have to be evaluated in the EV-dependent dialog between microglia and neurons.

miR828 and miR858 regulate VvMYB114 to promote anthocyanin and flavonol accumulation in grapes

MicroRNAs are a class of non-coding small RNAs involved in the negative regulation of gene expression, which play critical roles in developmental and metabolic pathways. Studies in several plants have identified a few microRNAs and other small RNAs that target regulators of the phenylpropanoid metabolic pathway called the MYB transcription factors. However, it is not well understood how sRNA-mediated regulation of MYBs influences the accumulation of specific secondary metabolites. Using sRNA sequencing, degradome analysis, mRNA sequencing, and proteomic analysis, we establish that grape lines with high anthocyanin content express two MYB-targeting microRNAs abundantly, resulting in the differential expression of specific MYB proteins. miR828 and miR858 target coding sequences of specific helix motifs in the mRNA sequences of MYB proteins. Targeting by miR828 caused MYB RNA decay and the production of a cascade of secondary siRNAs that depend on RNA-dependent RNA polymerase 6. MYB suppression and cascade silencing was more robust in grape lines with high anthocyanin content than in a flavonol-rich grape line. We establish that microRNA-mediated silencing targeted the repressor class of MYBs to promote anthocyanin biosynthesis in grape lines with high anthocyanins. We propose that this process regulates the expression of appropriate MYBs in grape lines to produce specific secondary metabolites.

RNAi-Mediated Resistance Against Viruses in Perennial Fruit Plants

Small RNAs (sRNAs) are 20–30-nucleotide-long, regulatory, noncoding RNAs that induce silencing of target genes at the transcriptional and posttranscriptional levels. They are key components for cellular functions during plant development, hormone signaling, and stress responses. Generated from the cleavage of double-stranded RNAs (dsRNAs) or RNAs with hairpin structures by Dicer-like proteins (DCLs), they are loaded onto Argonaute (AGO) protein complexes to induce gene silencing of their complementary targets by promoting messenger RNA (mRNA) cleavage or degradation, translation inhibition, DNA methylation, and/or histone modifications. This mechanism of regulating RNA activity, collectively referred to as RNA interference (RNAi), which is an evolutionarily conserved process in eukaryotes. Plant RNAi pathways play a fundamental role in plant immunity against viruses and have been exploited via genetic engineering to control disease. Plant viruses of RNA origin that contain double-stranded RNA are targeted by the RNA-silencing machinery to produce virus-derived small RNAs (vsRNAs). Some vsRNAs serve as an effector to repress host immunity by capturing host RNAi pathways. High-throughput sequencing (HTS) strategies have been used to identify endogenous sRNA profiles, the “sRNAome”, and analyze expression in various perennial plants. Therefore, the review examines the current knowledge of sRNAs in perennial plants and fruits, describes the development and implementation of RNA interference (RNAi) in providing resistance against economically important viruses, and explores sRNA targets that are important in regulating a variety of biological processes.

Identification of miRNAs and their target genes in Larix olgensis and verified of differential expression miRNAs

BACKGROUND

MiRNAs (microRNA) are 18-24 nt endogenous noncoding RNAs that regulate gene expression at the post-transcriptional level, including tissue-specific, developmental timing and evolutionary conservation gene expression.

RESULTS

This study used high-throughput sequencing technology for the first time in Larix olgensis, predicted 78 miRNAs, including 12,229,003 reads sRNA, screened differentially expressed miRNAs. Predicting target genes was helpful for understanding the miRNA regulation function and obtained 333 corresponding target genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation were analysed, mostly including nucleic acid binding, plant hormone signal transduction, pantothenate and CoA biosynthesis, and cellulose synthase. This study will lay the foundation for clarifying the complex miRNA-mediated regulatory network for growth and development. In view of this, spatio-temporal expression of miR396, miR950, miR164, miR166 and miR160 were analysed in Larix olgensis during the growth stages of not lignified, beginning of lignification, and completely lignified in different tissues (root, stem, and leaf) by quantitative real-time PCR (qRT-PCR). There were differences in the expression of miRNAs in roots, stems and leaves in the same growth period. At 60 days, miR160, miR166 and miR396-2 exhibited the highest expression in leaves. At 120 days, most miRNAs in roots and stems decreased significantly. At 180 days, miRNAs were abundantly expressed in roots and stems. Meanwhile, analysis of the expression of miRNAs in leaves revealed that miR396-2 was reduced as time went on, whereas other miRNAs increased initially and then decreased. On the other hand, in the stems, miR166-1 was increase, whereas other miRNAs, especially miR160, miR164, miR396 and miR950-1, first decreased and then increased. Similarly, in the roots, miR950-2 first decreased and then increased, whereas other miRNAs exhibited a trend of continuous increase.

CONCLUSIONS

The present investigation included rapid isolation and identification of miRNAs in Larix olgensis through construction of a sRNA library using Solexa and predicted 78 novel miRNAs, which showed differential expression levels in different tissues and stages. These results provided a theoretical basis for further revealing the genetic regulation mechanism of miRNA in the growth and development of conifers and the verification of function in target genes.

Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosarugosa Thunb

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles by regulating other genes. Rosa rugosa Thunb. is an important ornamental and edible plant, yet there are only a few studies on the miRNAs and their functions in R. rugosa.

RESULTS

We sequenced 10 samll RNA profiles from the roots, petals, pollens, stamens, and leaves and 4 RNA-seq profiles in leaves and petals to analysis miRNA, phasiRNAs and mRNAs in R. rugosa. In addition, we acquired a degradome sequencing profile from leaf of R. rugosa to identify miRNA and phasiRNA targets using the SeqTar algorithm. We have identified 321 conserved miRNA homologs including primary transcripts for 25 conserved miRNAs, and 22 novel miRNAs. We identified 592 putative targets of the conserved miRNAs or tasiRNAs that showed significant accumulations of degradome reads. We found differential expression patterns of conserved miRNAs in five different tissues of R. rugosa. We identified three hundred and thirty nine 21 nucleotide (nt) PHAS loci, and forty nine 24 nt PHAS loci, respectively. Our results suggest that miR482 triggers generations of phasiRNAs by targeting nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes in R. rugosa. Our results also suggest that the deregulated genes in leaves and petals are significantly enriched in GO terms and KEGG pathways related to metabolic processes and photosynthesis.

CONCLUSIONS

These results significantly enhanced our knowledge of the miRNAs and phasiRNAs, as well as their potential functions, in R. rugosa.

The Role of UV-B light on Small RNA Activity During Grapevine Berry Development

We explored the effects of ultraviolet B radiation (UV-B) on the developmental dynamics of microRNAs and phased small-interfering-RNA (phasi-RNAs)-producing loci by sequencing small RNAs in vegetative and reproductive organs of grapevine (Vitis vinifera L.). In particular, we tested different UV-B conditions in in vitro-grown plantlets (high-fluence exposition) and in berries from field-grown (radiation filtering) and greenhouse-grown (low- and high-fluence expositions) adult plants throughout fruit development and ripening. The functional significance of the observed UV-coordinated miRNA responses was supported by degradome evidences of ARGONAUTE (AGO)-programmed slicing of mRNAs. Co-expression patterns of the up-regulated miRNAs miR156, miR482, miR530, and miR828 with cognate target gene expressions in response to high-fluence UV-B was tested by q-RT-PCR. The observed UV-response relationships were also interrogated against two published UV-stress and developmental transcriptome datasets. Together, the dynamics observed between miRNAs and targets suggest that changes in target abundance are mediated transcriptionally and, in some cases, modulated post-transcriptionally by miRNAs. Despite the major changes in target abundance are being controlled primarily by those developmental effects that are similar between treatments, we show evidence for novel miRNA-regulatory networks in grape. A model is proposed where high-fluence UV-B increases miR168 and miR530 that target ARGONAUTE 1 (AGO1) and a Plus-3 domain mRNA, respectively, while decreasing miR403 that targets AGO2, thereby coordinating post-transcriptional gene silencing activities by different AGOs. Up-regulation of miR3627/4376 could facilitate anthocyanin accumulation by antagonizing a calcium effector, whereas miR395 and miR399, induced by micronutrient deficiencies known to trigger anthocyanin accumulation, respond positively to UV-B radiation. Finally, increases in the abundance of an anthocyanin-regulatory MYB-bHLH-WD40 complex elucidated in Arabidopsis, mediated by UV-B-induced changes in miR156/miR535, could contribute to the observed up-regulation of miR828. In turn, miR828 would regulate the AtMYB113-ortologues MYBA5, A6 and A7 (and thereby anthocyanins) via a widely conserved and previously validated auto-regulatory loop involving miR828 and phasi TAS4abc RNAs.

Banana sRNAome and degradome identify microRNAs functioning in differential responses to temperature stress

Background

Temperature stress is a major environmental factor affecting not only plant growth and development, but also fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play important roles in various biological processes. Harvested banana fruit can exhibit distinct symptoms in response to different temperature stresses, but the underlying miRNA-mediated regulatory mechanisms remained unknown.

Results

Here, we profiled temperature-responsive miRNAs in banana, using deep sequencing and computational and molecular analyses. In total 113 known miRNAs and 26 novel banana-specific miRNAs were identified. Of these miRNAs, 42 miRNAs were expressed differentially under cold and heat stresses. Degradome sequencing identified 60 target genes regulated by known miRNAs and half of these targets were regulated by 15 temperature-responsive miRNAs. The correlative expression patterns between several miRNAs and their target genes were further validated via qRT-PCR. Our data showed that miR535 and miR156 families may derive from a common ancestor during evolution and jointly play a role in fine-tuning SPL gene expression in banana. We also identified the miRNA-triggered phased secondary siRNAs in banana and found miR393-TIR1/AFB phasiRNA production displaying cold stress-specific enrichment.

Conclusions

Our results provide a foundation for understanding the miRNA-dependent temperature stress response in banana. The characterized correlations between miRNAs and their response to temperature stress could serve as markers in the breeding programs or tools for improving temperature tolerance of banana.

Electronic supplementary material

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

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Cis-directed cleavage and nonstoichiometric abundances of 21-nucleotide reproductive phased small interfering RNAs in grasses

Post-transcriptional gene silencing in plants results from independent activities of diverse small RNA types. In anthers of grasses, hundreds of loci yield noncoding RNAs that are processed into 21- and 24-nucleotide (nt) phased small interfering RNAs (phasiRNAs); these are triggered by miR2118 and miR2275. We characterized these 'reproductive phasiRNAs' from rice (Oryza sativa) panicles and anthers across seven developmental stages. Our computational analysis identified characteristics of the 21-nt reproductive phasiRNAs that impact their biogenesis, stability, and potential functions. We demonstrate that 21-nt reproductive phasiRNAs can function in cis to target their own precursors. We observed evidence of this cis regulatory activity in both rice and maize (Zea mays). We validated this activity with evidence of cleavage and a resulting shift in the pattern of phasiRNA production. We characterize biases in phasiRNA biogenesis, demonstrating that the Pol II-derived 'top' strand phasiRNAs are consistently higher in abundance than the bottom strand. The first phasiRNA from each precursor overlaps the miR2118 target site, and this impacts phasiRNA accumulation or stability, evident in the weak accumulation of this phasiRNA position. Additional influences on this first phasiRNA duplex include the sequence composition and length, and we show that these factors impact Argonaute loading.

Deletion of the miR172 target site in a TOE-type gene is a strong candidate variant for dominant double-flower trait in Rosaceae

Double flowers with supernumerary petals have been selected by humans for their attractive appearance and commercial value in several ornamental plants, including Prunus persica (peach), a recognized model for Rosaceae genetics and genomics. Despite the relevance of this trait, knowledge of the underlying genes is limited. Of two distinct loci controlling the double-flower phenotype in peach, we focused on the dominant Di2 locus. High-resolution linkage mapping in five segregating progenies delimited Di2 to an interval spanning 150 858 bp and 22 genes, including Prupe.6G242400 encoding an euAP2 transcription factor. Analyzing genomic resequencing data from single- and double-flower accessions, we identified a deletion spanning the binding site for miR172 in Prupe.6G242400 as a candidate variant for the double-flower trait, and we showed transcript expression for both wild-type and deleted alleles. Consistent with the proposed role in controlling petal number, Prupe.6G242400 is expressed in buds at critical times for floral development. The indelDi2 molecular marker designed on this sequence variant co-segregated with the phenotype in 621 progenies, accounting for the dominant inheritance of the Di2 locus. Further corroborating the results in peach, we identified a distinct but similar mutation in the ortholog of Prupe.6G242400 in double-flower roses. Phylogenetic analysis showed that these two genes belong to a TARGET OF EAT (TOE)-type clade not represented in Arabidopsis, indicating a divergence of gene functions between AP2-type and TOE-type factors in Arabidopsis and other species. The identification of orthologous candidate genes for the double-flower phenotype in two important Rosaceae species provides valuable information to understand the genetic control of this trait in other major ornamental plants.

Genome-wide identification and classification of MYB superfamily genes in peach

The MYB transcription factor superfamily is one of the largest superfamilies modulating various biological processes in plants. Over the past few decades, many MYB superfamily genes have been identified and characterized in some plant species. However, genes belonging to the MYB superfamily in peach (Prunus persica) have not been comprehensively identified and characterized although the genome sequences of peach were released several years ago. In total, this study yielded a set of 256 MYB superfamily genes that was divided into five subfamilies: the R2R3-MYB (2R-MYB), R1R2R3-MYB (3R-MYB), MYB-related (1R-MYB), 4R-MYB, and Atypical-MYB subfamilies. These subfamilies contained 128, 4, 109, 1, and 14 members, respectively. The 128 R2R3-MYB subfamily genes in peach were further clustered into 35 groups, and the 109 MYB-related subfamily genes were further clustered into 6 groups: the CCA1-like, CPC-like, TBP-like, I-box-binding-like, R-R-type, and Peach-specific groups. The motif compositions and exon/intron structures within each group within the R2R3-MYB or MYB-related subfamily in peach were highly conserved. The logo sequences of the R2 and R3 repeats of R2R3-MYB subfamily members were highly conserved with those in these repeats of several other plant species. Except for 48 novel peach-specific MYB genes, the remaining 208 out of 256 MYB genes in peach were conserved with the corresponding 198 MYB genes in A. thaliana. Additionally, the 256 MYB genes unevenly distributed on chromosomes 1 to 8 of the peach genome. Eighty-one orthologous pairs of peach/A. thaliana MYB genes were identified among 256 MYB genes in peach and 198 MYB genes in A. thaliana in this study. In addition, 146 pairs of paralogous MYB genes were identified on the eight chromosomes of peach. The expression levels of some of the 51 MYB genes selected for qRT-PCR analysis decreased or increased with red-fleshed fruit development, while the expression patterns of some genes followed no clear rules over the five developmental stages of fruits. This study laid the foundation for further functional analysis of MYB superfamily genes in peach and enriched the knowledge of MYB superfamily genes in plant species.

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.

Genome-wide analysis of basic helix-loop-helix superfamily members in peach

The basic helix-loop-helix (bHLH) transcription factor family is a superfamily found in all eukaryotes that plays important roles in regulating growth and development. Over the past several decades, many bHLH superfamily genes have been identified and characterized in herbaceous and woody plants. However, the genes belonging to the bHLH superfamily in peach (Prunus persica) have not yet been comprehensively identified and characterized. Here, we identified 95 members of the bHLH superfamily in the peach genome, and these genes were classified into 19 subfamilies based on a phylogenetic comparison with bHLH proteins from Arabidopsis. The members within each subfamily were highly conserved according to the analysis of motif compositions and exon/intron organizations. The 95 bHLH genes were unevenly distributed on chromosomes 1 to 8 of the peach genome. We identified 57 pairs of bHLH members that were orthologous between peach and Arabidopsis. Additionally, 48 pairs of paralogous bHLH genes were identified on the eight chromosomes of the peach genome. Coupled with relative expression analysis of bHLH genes in red-fleshed peach fruit at five developmental stages, we identified several bHLH genes that might be involved in fruit development and anthocyanin biosynthesis. This study provides insight into the molecular mechanisms through which these genes are involved in the regulation of biological and biochemical processes in peach and lays the foundation for further studies on these genes.

Coupling of microRNA-directed phased small interfering RNA generation from long noncoding genes with alternative splicing and alternative polyadenylation in small RNA-mediated gene silencing

MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play vital regulatory roles in plant growth and development. Little is known about these small RNAs in litchi (Litchi chinensis), an economically important fruit crop widely cultivated in Southeast Asia. We profiled the litchi small RNA population with various deep-sequencing techniques and in-depth bioinformatic analyses. The genome-wide identification of miRNAs, their target genes, and phasiRNA-generating (PHAS) genes/loci showed that the function of miR482/2118 has expanded, relative to its canonical function. We also discovered that, for 29 PHAS loci, miRNA-mediated phasiRNA production was coupled with alternative splicing (AS) and alternative polyadenylation (APA). Most of these loci encoded long noncoding RNAs. An miR482/2118 targeted locus gave rise to four main transcript isoforms through AS/APA, and diverse phasiRNAs generated from these isoforms appeared to target long terminal repeat (LTR) retrotransposons and other unrelated genes. This coupling enables phasiRNA production from different exons of noncoding PHAS genes and yields diverse phasiRNA populations, both broadening and altering the range of downstream phasiRNA-regulated genes. Our results reveal the diversity of miRNA and phasiRNA in litchi, and demonstrate AS/APA as a new layer of regulation in small RNA-mediated gene silencing.

Small RNA pathways responsible for non-cell-autonomous regulation of plant reproduction

In angiosperms, germline precursors and germ cells are always attached to or engulfed within somatic companion cells until just before fertilization. This is because sperm and egg cells develop as part of the multicellular gametophyte. Thus, the non-cell-autonomous regulation by somatic companions plays important roles in efficient reproduction, in addition to the cell-autonomous regulation. Epigenetic silencing of transposable elements is one of the central events by which the germline transmits the error-free genome to the next generation. This review focuses on small RNA-mediated epigenetic regulation of meiosis, spore formation and pollen development. Besides microRNA (miRNA) and small interfering RNA (siRNA), animals express PIWI-interacting RNA (piRNA), a germline-specific class of small RNAs. Plants lack piRNA-like RNAs and, instead, express unique classes of small RNAs: trans-acting siRNA (tasiRNA) and phased secondary siRNA (phasiRNA). Especially in grass species, 21- and 24-nucleotide phasiRNAs are abundant in anthers during premeiosis and meiosis. This review also describes recent progress in reproductive phasiRNA research.

Phased secondary small interfering RNAs in Panaxnotoginseng

Background

Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species.

Methods

Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5’-RACE to validate some of the identified phasiRNA targets.

Results

After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5’-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets.

Conclusions

We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Electronic supplementary material

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

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

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.

Identification and characterization of Prunus persica miRNAs in response to UVB radiation in greenhouse through high-throughput sequencing

Background

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression of target mRNAs involved in plant growth, development, and abiotic stress. As one of the most important model plants, peach (Prunus persica) has high agricultural significance and nutritional values. It is well adapted to be cultivated in greenhouse in which some auxiliary conditions like temperature, humidity, and UVB etc. are needed to ensure the fruit quality. However, little is known about the genomic information of P. persica under UVB supplement. Transcriptome and expression profiling data for this species are therefore important resources to better understand the biological mechanism of seed development, formation and plant adaptation to environmental change. Using a high-throughput miRNA sequencing, followed by qRT-PCR tests and physiological properties determination, we identified the responsive-miRNAs under low-dose UVB treatment and described the expression pattern and putative function of related miRNAs and target genes in chlorophyll and carbohydrate metabolism.

Results

A total of 164 known peach miRNAs belonging to 59 miRNA families and 109 putative novel miRNAs were identified. Some of these miRNAs were highly conserved in at least four other plant species. In total, 1794 and 1983 target genes for known and novel miRNAs were predicted, respectively. The differential expression profiles of miRNAs between the control and UVB-supplement group showed that UVB-responsive miRNAs were mainly involved in carbohydrate metabolism and signal transduction. UVB supplement stimulated peach to synthesize more chlorophyll and sugars, which was verified by qRT-PCR tests of related target genes and metabolites’ content measurement.

Conclusion

The high-throughput sequencing data provided the most comprehensive miRNAs resource available for peach study. Our results identified a series of differentially expressed miRNAs/target genes that were predicted to be low-dose UVB-responsive. The correlation between transcriptional profiles and metabolites contents in UVB supplement groups gave novel clues for the regulatory mechanism of miRNAs in Prunus. Low-dose UVB supplement could increase the chlorophyll and sugar (sorbitol) contents via miRNA-target genes and therefore improve the fruit quality in protected cultivation of peaches.

Electronic supplementary material

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

Global Identification, Classification, and Expression Analysis of MAPKKK genes: Functional Characterization of MdRaf5 Reveals Evolution and Drought-Responsive Profile in Apple

Mitogen-activated protein kinase kinase kinases (MAPKKKs) are pivotal components of Mitogen-activated protein kinase (MAPK) cascades, which play a significant role in many biological processes. Although genome-wide analysis of MAPKKKs has been conducted in many species, extant results in apple are scarce. In this study, a total of 72 putative MdMAPKKKs in Raf-like group, 11 in ZIK-like group and 37 in MEEK were identified in apple firstly. Predicted MdMAPKKKs were located in 17 chromosomes with diverse densities, and there was a high-level of conservation in and among the evolutionary groups. Encouragingly, transcripts of 12 selected MdMAPKKKs were expressed in at least one of the tested tissues, indicating that MdMAPKKKs might participate in various physiological and developmental processes in apple. Moreover, they were found to respond to drought stress in roots and leaves, which suggested a possible conserved response to drought stress in different species. Overexpression of MdRaf5 resulted in a hyposensitivity to drought stress, which was at least partially due to the regulation of stomatal closure and transpiration rates. To the best of our knowledge, this is the first genome-wide functional analysis of the MdMAPKKK genes in apple, and it provides valuable information for understanding MdMAPKKKs signals and their putative functions.

miR3954 is a trigger of phasiRNAs that affects flowering time in citrus

In plant, a few 22-nt miRNAs direct cleavages of their targets and trigger the biogenesis of phased small interfering RNAs (phasiRNAs) in plant. In this study, we characterized a miRNA triggering phasiRNAs generation, miR3954, and explored its downstream target genes and potential function. Our results demonstrated that miR3954 showed specific expression in the flowers of citrus species, and it targeted a NAC transcription factor (Cs7 g22460) and two non-coding RNA transcripts (lncRNAs, Cs1 g09600 and Cs1 g09635). The production of phasiRNAs was detected from transcripts targeted by miR3954, and was further verified in both sequencing data and transient expression experiments. PhasiRNAs derived from the two lncRNAs targeted not only miR3954-targeted NAC gene but also additional NAC homologous genes. No homologous genes of these two lncRNAs were found in plants other than citrus species, implying that this miR3954-lncRNAs-phasiRNAs-NAC pathway is likely citrus-specific. Transgenic analysis indicated that the miR3954-overexpressing lines showed decreased transcripts of lncRNA, elevated abundance of phasiRNAs and reduced expression of NAC genes. Interestingly, the overexpression of miR3954 leads to early flowering in citrus plants. In summary, our results illustrated a model of the regulatory network of miR3954-lncRNA-phasiRNAs-NAC, which may be functionally involved in flowering in citrus.

Genome-wide identification and characterization of phased small interfering RNA genes in response to Botrytis cinerea infection in Solanum lycopersicum

Phased small interfering RNAs (phasiRNAs) are encoded by a novel class of genes known as phasiRNA producing (PHAS) genes. These genes play important regulatory roles by targeting protein coding transcripts in plant species. In this study, 91 regions were identified as potential PHAS loci in tomato, with additional evidence that seven of them can be triggered by five miRNAs. Among the identified loci, 51 were located in genic regions, and the remaining 40 were located in intergenic regions. The transient overexpression of PHAS15 and PHAS26 demonstrated that phasiRNAs predicted by PhaseTank were indeed generated from their respective PHAS loci. Using sRNA-seq data from B. cinerea-infected tomato leaves, we identified 50 B. cinerea-responsive phasiRNAs with increased abundance and five with decreased abundance. Moreover, 164 targets of these differentially expressed phasiRNAs were predicted, and 94 of them were confirmed experimentally using degradome data. Gene ontology analysis of the targets revealed an enrichment of genes with functions related to defense responses and signaling regulation. These results suggest that a large number of endogenous siRNAs, such as phasiRNAs, have not yet been identified in tomato and underscore the urgent need to systematically identify and functionally analyze siRNAs in tomato.

The legume miR1514a modulates a NAC transcription factor transcript to trigger phasiRNA formation in response to drought

Recent studies have identified microRNAs as post-transcriptional regulators involved in stress responses. miR1514a is a legume microRNA that is induced in response to drought stress in Phaseolus vulgaris (common bean) and shows differential accumulation levels in roots during water deficit in two cultivars with different drought tolerance phenotypes. A recent degradome analysis revealed that miR1514a targets the transcripts of two NAC transcription factors (TFs), Phvul.010g121000 and Phvul.010g120700. Furthermore, expression studies and small RNA-seq data indicate that only Phvul.010g120700 generates phasiRNAs, which also accumulate under water deficit conditions. To confirm these results, we over-expressed miR1514a in transgenic hairy roots, and observed a reduced accumulation of Phvul.010g120700 and an increase in NAC-derived phasiRNAs; inhibition of miR1514a activity resulted in the opposite effect. Moreover, we determined that a NAC-derived phasiRNA associates with ARGONAUTE 1 (AGO1), suggesting that it is functional. In addition, a transcriptome analysis of transgenic hairy roots with reduced miR1514a levels revealed several differentially expressed transcripts, mainly involved in metabolism and stress responses, suggesting they are regulated by the NAC TF and/or by phasiRNAs. This work therefore demonstrates the participation of miR1514 in the regulation of a NAC transcription factor transcript through phasiRNA production during the plant response to water deficit.