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

Research advances of coloring mechanism regulated by MicroRNAs in plants

In plants, microRNAs (miRNAs) are a class of important small RNAs involved in their growth and development, and play a very significant role in regulating their tissue coloring. In this paper, the mechanisms on miRNA regulation of plant coloring are mainly reviewed from three aspects: macroscopic physiological and molecular foundations related to tissue coloring, miRNA biosynthesis and function, and specific analysis of miRNA regulation studies on leaf color, flower color, fruit color, and other tissue color formation in plants. Furthermore, we also systematically summarize the miRNA regulatory mechanisms identified on pigments biosynthesis and color formation in plants, and the regulatory mechanisms of these miRNAs mentioned on the existing researches can be divided into four main categories: directly targeting the related transcription factors, directly targeting the related structural genes, directly targeting the related long noncoding RNAs (LncRNAs) and miRNA-mediated production of trans-acting small interfering RNAs (ta-siRNAs). Together, these research results aim to provide a theoretical reference for the in-depth study of plant coloring mechanism and molecular breeding study of related plants in the future.

Identification and Profiling Analysis of microRNAs in Guava Fruit (Psidium guajava L.) and Their Role during Ripening

The guava (Psidium guajava L.) is a climacteric fruit with an accelerated post-harvest overripening. miRNAs are small RNA sequences that function as gene regulators in eukaryotes and are essential for their survival and development. In this study, miRNA libraries were constructed, sequenced and analyzed from the breaker and ripe stages of guava fruit cv. Siglo XXI. One hundred and seventy-four mature miRNA sequences from 28 miRNA families were identified. The taxonomic distribution of the guava miRNAs showed a high level of conservation among the dicotyledonous plants. Most of the predicted miRNA target genes were transcription factors and genes involved in the metabolism of phytohormones such as abscisic acid, auxins, and ethylene, as revealed through an ontology enrichment analysis. The miRNA families miR168, miR169, miR396, miR397, and miR482 were classified as being directly associated with maturation, whereas the miRNA families miR160, miR165, miR167, miR3930, miR395, miR398, and miR535 were classified as being indirectly associated. With this study, we intended to increase our knowledge and understanding of the regulatory process involved in the ripening process, thereby providing valuable information for future research on the ripening of guava fruit.

Metabolomic and transcriptomic analyses of the flavonoid biosynthetic pathway in blueberry (Vaccinium spp.)

As a highly economic small fruit crop, blueberry is enjoyed by most people in terms of color, taste, and rich nutrition. To better understand its coloring mechanism on the process of ripening, an integrative analysis of the metabolome and transcriptome profiles was performed in three blueberry varieties at three developmental stages. In this study, 41 flavonoid metabolites closely related to the coloring in blueberry samples were analyzed. It turned out that the most differential metabolites in the ripening processes were delphinidin-3-O-arabinoside (dpara), peonidin-3-O-glucoside (pnglu), and delphinidin-3-O-galactoside (dpgal), while the most differential metabolites among different varieties were flavonols. Furthermore, to obtain more accurate and comprehensive transcripts of blueberry during the developmental stages, PacBio and Illumina sequencing technology were combined to obtain the transcriptome of the blueberry variety Misty, for the very first time. Finally, by applying the gene coexpression network analysis, the darkviolet and bisque4 modules related to flavonoid synthesis were determined, and the key genes related to two flavonoid 3', 5'-hydroxylase (F3'5'H) genes in the darkviolet module and one bHLH transcription factor in the bisque4 module were predicted. It is believed that our findings could provide valuable information for the future study on the molecular mechanism of flavonoid metabolites and flavonoid synthesis pathways in blueberries.

Coordinating Diverse Functions of miRNA and lncRNA in Fleshy Fruit

Non-coding RNAs play vital roles in the diverse biological processes of plants, and they are becoming key topics in horticulture research. In particular, miRNAs and long non-coding RNAs (lncRNAs) are receiving increased attention in fruit crops. Recent studies in horticulture research provide both genetic and molecular evidence that miRNAs and lncRNAs regulate biological function and stress responses during fruit development. Here, we summarize multiple regulatory modules of miRNAs and lncRNAs and their biological roles in fruit sets and stress responses, which would guide the development of molecular breeding techniques on horticultural crops.

Genome-wide identification and comprehensive analysis reveal potential roles of long non-coding RNAs in fruit development of southern highbush blueberry (Vaccinium corymbosum L.)

Introduction

Blueberries have a high antioxidant content and are produced as healthy food worldwide. Long non-coding RNAs (lncRNAs) are a type of regulatory RNAs that play a variety of roles in plants. Nonetheless, information on lncRNAs and their functions during blueberry fruit development is scarce in public databases.

Methods

In the present study, we performed genome-wide identification of lncRNAs in a southern highbush blueberry using strand-specific RNA sequencing (ssRNA-Seq). Differentially expressed lncRNAs (DE-lncRNAs) and their potential target genes were analyzed at four stages of fruit development. Cis-regulatory DE-lncRNAs were predicted using co-localization analysis.

Results

These findings included a total of 25,036 lncRNAs from 17,801 loci. Blueberry lncRNAs had shorter transcript lengths, smaller open reading frame (ORF) sizes, fewer exons, and fewer isoforms than protein-coding RNAs, as well as lower expression levels and higher stage-specificity during fruit development. A total of 105 DE-lncRNAs were identified among the comparison group of PAD vs. CUP, 443 DE-lncRNAs were detected when comparing CUP with PINK fruits, and 285 DE-lncRNAs were revealed when comparing PINK and BLUE fruits. According to Kyoto Encyclopedia of Genes and Genomes annotation, target genes of DE-lncRNAs were primarily enriched in the "Autophagy-other", "DNA replication", "Endocytosis", 'photosynthesis' and 'chlorophyll metabolism' pathways, suggesting that lncRNAs may pay potential roles in fruit expansion and ripening. Moreover, several lncRNAs have been proposed as cis-regulators of the key genes involved in flavonoid biosynthesis. MSTRG.107242.6, and its putative target gene, BTB/POZ and TAZ domain-containing protein, might play critical roles in anthocyanin accumulation in blueberries.

Discussion

These findings highlight the regulatory function of lncRNAs and aid in elucidating the molecular mechanism underlying blueberry fruit growth.

MicroRNAs: emerging regulators in horticultural crops

Horticulture is one of the oldest agricultural practices with great popularity throughout the world. Horticultural crops include fruits, vegetables, ornamental plants, as well as medicinal and beverage plants. They are cultivated for food, specific nutrition, and medical use, or for aesthetic pleasure. MicroRNAs (miRNAs), which constitute a major class of endogenous small RNAs in plants, affect a multitude of developmental and physiological processes by imparting sequence specificity to gene regulation. Over the past decade, tens of thousands of miRNAs have been identified in more than 100 horticultural crops and their critical roles in regulating quality development of diverse horticultural crops have been demonstrated. Here, we review how miRNAs have emerged as important regulators and promising tools for horticultural crop improvement.

Small RNA and degradome sequencing reveal the role of blackberry miRNAs in flavonoid and anthocyanin synthesis during fruit ripening

Blackberry shrubs are economically important for their production of small, pulped fruits. This species has attracted much attention because of the unique flavor of its fruits and their rich nutritional and medicinal value. In this study, microRNAs (miRNAs) and their target genes related to flavonoids and anthocyanins in blackberry fruits during ripening were analyzed and identified by small RNA and degradome sequencing technology, and the expression levels of key miRNAs in unripe and ripe blackberry fruits were verified via the RT-qPCR. A total of 258 known miRNAs were identified. Eighty differentially expressed miRNAs (DEMs) were detected in the fruits of the ripe group compared with those of the unripe group. Differentially expressed miR828-x/miR828-z and unigene0086056 (unknown function) were coexpressed. Moreover, miR858 had the most target genes for the synthesis of flavonoids and anthocyanins. Taken together, these results provide important value for improving the quality of blackberry fruits and breeding blackberry plants that produce high-flavonoid fruits for the pharmaceutical industry.

Conservation and Diversity of miR166 Family Members From Highbush Blueberry (Vaccinium corymbosum) and Their Potential Functions in Abiotic Stress

MicroRNA166 (miR166) is highly conserved and has diverse functions across plant species. The highbush blueberry (Vaccinium corymbosum) genome is thought to harbor 10 miRNA166 loci (Vco-miR166), but the extent of their evolutionary conservation or functional diversification remains unknown. In this study, we identified six additional Vco-miR166 loci based on conserved features of the miR166 family. Phylogenetic analyses showed that mature Vco-miR166s and their precursor cluster in several clades are evolutionary conserved with diverse species. The cis-regulatory elements in the Vco-miR166 promoters indicated functions related to different phytohormones and defense responses. We also identified putative targets of vco-miR166s, which targeted the same gene families, suggesting the functional conservation and diversification of Vco-miR166 family members. Furthermore, we examined the accumulation patterns of six mature Vco-miR166s in response to abiotic stresses by stem-loop reverse RT-qPCR, which revealed their upregulation under freezing, cold, and heat stress, while they were downregulated by drought compared to control growth conditions. However, Vco-miR166 members showed different expression patterns when exposed to salt stress. These results showed that conserved Vco-miR166 family members display functional diversification but also coordinately influence plant responses to abiotic stress.

Identification of ARF family in blueberry and its potential involvement of fruit development and pH stress response

Background

Auxin responsive factor (ARF) family is one of core components in auxin signalling pathway, which governs diverse developmental processes and stress responses. Blueberry is an economically important berry-bearing crop and prefers to acidic soil. However, the understandings of ARF family has not yet been reported in blueberry.

Results

In the present study, 60 ARF genes (VcARF) were identified in blueberry, and they showed diverse gene structures and motif compositions among the groups and similar within each group in the phylogenetic tree. Noticeably, 9 digenic, 5 trigenic and 6 tetragenic VcARF pairs exhibited more than 95% identity to each other. Computational analysis indicated that 23 VcARFs harbored the miRNA responsive element (MRE) of miR160 or miR167 like other plant ARF genes. Interestingly, the MRE of miR156d/h-3p was observed in the 5’UTR of 3 VcARFs, suggesting a potentially novel post-transcriptional control. Furthermore, the transcript accumulations of VcARFs were investigated during fruit development, and three categories of transcript profiles were observed, implying different functional roles. Meanwhile, the expressions of VcARFs to different pH conditions (pH4.5 and pH6.5) were surveyed in pH-sensitive and tolerant blueberry species, and a number of VcARFs showed different transcript accumulations. More importantly, distinct transcriptional response to pH stress (pH6.5) were observed for several VcARFs (such as VcARF6s and VcARF19-3/19–4) between pH-sensitive and tolerant species, suggesting their potential roles in adaption to pH stress.

Conclusions

Sixty VcARF genes were identified and characterized, and their transcript profiles were surveyed during fruit development and in response to pH stress. These findings will contribute to future research for eliciting the functional roles of VcARFs and regulatory mechanisms, especially fruit development and adaption to pH stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08556-y.

Transcriptomic and Physiological Analysis Reveals the Responses to Auxin and Abscisic Acid Accumulation During Vaccinium corymbosum Flower Bud and Fruit Development

Blueberry (Vaccinium corymbosum) is reputed as a rich source of health-promoting phytonutrients, which contributes to its burgeoning consumer demand and production. However, blueberries are much smaller and have lower yields than most domesticated berries, and the inherent regulatory mechanisms remain elusive. In this study, the cytological and physiological changes, as well as comparative transcriptomic analysis throughout flower and fruit development in the southern highbush blueberry cultivar 'O'Neal' were performed. 'O'Neal' hypanthium and fruit exhibited a distinctive cell proliferation pattern, and auxin accumulation was unusual throughout development, while abscisic acid (ABA) levels rapidly increased in association with anthocyanin accumulation, total phenolic reduction and fruit maturation. Transcriptomic data showed that many differentially expressed genes (DEGs) were specifically expressed at each flower bud and fruit developmental stage. Further weighted gene co-expression network analysis (WGCNA) revealed numerous DEGs that correlated with the cell numbers of outer mesocarp and columella, showed two distinctive expression patterns. Most of the DEGs involved in auxin biosynthesis, transportation and signal transduction were upregulated, and this upregulation was accompanied by cell expansion, and flower bud and fruit development. However, individual members of VcSAUR50 and VcIAA9 families might be insensitive to auxin, suggesting that these genes play a distinctive role in the growth and development of blueberry fruits. These results will support future research to better understand the flower and fruit development of southern highbush blueberry.

Comprehensive Identification and Profiling of miRNAs Involved in Terpenoid Synthesis of Gleditsia sinensis Lam

Gleditsia sinensis Lam. is a tree with worldwide distribution and important economic and medicinal values; its pods contain terpenoids including gleditsioside, thiamine, and brassinosteroids. However, thus far, there are few studies on the terpenoid regulation of G. sinensis at the molecular level. microRNA (miRNA) is a class of small RNAs with conserved and crucial roles in the regulation of diverse biological processes during plant growth and development. To identify the miRNAs of G. sinensis and evaluate their involvement in terpenoid synthesis, this investigation quantified the content changes in saponins in pods at three developmental stages: May (pod-setting stage), July (elongation stage), and September (browning stage), and then we performed genome-wide miRNA profiles during the three development stages of the G. sinensis pods. A total of 351 conserved miRNAs belonging to 216 families were identified, among which 36 conserved miRNAs exist specifically in legumes. Through target analysis, 708 unigenes were predicted to be candidate targets of 37 differentially expressed miRNAs. The targets of miR838-3p and miR2093-5p were involved in the derived branches of monoterpenes and gleditsioside, in brassinosteroid biosynthesis (BRB), and in indole alkaloid biosynthesis (IAB). Intriguingly, the targets of miR829-3p.1 were predicted to take part in thiamine biosynthesis, and the targets of miR4414b and miR5037a were involved in the main process of cytokinin synthesis. The corresponding targets participated in BRB, IAB, and terpenoid backbone biosynthesis, which were enriched significantly, suggesting that miR2093-5p, miR4414b, miR5037a, miR829-3p.1, and miR838-3p play indispensable roles in the regulation of triterpenoid saponin and monoterpenoid biosynthesis. To date, this is the first report of miRNA identification in G. sinensis and miRNA expression profiles at different developmental stages of G. sinensis pods, which provides a basis for further uncovering the molecular regulation of terpenoid synthesis in G. sinensis and new insights into the role of miRNAs in legumes.

Dynamic regulation of small RNAs in anthocyanin accumulation during blueberry fruit maturation

Blueberry is rich in anthocyanins which accumulate during fruit maturation. Previous studies mostly focus on their translational/transcriptional regulation, but usually underestimate their post-transcriptional regulation, e.g. small RNAs. This study aimed to identify sRNAs and their potential pathways associated with anthocyanin biosynthesis. During three typical phases of fruit maturation (green, pink, and blue), we investigated dynamic changes of sRNA by deep sequencing sRNA and examined the interaction of sRNAs with their target genes by degradome and RLM-PCR. During maturation, up-regulation of VcmiRNA156 and VcmiR393 resulted in down-regulation of VcSPLs and VcTIR1/AFBs, respectively. An important gene of anthocyanin biosynthesis, VcDFR, was substantially down-regulated at both the mRNA and protein levels, and potentially responded to regulation of VcSPLs and VcTIR1/AFBs. Additionally, indole acetic acid (IAA) and abscisic acid (ABA) were involved in the regulation of anthocyanin biosynthesis by interacting with VcmiR393-TIR1/AFBs and VcmiRNA319-VcMYBs respectively. This information provides another insight into blueberry anthocyanin biosynthesis.

Comprehensive Analysis of the SBP Family in Blueberry and Their Regulatory Mechanism Controlling Chlorophyll Accumulation

SQUAMOSA Promoter Binding Protein (SBP) family genes act as central players to regulate plant growth and development with functional redundancy and specificity. Addressing the diversity of the SBP family in crops is of great significance to precisely utilize them to improve agronomic traits. Blueberry is an important economic berry crop. However, the SBP family has not been described in blueberry. In the present study, twenty VcSBP genes were identified through data mining against blueberry transcriptome databases. These VcSBPs could be clustered into eight groups, and the gene structures and motif compositions are divergent among the groups and similar within each group. The VcSBPs were differentially expressed in various tissues. Intriguingly, 10 VcSBPs were highly expressed at green fruit stages and dramatically decreased at the onset of fruit ripening, implying that they are important regulators during early fruit development. Computational analysis showed that 10 VcSBPs were targeted by miR156, and four of them were further verified by degradome sequencing. Moreover, their functional diversity was studied in Arabidopsis. Noticeably, three VcSBPs significantly increased chlorophyll accumulation, and qRT-PCR analysis indicated that VcSBP13a in Arabidopsis enhanced the expression of chlorophyll biosynthetic genes such as AtDVR, AtPORA, AtPORB, AtPORC, and AtCAO. Finally, the targets of VcSBPs were computationally identified in blueberry, and the Y1H assay showed that VcSBP13a could physically bind to the promoter region of the chlorophyll-associated gene VcLHCB1. Our findings provided an overall framework for individually understanding the characteristics and functions of the SBP family in blueberry.

Identification of miRNAs and Their Targets Involved in Flower and Fruit Development across Domesticated and Wild Capsicum Species

MicroRNAs (miRNAs) are regulators of the post-transcription stage of gene activity documented to play central roles in flower and fruit development in model plant species. However, little is known about their roles and differences in domesticated and wild Capsicum species. In this study, we used high-throughput sequencing to analyze the miRNA content at three developmental stages (flower, small fruit, and middle fruit) from two cultivated (C. baccatum and C. annuum) and two wild (C. chacoense and C. eximium) pepper species. This analysis revealed 22 known and 27 novel miRNAs differentially expressed across species and tissues. A number of stage- and species-specific miRNAs were identified, and Gene Ontology terms were assigned to 138 genes targeted by the miRNAs. Most Gene Ontology terms were for the categories "genetic information processing", "signaling and cellular processes", "amino acid metabolism", and "carbohydrate metabolism". Enriched KEGG analysis revealed the pathways amino acids, sugar and nucleotide metabolism, starch and sucrose metabolism, and fructose-mannose metabolism among the principal ones regulated by miRNAs during pepper fruit ripening. We predicted miRNA-target gene interactions regulating flowering time and fruit development, including miR156/157 with SPL genes, miR159 with GaMYB proteins, miR160 with ARF genes, miR172 with AP2-like transcription factors, and miR408 with CLAVATA1 gene across the different Capsicum species. In addition, novel miRNAs play an important role in regulating interactions potentially controlling plant pathogen defense and fruit quality via fructokinase, alpha-L-arabinofuranosidase, and aromatic and neutral amino acid transporter. Overall, the small RNA-sequencing results from this study represent valuable information that provides a solid foundation for uncovering the miRNA-mediated mechanisms of flower and fruit development between domesticated and wild Capsicum species.

A blueberry MIR156a-SPL12 module coordinates the accumulation of chlorophylls and anthocyanins during fruit ripening

Color change is an important event during fruit maturation in blueberry, usually depending on chlorophyll degradation and anthocyanin accumulation. MicroRNA156 (miR156)-SPL modules are an important group of regulatory hubs involved in the regulation of anthocyanin biosynthesis. However, little is known regarding their roles in blueberry or in chlorophyll metabolism during color change. In this study, a MIR156 gene (VcMIR156a) was experimentally identified in blueberry (Vaccinium corymbosum). Overexpression of VcMIR156a in tomato (Solanum lycopersicum) enhanced anthocyanin biosynthesis and chlorophyll degradation in the stem by altering pigment-associated gene expression. Further investigation indicated that the VcSPL12 transcript could be targeted by miR156, and showed the reverse accumulation patterns during blueberry fruit development and maturation. Noticeably, VcSPL12 was highly expressed at green fruit stages, while VcMIR156a transcripts mainly accumulated at the white fruit stage when expression of VcSPL12 was dramatically decreased, implying that VcMIR156a-VcSPL12 is a key regulatory hub during fruit coloration. Moreover, VcSPL12 decreased the expression of several anthocyanin biosynthetic and regulatory genes, and a yeast two-hybrid assay indicated that VcSPL12 interacted with VcMYBPA1. Intriguingly, expression of VcSPL12 significantly enhanced chlorophyll accumulation and altered the expression of several chlorophyll-associated genes. Additionally, the chloroplast ultrastructure was altered by the expression of VcMIR156a and VcSPL12. These findings provide a novel insight into the functional roles of miR156-SPLs in plants, especially in blueberry fruit coloration.

Cross-Kingdom Regulation by Plant microRNAs Provides Novel Insight into Gene Regulation

microRNAs (miRNAs) are well known as major players in mammalian and plant genetic systems that act by regulating gene expression at the post-transcriptional level. These tiny molecules can regulate target genes (mRNAs) through either cleavage or translational inhibition. Recently, the discovery of plant-derived miRNAs showing cross-kingdom abilities to regulate mammalian gene expression has prompted exciting discussions among researchers. After being acquired orally through the diet, plant miRNAs can survive in the digestive tract, enter the circulatory system, and regulate endogenous mRNAs. Here, we review current knowledge regarding the cross-kingdom mechanisms of plant miRNAs, related controversies, and potential applications of these miRNAs in dietary therapy, which will provide new insights for plant miRNA investigations related to health issues in humans.

Single-base cytosine methylation analysis in fruits of three Capsicum species

Single-base cytosine methylation analysis across fruits of Capsicum annuum, C. chinense and C. frutescens showed global average methylation ranging from 82.8-89.1%, 77.6-83.9%, and 22.4-25% at CG, CHG and CHH contexts, respectively. High gene-body methylation at CG and CHG was observed across Capsicum species. The C. annuum showed the highest proportion (>80%) of mCs at different genomic regions compared to C. chinense and C. frutescens. Cytosine methylation for transposable-elements were lower in C. frutescens compared to C. annuum and C. chinense. A total of 510,165 CG, 583112 CHG and 277,897 CHH DMRs were identified across three Capsicum species. The differentially methylated regions (DMRs) distribution analysis revealed C. frutescens as more hypo-methylated compared to C. annuum and C. chinense, and also the presence of more intergenic DMRs in Capsicum genome. At CG and CHG context, gene expression and promoter methylation showed inverse correlations. Furthermore, the observed correlation between methylation and expression of genes suggested the potential role of methylation in Capsicum fruit development/ripening.

Roles of transcription factor SQUAMOSA promoter binding protein-like gene family in papaya (Carica papaya) development and ripening

SQUAMOSA promoter binding protein-like (SPL) family plays vital regulatory roles in plant growth and development. The SPL family in climacteric fruit Carica papaya has not been reported. This study identified 14 papaya SPLs (CpSPL) from papaya genome and analyzed their sequence features, phylogeny, intron/exon structure, conserved motif, miR156-mediated posttranscriptional regulation, and expression patterns. 14 CpSPLs were clustered into 8 groups, and two distinct expression patterns were revealed for miR156-targeted and nontargeted CpSPLs in different tissues and fruit development stages. The expression changes of CpSPLs in ethephon and 1-MCP treated fruit during ripening suggested that the CpSPLs guided by CpmiR156 play crucial roles in ethylene signaling pathway. This study sheds light on the new function of SPL family in fruit development and ripening, providing insights on understanding evolutionary divergence of the members of SPL family among plant species.

Genome-Wide Identification of lncRNAs and Analysis of ceRNA Networks During Tomato Resistance to Phytophthora infestans

Our previous studies have revealed the function of long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) in tomato in response to Phytophthora infestans infection. However, the interaction relationships between lncRNAs and miRNAs during tomato resistance to P. infestans infection are unknown. In this study, 9,011 lncRNAs were identified from tomato plants, including 115 upregulated and 81 downregulated lncRNAs. Among these, 148 were found to be differentially expressed and might affect the expression of 771 genes, which are composed of 887 matched lncRNA-mRNA pairs. In total, 88 lncRNAs were identified as endogenous RNAs (ceRNAs) and predicted to decoy 46 miRNAs. Degradome sequencing revealed that 11 miRNAs that were decoyed by 20 lncRNAs could target 30 genes. These lncRNAs, miRNAs, and target genes were predicted to form 10 regulatory modules. Among them, lncRNA42705/lncRNA08711, lncRNA39896, and lncRNA11265/lncRNA15816 might modulate MYB, HD-Zip, and NAC transcription factors by decoying miR159, miR166b, and miR164a-5p, respectively. Upon P. infestans infection, the expression levels of lncRNA42705 and lncRNA08711 displayed a negative correlation with the expression level of miR159 and a positive correlation with the expression levels of MYB genes. Tomato plants in which lncRNA42705 and lncRNA08711 were silenced displayed increased levels of miR159 and decreased levels of MYB, respectively. The result demonstrated that lncRNAs might function as ceRNAs to decoy miRNAs and affect their target genes in tomato plants, increasing resistance to disease.

Identification and functional characterization of the Aux/IAA gene VcIAA27 in blueberry

Aux/IAA genes are an important class of players in diverse developmental processes in plants, which generally exert their functions through the auxin signaling pathway. Blueberry is an economically and nutritionally important berry-bearing crop. However, Aux/IAA genes remain unknown in blueberry. In the present study, an Aux/IAA gene (VcIAA27) was identified and characterized in blueberry, and it is most closely related to IAA27 in other plant species. Expression analysis indicated that VcIAA27 transcripts accumulate highly in shoot, flower and fruit. Interestingly, VcIAA27 was highly expressed at early fruit developmental stages, and dramatically decreased from the onset of fruit ripening, implying that VcIAA27 possibly plays important roles during fruit enlargement. Meanwhile, the analysis of promoter activity in Arabidopsis showed that strong GUS signal was detected in the trichome and hydathodes of leaves, receptacle of silique, and lateral roots of seedling. Overexpression of VcIAA27 in Arabidopsis leads to auxin-related defects such as downward-curled leaves, short or sterile siliques, shorter stature, and more shoot branches. Moreover, qPCR analysis indicated that VcIAA27 is able to alter the expression patterns of the auxin-related genes BRU6, SAG13, SAUR26 in Arabidopsis, suggesting that VcIAA27 might be negatively involved in the auxin signaling pathway. The findings will greatly contribute to future investigation of Aux/IAA-mediated mechanisms that control blueberry development, especially fruit development and ripening.

Network analysis of noncoding RNAs in pepper provides insights into fruit ripening control

Pepper is an important vegetable worldwide and is a model plant for nonclimacteric fleshy fruit ripening. Drastic visual changes and internal biochemical alterations are involved in fruit coloration, flavor, texture, aroma, and palatability to animals during the pepper fruit ripening process. To explore the regulation of bell pepper fruit ripening by noncoding RNAs (ncRNAs), we examined their expression profiles; 43 microRNAs (miRNAs), 125 circular RNAs (circRNAs), 366 long noncoding RNAs (lncRNAs), and 3266 messenger RNAs (mRNAs) were differentially expressed (DE) in mature green and red ripe fruit. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the targets of the DE ncRNAs and DE mRNAs included several kinds of transcription factors (TFs) (ERF, bHLH, WRKY, MYB, NAC, bZIP, and ARF), enzymes involved in cell wall metabolism (beta-galactosidase, beta-glucosidase, beta-amylase, chitinase, pectate lyase (PL), pectinesterase (PE) and polygalacturonase (PG)), enzymes involved in fruit color accumulation (bifunctional 15-cis-phytoene synthase, 9-cis-epoxycarotenoid dioxygenase, beta-carotene hydroxylase and carotene epsilon-monooxygenase), enzymes associated with fruit flavor and aroma (glutamate-1-semialdehyde 2,1-aminomutase, anthocyanin 5-aromatic acyltransferase, and eugenol synthase 1) and enzymes involved in the production of ethylene (ET) (ACO1/ACO4) as well as other plant hormones such as abscisic acid (ABA), auxin (IAA), and gibberellic acid (GA). Based on accumulation profiles, a network of ncRNAs and mRNAs associated with bell pepper fruit ripening was developed that provides a foundation for further developing a more refined understanding of the molecular biology of fruit ripening.

VvmiR160s/VvARFs interaction and their spatio-temporal expression/cleavage products during GA-induced grape parthenocarpy

BACKGROUND

Grape (Vitis vinifera) is highly sensitive to gibberellin (GA), which effectively induce grape parthenocarpy. Studies showed that miR160s and their target AUXIN RESPONSIVE FACTOR (ARF) responding hormones are indispensable for various aspects of plant growth and development, but their functions in GA-induced grape parthenocarpy remain elusive.

RESULTS

In this study, the morphological changes during flower development in response to GA treatments were examined in the 'Rosario Bianco' cultivar. The precise sequences of VvmiR160a/b/c/d/e and their VvARF10/16/17 target genes were cloned, sequenced and characterized. The phylogenetic relationship and intron-exon structure of VvARFs and other ARF family members derived from different species were investigated. All VvmiR160s (except VvmiR160b) and VvARF10/16/17 had the common cis-elements responsive to GA, which support their function in GA-mediated grape parthenocarpy. The cleavage role of VvmiR160s-mediated VvARF10/16/17 was verified in grape flowers. Moreover, spatio-temporal expression analysis demonstrated that among VvmiR160 family, VvmiR160a/b/c highly expressed at late stage of flower/berry development, while VvARF10/16/17showed a reverse expression trend. Interestingly, GA exhibited a long-term effect through inducing the expression of VvmiR160a/b/c/e to increase their cleavage product accumulations from 5 to 9 days after treatment, but GA enhanced the expressions of VvARF10/16/17 only at short term. Pearson correlation analysis based on expression data revealed a negative correlation between VvmiR160a/b/c and VvARF10/16/17 in flowers not berries during GA-induced grape parthenocarpy.

CONCLUSIONS

This work demonstrated that the negative regulation of VvARF10/16/17 expression by VvmiR160a/b/c as key regulatory factors is critical for GA-mediated grape parthenocarpy, and provide significant implications for molecular breeding of high-quality seedless berry.

Identification of Blueberry miRNAs and Their Targets Based on High-Throughput Sequencing and Degradome Analyses

miRNAs are important regulators of plant gene expression. To better characterize their functions, we applied high-throughput sequencing and degradome analyses to investigate three blueberry (Vaccinium ashei) tissues. A total of 127 known and 101 novel miRNAs were identified. Moreover, 141 targets for 42 known and 19 novel miRNAs were experimentally validated by degradome sequencing. A functional analysis of these miRNA targets revealed they were associated with diverse biological activities and several pathways, e.g., anthocyanin biosynthesis and cytokinin signal transduction. The data presented herein expand our understanding of the regulation of blueberry miRNAs during floral and fruit development stages. They may also provide new insights into the roles of miRNAs during anthocyanin biosynthesis in blueberry fruits.