Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis

miR394 modulates brassinosteroid signaling to regulate hypocotyl elongation in Arabidopsis

The interplay between microRNAs (miRNAs) and phytohormones allows plants to integrate multiple internal and external signals to optimize their survival of different environmental conditions. Here, we report that miR394 and its target gene LEAF CURLING RESPONSIVENESS (LCR), which are transcriptionally responsive to BR, participate in BR signaling to regulate hypocotyl elongation in Arabidopsis thaliana. Phenotypic analysis of various transgenic and mutant lines revealed that miR394 negatively regulates BR signaling during hypocotyl elongation, whereas LCR positively regulates this process. Genetically, miR394 functions upstream of BRASSINOSTEROID INSENSITIVE2 (BIN2), BRASSINAZOLEs RESISTANT1 (BZR1), and BRI1-EMS-SUPPRESSOR1 (BES1), but interacts with BRASSINOSTEROID INSENSITIVE1 (BRI1) and BRI1 SUPRESSOR PROTEIN (BSU1). RNA-sequencing analysis suggested that miR394 inhibits BR signaling through BIN2, as miR394 regulates a significant number of genes in common with BIN2. Additionally, miR394 increases the accumulation of BIN2 but decreases the accumulation of BZR1 and BES1, which are phosphorylated by BIN2. MiR394 also represses the transcription of PACLOBUTRAZOL RESISTANCE1/5/6 and EXPANSIN8, key genes that regulate hypocotyl elongation and are targets of BZR1/BES1. These findings reveal a new role for a miRNA in BR signaling in Arabidopsis.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Dynamic Regulation of Grapevine's microRNAs in Response to Mycorrhizal Symbiosis and High Temperature

MicroRNAs (miRNAs) are non-coding small RNAs that play crucial roles in plant development and stress responses and can regulate plant interactions with beneficial soil microorganisms such as arbuscular mycorrhizal fungi (AMF). To determine if root inoculation with distinct AMF species affected miRNA expression in grapevines subjected to high temperatures, RNA-seq was conducted in leaves of grapevines inoculated with either Rhizoglomus irregulare or Funneliformis mosseae and exposed to a high-temperature treatment (HTT) of 40 °C for 4 h per day for one week. Our results showed that mycorrhizal inoculation resulted in a better plant physiological response to HTT. Amongst the 195 identified miRNAs, 83 were considered isomiRs, suggesting that isomiRs can be biologically functional in plants. The number of differentially expressed miRNAs between temperatures was higher in mycorrhizal (28) than in non-inoculated plants (17). Several miR396 family members, which target homeobox-leucine zipper proteins, were only upregulated by HTT in mycorrhizal plants. Predicted targets of HTT-induced miRNAs in mycorrhizal plants queried to STRING DB formed networks for Cox complex, and growth and stress-related transcription factors such as SQUAMOSA promoter-binding-like-proteins, homeobox-leucine zipper proteins and auxin receptors. A further cluster related to DNA polymerase was found in R. irregulare inoculated plants. The results presented herein provide new insights into miRNA regulation in mycorrhizal grapevines under heat stress and can be the basis for functional studies of plant-AMF-stress interactions.

Integrated High-Throughput Sequencing, Microarray Hybridization and Degradome Analysis Uncovers MicroRNA-Mediated Resistance Responses of Maize to Pathogen Curvularia lunata

Curvularia lunata (Wakker) Boed, the causal agent of leaf spot in maize, is prone to mutation, making it difficult to control. RNAi technology has proven to be an important tool of genetic engineering and functional genomics aimed for crop improvement. MicroRNAs (miRNAs), which act as post-transcriptional regulators, often cause translational repression and gene silencing. In this article, four small RNA (sRNA) libraries were generated from two maize genotypes inoculated by C. lunata; among these, ltR1 and ltR2 were from the susceptible variety Huangzao 4 (HZ), ltR3 and ltR4, from the resistant variety Luyuan (LY), and 2286, 2145, 1556 and 2504 reads were annotated as miRNA in these four sRNA libraries, respectively. Through the combined analysis of high-throughput sequencing, microarray hybridization and degradome, 48 miRNAs were identified as being related to maize resistance to C. lunata. Among these, PC-732 and PC-169, two new maize miRNAs discovered, were predicted to cleave mRNAs of metacaspase 1 (AMC1) and thioredoxin family protein (Trx), respectively, possibly playing crucial roles in the resistance of maize to C. lunata. To further confirm the role of PC-732 in the interaction of maize and C. lunata, the miRNA was silenced through STTM (short tandem target mimic) technology, and we found that knocking down PC-732 decreased the susceptibility of maize to C. lunata. Precisely speaking, the target gene of PC-732 might inhibit the expression of disease resistance-related genes during the interaction between maize and C. lunata. Overall, the findings of this study indicated the existence of miRNAs involved in the resistance of maize to C. lunata and will contribute to rapidly clarify the resistant mechanism of maize to C. lunata.

Cultivar-specific miRNA-mediated RNA silencing in grapes

In-depth comparative degradome analysis of two domesticated grape cultivars with diverse secondary metabolite accumulation reveals differential miRNA-mediated targeting. Small (s)RNAs such as micro(mi)RNAs and secondary small interfering (si) often work as negative switches of gene expression. In plants, it is well known that miRNAs target and cleave mRNAs that have high sequence complementarity. However, it is not known if there are variations in miRNA-mediated targeting between subspecies and cultivars that have been subjected to vast genetic modifications through breeding and other selections. Here, we have used PAREsnip2 tool for analysis of degradome datasets derived from two contrasting domesticated grape cultivars having varied fruit color, habit and leaf shape. We identified several interesting variations in sRNA targeting using degradome and 5'RACE analysis between two contrasting grape cultivars that was further correlated using RNA-seq analysis. Several of the differences we identified are associated with secondary metabolic pathways. We propose possible means by which sRNAs might contribute to diversity in secondary metabolites and other development pathways between two domesticated cultivars of grapes.

MicroRNA miR394 regulates flowering time in Arabidopsis thaliana

miR394 regulates Arabidopsis flowering time in a LCR-independent manner. Arabidopsis plants harboring mutations in theMIR394 genes exhibit early flowering, lower expression of floral repressor FLC and higher expression of floral integrators FT and SOC1. Plant development occurs throughout its entire life cycle and involves a phase transition between vegetative and reproductive phases, leading to the flowering process, fruit formation and ultimately seed production. It has been shown that the microRNA394 (miR394) regulates the accumulation of the transcript coding for LEAF CURLING RESPONSIVENESS, a member of a family of F-Box proteins. The miR394 pathway regulates several processes including leaf morphology and development of the shoot apical meristem during embryogenesis, as well as having been assigned a role in the response to biotic and abiotic stress in Arabidopsis thaliana and other species. Here, we characterized plants harboring mutations in MIR394 precursor genes and demonstrate that mir394a mir394b double mutants display an early flowering phenotype which correlates with a lower expression of FLOWERING LOCUS C earlier in development and higher expression of the floral integrators FLOWERING LOCUS T and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1. Consequently, mutant plants produce fewer branches and exhibit lower seed production. Our work reveals previously unknown developmental aspects regulated by the miR394 pathway, in an LCR-independent manner, contributing to the characterization of the multiple roles of this versatile plant regulatory miRNA.

Exploring the Effect of Methyl Jasmonate on the Expression of microRNAs Involved in Biosynthesis of Active Compounds of Rosemary Cell Suspension Cultures through RNA-Sequencing

Our aim in the experiment was to study the effects of methyl jasmonates (MeJA) on the active compounds of rosemary suspension cells, the metabolites' change of contents under different concentrations of MeJA, including 0 (CK), 10 (M10), 50 (M50) and 100 μM MeJA (M100). The results demonstrated that MeJA treatments promoted the accumulation of rosmarinic acid (RA), carnosic acid (CA), flavonoids, jasmonate (JA), gibberellin (GA), and auxin (IAA); but reduced the accumulations of abscisic acid (ABA), salicylic acid (SA), and aspartate (Asp). In addition, 50 and 100 μM MeJA promoted the accumulation of alanine (Ala) and glutamate (Glu), and 50 μM MeJA promoted the accumulation of linoleic acid and alpha-linolenic acid in rosemary suspension cells. Comparative RNA-sequencing analysis of different concentrations of MeJA showed that a total of 30, 61, and 39 miRNAs were differentially expressed in the comparisons of CKvsM10, CKvsM50, CKvsM100, respectively. The analysis of the target genes of the differentially expressed miRNAs showed that plant hormone signal transduction, linoleic acid, and alpha-linolenic acid metabolism-related genes were significantly enriched. In addition, we found that miR160a-5p target ARF, miR171d_1 and miR171f_3 target DELLA, miR171b-3p target ETR, and miR156a target BRI1, which played a key role in rosemary suspension cells under MeJA treatments. qRT-PCR of 12 differentially expressed miRNAs and their target genes showed a high correlation between the RNA-seq and the qRT-PCR result. Amplification culture of rosemary suspension cells in a 5 L stirred bioreactor showed that cell biomass accumulation in the bioreactor was less than that in the shake flask under the same conditions, and the whole cultivation period was extended to 14 d. Taken together, MeJA promoted the synthesis of the active compounds in rosemary suspension cells in a wide concentration range via concentration-dependent differential expression patterns. This study provided an overall view of the miRNAs responding to MeJA in rosemary.

Characterization and regulatory mechanism analysis of VvmiR156a-VvAGL80 pair during grapevine flowering and parthenocarpy process induced by gibberellin

MicroRNA156 (miR156) is an important conserved miRNA family in plants. Recently, we revealed VvmiR156a could involve in the modulation of gibberellin (GA)-mediated flower and berry development process of grapevine (Vitis vinifera L.). However, how to manipulate this process is unclear. For this, we used the GA-induced grapevine parthenocarpy system to investigate the regulatory roles of VvmiR156a during this process. Here, we cloned the mature and precursor sequences of VvmiR156a in Wink grape and identified its potential target gene VvAGL80, which belongs to the MADS-box gene family. Moreover, using RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-RACE) and poly(A)polymerase-mediated 3' rapid amplification of cDNA (PPM-RACE) technologies, it confirmed that VvAGL80 was the true target gene of VvmiR159a. Analysis of promoter cis-elements and β-glucuronidase (GUS) staining showed that both VvmiR156a and VvAGL80 contained GA-responsive elements and could respond to GA treatments. Quantitative real-time-polymerase chain reaction (qRT-PCR) analysis exhibited the VvmiR156a and VvAGL80 showed opposite expression trends during grapevine flower and berry development, indicating that VvmiR156a negatively regulated the expression of VvAGL80 during this process. After GA treatment, the expression of miR156 in flowers was downregulated significantly, while that of VvAGL80 was upregulated, thereby accelerating grapevine flowering. Furthermore, GA treatment enhanced the negative regulation of VvmiR156a on VvAGL80 in seed, especially at the seed-coat hardening stage, which was the key period of seed growth and development. Our findings enriched the knowledge of the regulatory mechanism of the miRNA-mediated grapevine parthenocarpy process.

Integrated Four Comparative-Omics Reveals the Mechanism of the Terpenoid Biosynthesis in Two Different Overwintering Cryptomeria fortunei Phenotypes

Chinese cedar (Cryptomeria fortunei) is a tree species with important ornamental, medicinal, and economic value. Terpenoids extracted from the essential oil of C. fortunei needles have been considered valuable ingredients in the pharmaceutical and cosmetic industries. However, the possible gene regulation mechanisms that limit terpenoid biosynthesis in this genus are poorly understood. Here, we adopted integrated metabolome analysis, transcriptome, small-RNA (sRNA), and degradome sequencing to analyze the differences in terpenoid regulatory mechanisms in two different overwintering C. fortunei phenotypes (wild-type and an evergreen mutant). A total of 1447/6219 differentially synthesized metabolites (DSMs)/unigenes (DEGs) were detected through metabolome/transcriptome analyses, and these DSMs/DEGs were significantly enriched in flavonoid and diterpenoid biosynthesis pathways. In C. fortunei needles, 587 microRNAs (miRNAs), including 67 differentially expressed miRNAs (DERs), were detected. Among them, 8346 targets of 571 miRNAs were predicted using degradome data, and a 72-miRNA-target regulatory network involved in the metabolism of terpenoids and polyketides was constructed. Forty-one targets were further confirmed to be involved in terpenoid backbone and diterpenoid biosynthesis, and target analyses revealed that two miRNAs (i.e., aly-miR168a-5p and aof-miR396a) may be related to the different phenotypes and to differential regulation of diterpenoid biosynthesis. Overall, these results reveal that C. fortunei plants with the evergreen mutation maintain high terpenoid levels in winter through miRNA-target regulation, which provides a valuable resource for essential oil-related bioengineering research.

Movement of small RNAs in and between plants and fungi

RNA interference is a biological process whereby small RNAs inhibit gene expression through neutralizing targeted mRNA molecules. This process is conserved in eukaryotes. Here, recent work regarding the mechanisms of how small RNAs move within and between organisms is examined. Small RNAs can move locally and systemically in plants through plasmodesmata and phloem, respectively. In fungi, transportation of small RNAs may also be achieved by septal pores and vesicles. Recent evidence also supports bidirectional cross-kingdom communication of small RNAs between host plants and adapted fungal pathogens to affect the outcome of infection. We discuss several mechanisms for small RNA trafficking and describe evidence for transport through naked form, combined with RNA-binding proteins or enclosed by vesicles.

MicroRNAs in Vitis vinifera cv. Chardonnay Are Differentially Expressed in Response to Diaporthe Species

Diaporthe species are important pathogens, saprobes, and endophytes on grapevines. Several species are known, either as agents of pre- or post-harvest infections, as causal agents of many relevant diseases, including swelling arm, trunk cankers, leaf spots, root and fruit rots, wilts, and cane bleaching. A growing body of evidence exists that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation, during plant development and responses to biotic and abiotic stresses. In this study, we explored differentially expressed miRNAs in response to Diaporthe eres and Diaporthe bohemiae infection in Vitis vinifera cv. Chardonnay under in vitro conditions. We used computational methods to predict putative miRNA targets in order to explore the involvement of possible pathogen response pathways. We identified 136 known and 41 new miRNA sequence variants, likely generated through post-transcriptional modifications. In the Diaporthe eres treatment, 61 known and 17 new miRNAs were identified while in the Diaporthe bohemiae treatment, 101 known and 21 new miRNAs were revealed. Our results contribute to further understanding the role miRNAs play during plant pathogenesis, which is possibly crucial in understanding disease symptom development in grapevines infected by D. eres and D. bohemiae.

Characterization of genome-wide microRNAs and their roles in development and biotic stress in pear

Using a genome-wide analysis of miRNAs in 'Yali' pear (Pyrus bretschneideri) via the next-generation high-throughput sequencing of small RNAs with a bioinformatics analysis, we found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 and their target genes function in viral defense in 'Duli' and 'Hongbaoshi'. pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319 and their targets function in auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Successful fruit production in pear (Pyrus spp.) depends on the use of optimal combinations of rootstocks and scions. Deciphering plant-pathogen defense mechanisms and hormone signaling pathways is an important step towards developing pear rootstocks and varieties with improved qualities. In the current study, we combined next-generation sequencing of small RNAs with a bioinformatics analysis to systematically identify and characterize 298 miRNAs in the pear scion cultivar 'Yali' (Pyrus bretschneideri). We also analyzed miRNAs in three rootstock varieties ('Duli', 'Zhongai 4', and 'Zhongai 5') and one scion cultivar ('Hongbaoshi'). We found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 are induced following infection with the pear virus Apple stem pitting virus (ASPV), and identified their target genes (pbRPS6, pbNAC, pbTLR, and pbRX-CC, respectively), which participate in viral defense pathways in 'Duli' and 'Hongbaoshi'. Furthermore, we identified pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319, and found that the production of these miRNAs was suppressed under low levels of synthetic auxin. The targets of these miRNAs (pbARF, pbAEC, pbSCL, and pbTCP4) respond to auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Our results lay the foundation for breeding improved pear cultivars.

Selection of reference genes for miRNA qRT-PCR under abiotic stress in grapevine

Grapevine is among the fruit crops with high economic value, and because of the economic losses caused by abiotic stresses, the stress resistance of Vitis vinifera has become an increasingly important research area. Among the mechanisms responding to environmental stresses, the role of miRNA has received much attention recently. qRT-PCR is a powerful method for miRNA quantitation, but the accuracy of the method strongly depends on the appropriate reference genes. To determine the most suitable reference genes for grapevine miRNA qRT-PCR, 15 genes were chosen as candidate reference genes. After eliminating 6 candidate reference genes with unsatisfactory amplification efficiency, the expression stability of the remaining candidate reference genes under salinity, cold and drought was analysed using four algorithms, geNorm, NormFinder, deltaCt and Bestkeeper. The results indicated that U6 snRNA was the most suitable reference gene under salinity and cold stresses; whereas miR168 was the best for drought stress. The best reference gene sets for salinity, cold and drought stresses were miR160e + miR164a, miR160e + miR168 and ACT + UBQ + GAPDH, respectively. The selected reference genes or gene sets were verified using miR319 or miR408 as the target gene.

Identification of Submergence-Responsive MicroRNAs and Their Targets Reveals Complex MiRNA-Mediated Regulatory Networks in Lotus (Nelumbo nucifera Gaertn)

MicroRNAs (miRNAs) are endogenous non-coding RNAs with important regulatory functions in plant development and stress responses. However, their population abundance in lotus (Nelumbo nucifera Gaertn) has so far been poorly described, particularly in response to stresses. In this work, submergence-related miRNAs and their target genes were systematically identified, compared, and validated at the transcriptome-wide level using high-throughput sequencing data of small RNA, Mrna, and the degradome. A total of 128 known and 20 novel miRNAs were differentially expressed upon submergence. We identified 629 target transcripts for these submergence-responsive miRNAs. Based on the miRNA expression profiles and GO and KEGG annotation of miRNA target genes, we suggest possible molecular responses and physiological changes of lotus in response to submergence. Several metabolic, physiological and morphological adaptations-related miRNAs, i.e., NNU_far-miR159, NNU_gma-miR393h, and NNU_aly-miR319c-3p, were found to play important regulatory roles in lotus response to submergence. This work will contribute to a better understanding of miRNA-regulated adaption responses of lotus to submergence stress.

Identification of ARGONAUTE/Small RNA Cleavage Sites by Degradome Sequencing

The method described here enables the high-throughput identification of cleaved mRNA targets of ARGONAUTE/small RNA complexes. The protocol is based on a modified 5'-rapid amplification of cDNA ends combined with deep sequencing of 3' cleavage products of mRNAs. Poly(A) RNA is purified from the tissue of interest which is followed by a 5'-RNA adapter ligation. The ligated products are then reverse transcribed, amplified, and digested with MmeI. After gel separation, a 3' double-stranded DNA adapter is ligated to the fragments, which are then amplified and index labeled for the high-throughput sequencing of pooled degradome libraries. Sequencing datasets from pooled libraries can be analyzed with different bioinformatic approaches.

MicroRNA-mediated regulation of gene expression in the response of rice plants to fungal elicitors

MicroRNAs (miRNAs) are small non-coding RNAs that have important regulatory functions in plant growth, development, and response to abiotic stress. Increasing evidence also supports that plant miRNAs contribute to immune responses to pathogens. Here, we used deep sequencing of small RNA libraries for global identification of rice miRNAs that are regulated by fungal elicitors. We also describe 9 previously uncharacterized miRNAs in rice. Combined small RNA and degradome analyses revealed regulatory networks enriched in elicitor-regulated miRNAs supported by the identification of their corresponding target genes. Specifically, we identified an important number of miRNA/target gene pairs involved in small RNA pathways, including miRNA, heterochromatic and trans-acting siRNA pathways. We present evidence for miRNA/target gene pairs implicated in hormone signaling and cross-talk among hormone pathways having great potential in regulating rice immunity. Furthermore, we describe miRNA-mediated regulation of Conserved-Peptide upstream Open Reading Frame (CPuORF)-containing genes in rice, which suggests the existence of a novel regulatory network that integrates miRNA and CPuORF functions in plants. The knowledge gained in this study will help in understanding the underlying regulatory mechanisms of miRNAs in rice immunity and develop appropriate strategies for rice protection.

Identification of microRNAs Involved in Regeneration of the Secondary Vascular System in Populus tomentosa Carr

Wood formation is a complex developmental process primarily controlled by a regulatory transcription network. MicroRNAs (miRNAs) can modulate the expression of target genes involved in plant growth and development by inducing mRNA degradation and translational repression. In this study, we used a model of secondary vascular system regeneration established in Populus tomentosa to harvest differentiating xylem tissues over time for high-throughput sequencing of small RNAs. Analysis of the sequencing data identified 209 known and 187 novel miRNAs during this regeneration process. Degradome sequencing analysis was then performed, revealing 157 and 75 genes targeted by 21 known and 30 novel miRNA families, respectively. Gene ontology enrichment of these target genes revealed that the targets of 15 miRNAs were enriched in the auxin signaling pathway, cell differentiation, meristem development, and pattern specification process. The major biological events during regeneration of the secondary vascular system included the sequential stages of vascular cambium initiation, formation, and differentiation stages in sequence. This study provides the basis for further analysis of these miRNAs to gain greater insight into their regulatory roles in wood development in trees.

Comparative Profiling of miRNAs and Target Gene Identification in Distant-Grafting between Tomato and Lycium (Goji Berry)

Local translocation of small RNAs between cells is proved. Long distance translocation between rootstock and scion is also well documented in the homo-grafting system, but the process in distant-grafting is widely unexplored where rootstock and scion belonging to different genera. Micro RNAs are a class of small, endogenous, noncoding, gene silencing RNAs that regulate target genes of a wide range of important biological pathways in plants. In this study, tomato was grafted onto goji (Lycium chinense Mill.) to reveal the insight of miRNAs regulation and expression patterns within a distant-grafting system. Goji is an important traditional Chinese medicinal plant with enriched phytochemicals. Illumina sequencing technology has identified 68 evolutionary known miRNAs of 37 miRNA families. Moreover, 168 putative novel miRNAs were also identified. Compared with control tomato, 43 (11 known and 32 novels) and 163 (33 known and 130 novels) miRNAs were expressed significantly different in shoot and fruit of grafted tomato, respectively. The fruiting stage was identified as the most responsive in the distant-grafting approach and 123 miRNAs were found as up-regulating in the grafted fruit which is remarkably higher compare to the grafted shoot tip (28). Potential targets of differentially expressed miRNAs were found to be involved in diverse metabolic and regulatory pathways. ADP binding activities, molybdopterin synthase complex and RNA helicase activity were found as enriched terms in GO (Gene Ontology) analysis. Additionally, "metabolic pathways" was revealed as the most significant pathway in KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis. The information of the small RNA transcriptomes that are obtained from this study might be the first miRNAs elucidation for a distant-grafting system, particularly between goji and tomato. The results from this study will provide the insights into the molecular aspects of miRNA-mediated regulation in the medicinal plant goji, and in grafted tomato. Noteworthy, it would provide a basis how miRNA signals could exchange between rootstock and scion, and the relevance to diverse biological processes.

MicroRNA Regulation of Abiotic Stress Response in 7B-1 Male-Sterile Tomato Mutant

The 7B-1 tomato (Solanum lycopersicum L. 'Rutgers') is a male-sterile mutant with enhanced tolerance to abiotic stress in a blue-light (BL) specific manner compared with its wild-type (WT). This makes the 7B-1 a potential candidate for hybrid seed breeding and stress engineering. To identify small RNAs (sRNAs) linked to stress tolerance of 7B-1, two sRNA libraries from BL-grown 7B-1 and WT seedlings treated simultaneously with abscisic acid (ABA) and mannitol were sequenced, and sRNA profiles were compared. Twenty nine families of known microRNAs (miRNAs) and 27 putative novel miRNAs were identified from the two libraries. MiR5300, miR5301, miR2916, and a novel miRNA denoted miR#C were upregulated, while miR159, miR166, miR472, miR482, and two novel miRNAs, miR#A and miR#D, were downregulated in stress-treated 7B-1 seedlings. MiRNA targets with potential roles in stress regulation were validated by rapid amplification of 5' complementary DNA ends (5'-RACE) analysis. Expression of miR159, miR166, miR472, miR482, miR#A, and miR#D together with their targets were further investigated in response to ABA, mannitol, NaCl, and cold treatments and a strong negative correlation was observed between the levels of these miRNAs and expression of their targets. Only miR159 and miR166 responded to cold treatment. MiR#A and its target were regulated by ABA and mannitol as early as 0.5 h after the treatments, while other miRNAs and targets were regulated only after 2 h. This suggests a role in early response to stress for miR#A. Our data suggests that miR159, miR166, miR472, miR482, miR#A, and miR#D are likely to facilitate the BL-specific enhanced tolerance of 7B-1 to abiotic stress.

Characterization of Expression Patterns of Grapevine MicroRNA Family Members using MicroRNA Rapid Amplification of Complementary DNA Ends

Grapevine (Vitis vinifera L.), with important nutritional values and health benefits, is one of the most economically fruit crop worldwide. In the present study, real-time quantitative polymerase chain reaction (qRT-PCR) and microRNA rapid amplification of cDNA ends (miR-RACE) techniques were used to characterize the expression and diversification patterns of various grapevine microRNAs (Vv-miRNAs) and their family members in grapevine. Based on our results, eight different grapevine miRNAs (miR159, miR164, miR167, miR172, miR319, miR393, miR396, and miR398) and their family members were expressed in different tissues at various developmental stages. The qRT-PCR results showed that the expression levels of Vv-miRNAs during grapevine development were dynamic. Furthermore, based on miR-RACE analysis and polymerase chain reaction (PCR) product sequencing results, different members within the same miRNA family were also expressed at different levels. Comparing the spatiotemporal expression levels of different members in the same miRNA family indicated that some miRNA families might have a key miRNA member that played the prominent role in regulation of their subsequent common target genes. In conclusion, our results showed that miR-RACE is a powerful technique to analyze the expression patterns of different members in the same miRNA family in terms of reverse-transcription (RT) efficiency and specificity. The findings of the expression diversification among Vv-miRNA family members and the existence of some Vv-miRNAs playing the key role could add to our understanding about the regulatory role of miRNAs in grapevine.

Genome-Wide Mapping and Analysis of Grapevine MicroRNAs and Their Potential Target Genes

MicroRNAs (miRNAs) are single-stranded, nonprotein-coding, endogenously expressed, small RNAs 19 to 25 nucleotides in length. Recognizing the lack of specific and systematic studies on genome-wide mapping of grapevine (Vitis vinifera L.) miRNAs, we conducted genome-wide mapping of Vv-miRNAs (V. vinifera miRNAs), SB-miRNAs (V. vinifera L. 'Summer Black' miRNAs), and Va-miRNAs (V. amurensis Rupr. miRNAs). The mapping results revealed that many of miRNAs located within the intergenic region had independent transcription units. To further validate the mapping results and existence of miRNAs, 12 randomly selected precursors of miRNAs (pre-miRNAs) were successfully cloned and sequenced. Subsequently, 15 conserved and 29 nonconserved intragenic (intronic, exonic) Vv-miRNA genes, 24 nonconserved intragenic SB-miRNA genes, and 23 nonconserved intragenic Va-miRNA genes were labeled on the basis of their locations in host genes, and 15 MIRNA clusters were detected. Interestingly, five miRNA pairs, namely, Vv-MIR395b and Vv-MIR395c, Vv-MIR482 and Vv-MIRC13, Vv-MIR172a and Va-MIR057, SB-MIR024 and Vv-MIRC35, and Vv-MIRC36 and Va-MIR073 were clustered in the host genes GSVIVT01011558001, GSVIVT01008132001, GSVIVT01031524001, GSVIVT01028156001, and GSVIVT01024516001, respectively. To validate the existence of target genes and miRNA-guided cleavage sites, 3'-end product of four predicted target messenger RNAs were amplified by RNA ligase-mediated 5' rapid amplification of cDNA ends. In addition, we also conducted contrastive analysis on the genomic location of miRNAs and their potential target genes. Results showed that the order of priority of miRNA-target interaction may be less closely related with their genomic location. These findings could benefit some further study on grapevine functional genomics and will provide new insights into the regulatory mechanisms and evolution of miRNAs in Vitis species.

High-Throughput Sequencing Reveals Diverse Sets of Conserved, Nonconserved, and Species-Specific miRNAs in Jute

MicroRNAs play a pivotal role in regulating a broad range of biological processes, acting by cleaving mRNAs or by translational repression. A group of plant microRNAs are evolutionarily conserved; however, others are expressed in a species-specific manner. Jute is an agroeconomically important fibre crop; nonetheless, no practical information is available for microRNAs in jute to date. In this study, Illumina sequencing revealed a total of 227 known microRNAs and 17 potential novel microRNA candidates in jute, of which 164 belong to 23 conserved families and the remaining 63 belong to 58 nonconserved families. Among a total of 81 identified microRNA families, 116 potential target genes were predicted for 39 families and 11 targets were predicted for 4 among the 17 identified novel microRNAs. For understanding better the functions of microRNAs, target genes were analyzed by Gene Ontology and their pathways illustrated by KEGG pathway analyses. The presence of microRNAs identified in jute was validated by stem-loop RT-PCR followed by end point PCR and qPCR for randomly selected 20 known and novel microRNAs. This study exhaustively identifies microRNAs and their target genes in jute which will ultimately pave the way for understanding their role in this crop and other crops.

Identification of cold-inducible microRNAs in grapevine

Low temperature is one of the most important environmental factors that limits the geographical distribution and productivity of grapevine. However, the molecular mechanisms on how grapevine responds to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although miRNAs and their targets have been identified in several Vitis species, their participation during cold accumulation in grapevine remains unknown. In this study, two small RNA libraries were generated from micropropagated 'Muscat Hamburg' (V. vinifera) plantlets under normal and low temperatures (4°C). A total of 163 known miRNAs and 67 putative novel miRNAs were detected from two small RNA libraries by Solexa sequencing. Forty-four cold-inducible miRNAs were identified through differentially expressed miRNAs (DEMs) analysis; among which, 13 belonged to upregulated DEMs while 31 belonged downregulated DEMs. The expression patterns of the 13 DEMs were verified by real-time RT-PCR analysis. The prediction of the target genes for DEMs indicated that miRNA may regulate transcription factors, including AP2, SBP, MYB, bHLH, GRAS, and bZIP under cold stress. The 5'-RLM RACE were conducted to verify the cleavage site of predicted targets. Seven predicted target genes for four known and three novel vvi-miRNAs showed specific cleavage sites corresponding to their miRNA complementary sequences. The expression pattern of these seven target genes revealed negative correlation with the expression level of the corresponding vvi-miRNAs. Our results indicated that a diverse set of miRNAs in V. vinifera are cold-inducible and may play an important role in cold stress response.

Discovery of Novel Leaf Rust Responsive microRNAs in Wheat and Prediction of Their Target Genes

MicroRNAs are endogenous small noncoding RNAs which play critical roles in gene regulation. Few wheat (Triticum aestivum L.) miRNA sequences are available in miRBase repertoire and knowledge of their biological functions related to biotic stress is limited. We identified 52 miRNAs, belonging to 19 families, from next-generation transcriptome sequence data based on homology search. One wheat specific novel miRNA was identified but could not be ascribed or assigned to any known miRNA family. Differentially expressed 22 miRNAs were found between susceptible and resistant wheat near-isogenic lines inoculated with leaf rust pathogen Puccinia triticina and compared with mock inoculated controls. Most miRNAs were more upregulated in susceptible NIL compared to resistant NIL. We identified 1306 potential target genes for these 52 miRNAs with vital roles in response to stimuli, signaling, and diverse metabolic and cellular processes. Gene ontology analysis showed 66, 20, and 35 target genes to be categorized into biological process, molecular function, and cellular component, respectively. A miRNA-mediated regulatory network revealed relationships among the components of the targetome. The present study provides insight into potential miRNAs with probable roles in leaf rust pathogenesis and their target genes in wheat which establish a foundation for future studies.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

Comparative profiling of miRNA expression in developing seeds of high linoleic and high oleic safflower (Carthamus tinctorius L.) plants

Vegetable oils high in oleic acid are considered to be advantageous because of their better nutritional value and potential industrial applications. The oleic acid content in the classic safflower oil is normally 10-15% while a natural mutant (ol) accumulates elevated oleic acid up to 70% in seed oil. As a part of our investigation into the molecular features of the high oleic (HO) trait in safflower we have profiled the microRNA (miRNA) populations in developing safflower seeds expressing the ol allele in comparison to the wild type high linoleic (HL) safflower using deep sequencing technology. The small RNA populations of the mid-maturity developing embryos of homozygous ol HO and wild type HL safflower had a very similar size distribution pattern, however, only ~16.5% of the unique small RNAs were overlapping in these two genotypes. From these two small RNA populations we have found 55 known miRNAs and identified two candidate novel miRNA families to be likely unique to the developing safflower seeds. Target genes with conserved as well as novel functions were predicted for the conserved miRNAs. We have also identified 13 miRNAs differentially expressed between the HO and HL safflower genotypes. The results may lay a foundation for unraveling the miRNA-mediated molecular processes that regulate oleic acid accumulation in the HO safflower mutant and developmental processes in safflower embryos in general.