Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening

Noncoding RNAs in regulation of plant secondary metabolism

Plant secondary metabolites (PSMs) are a large class of structurally diverse molecules, mainly consisting of terpenoids, phenolic compounds, and nitrogen-containing compounds, which play active roles in plant development and stress responses. The biosynthetic processes of PSMs are governed by a sophisticated regulatory network at multiple levels. Noncoding RNAs (ncRNAs) such as microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs) may serve as post-transcriptional regulators for plant secondary metabolism through acting on genes encoding either transcription factors or participating enzymes in relevant metabolic pathways. High-throughput sequencing technologies have facilitated the large-scale identifications of ncRNAs potentially involved in plant secondary metabolism in model plant species as well as certain species with enriched production of specific types of PSMs. Moreover, a series of miRNA-target modules have been functionally characterized to be responsible for regulating PSM biosynthesis and accumulation in plants under abiotic or biotic stresses. In this review, we will provide an overview of current findings on the ncRNA-mediated regulation of plant secondary metabolism with special attention to its participation in plant stress responses, and discuss possible issues to be addressed in future fundamental research and breeding practice.

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Integrated Analysis of lncRNA–mRNA Regulatory Networks Related to Lipid Metabolism in High-Oleic-Acid Rapeseed

A high oleic acid content is considered an essential characteristic in the breeding of high-quality rapeseed in China. Long-chain non-coding RNA (lncRNA) molecules play an important role in the plant’s growth and its response to stress. To better understand the role of lncRNAs in regulating plant reproductive development, we analyzed whole-transcriptome and physiological data to characterize the dynamic changes in lncRNA expression during the four representative times of seed development of high- and low-oleic-acid rapeseed in three regions. We identified 21 and 14 lncRNA and mRNA modules, respectively. These modules were divided into three types related to region, development stages, and material. Next, we analyzed the key modules related to the oil content and the oleic acid, linoleic acid, and linolenic acid contents with physiological data and constructed the key functional network analysis on this basis. Genes related to lipid metabolism, such as 3-ketoacyl-CoA synthase 16 (KCS16) and acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1), were present in the co-expression network, suggesting that the effect of these genes on lipid metabolism might be embodied by the expression of these lncRNAs. Our results provide a fresh insight into region-, development-stage-, and material-biased changes in lncRNA expression in the seeds of Brassica napus. Some of these lncRNAs may participate in the regulatory network of lipid accumulation and metabolism, together with regulated genes. These results may help elucidate the regulatory system of lncRNAs in the lipid metabolism of high-oleic-acid rapeseed seeds.

Promoter Variation of the Key Apple Fruit Texture Related Gene MdPG1 and the Upstream Regulation Analysis

MdPG1 encoding polygalacturonase in apple (Malus × domestica) is a key gene associated with fruit firmness and texture variations among apple cultivars. However, the causative variants of MdPG1 are still not known. In this study, we identified a SNP^A/C variant within an ERF-binding element located in the promoter region of MdPG1. The promoter containing the ERF-binding element with SNP^A, rather than the SNP^C, could be strongly bound and activated by MdCBF2, a member of the AP2/ERF transcription factor family, as determined by yeast-one-hybrid and dual-luciferase reporter assays. We also demonstrated that the presence of a novel long non-coding RNA, lncRNA_PG1, in the promoter of MdPG1 was a causative variant. lncRNA_PG1 was specifically expressed in fruit tissues postharvest. lncRNA_PG1 could reduce promoter activity when it was fused to the promoter of MdPG1 and a tobacco gene encoding Mg-chelatase H subunit (NtCHLH) in transgenic tobacco cells but could not reduce promoter activity when it was supplied in a separate gene construct, indicating a cis-regulatory effect. Our results provide new insights into genetic regulation of MdPG1 allele expression and are also useful for the development of elite apple cultivars.

Targeted multi-platform metabolome analysis and enzyme activity analysis of kiwifruit during postharvest ripening

Kiwifruit is a climacteric fruit, in which the accumulation of flavor substances mainly occurs at the postharvest ripening stage. However, the dynamic changes in metabolite composition remain poorly understood. Here, targeted multi-platform metabolome analysis based on GC-MS and UPLC-MS/MS and enzyme activity analysis were performed at different postharvest ripening stages of kiwifruit. A total of 12 soluble sugars and 31 organic acids were identified. The main soluble sugars are sucrose, glucose and fructose, which exhibited similar variation tendencies along with the extension of ripening. The main organic acids are citric acid, quinic acid and malic acid, which showed different variation patterns. A total of 48 energy metabolites were identified, which were classified into two groups based on the content variation. The content of substances related to the respiratory metabolic pathway decreased gradually along with postharvest ripening, and there was obvious accumulation of downstream products such as amino acids at the late ripening stage. A total of 35 endogenous hormones were identified, among which seven cytokinins were highly accumulated at the later stage of softening. We further investigated the dynamic changes in the activities of 28 ripening-related enzymes. As a result, the activities of 13 enzymes were highly correlated with changes in starch, total pectin, and soluble sugars, and those of seven enzymes were closely associated with the change in firmness. In conclusion, this study comprehensively describes the dynamic changes in soluble sugars, organic acids, hormones, energy substances, and ripening-related enzyme activities during kiwifruit postharvest ripening, and provides a theoretical basis for the postharvest quality improvement of kiwifruit.

Identification and Characterization of Long Non-Coding RNAs: Implicating Insights into Their Regulatory Role in Kiwifruit Ripening and Softening during Low-Temperature Storage

Long non-coding RNAs (lncRNAs) are crucial players regulating many biological processes in plants. However, limited knowledge is available regarding their roles in kiwifruit ripening and softening. In this study, using lncRNA-seq technology, 591 differentially expressed (DE) lncRNAs (DELs) and 3107 DE genes (DEGs) were identified from kiwifruit stored at 4 °C for 1, 2, and 3 weeks in comparison with non-treated control fruits. Of note, 645 DEGs were predicted to be targets of DELs (DEGTLs), including some DE protein-coding genes (such as β-amylase and pectinesterase). DEGTL-based GO enrichment analysis revealed that these genes were significantly enriched in cell wall modification and pectinesterase activity in 1 W vs. CK and 3 W vs. CK, which might be closely related to the fruit softening during low-temperature storage. Moreover, KEGG enrichment analysis revealed that DEGTLs were significantly associated with starch and sucrose metabolism. Our study revealed that lncRNAs play critical regulatory roles in kiwifruit ripening and softening under low-temperature storage, mainly by mediating the expression of starch and sucrose metabolism and cell wall modification related genes.

Advances in Plant Epigenome Editing Research and Its Application in Plants

Plant epistatic regulation is the DNA methylation, non-coding RNA regulation, and histone modification of gene sequences without altering the genome sequence, thus regulating gene expression patterns and the growth process of plants to produce heritable changes. Epistatic regulation in plants can regulate plant responses to different environmental stresses, regulate fruit growth and development, etc. Genome editing can effectively improve plant genetic efficiency by targeting the design and efficient editing of genome-specific loci with specific nucleases, such as zinc finger nucleases (ZFNs), transcription activator-like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9). As research progresses, the CRISPR/Cas9 system has been widely used in crop breeding, gene expression, and epistatic modification due to its high editing efficiency and rapid translation of results. In this review, we summarize the recent progress of CRISPR/Cas9 in epigenome editing and look forward to the future development direction of this system in plant epigenetic modification to provide a reference for the application of CRISPR/Cas9 in genome editing.

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.

Genome-wide identification and characterization of long noncoding RNAs during peach (Prunus persica) fruit development and ripening

LncRNAs represent a class of RNA transcripts of more than 200 nucleotides (nt) in length without discernible protein-coding potential. The expression levels of lncRNAs are significantly affected by stress or developmental cues. Recent studies have shown that lncRNAs participate in fruit development and ripening processes in tomato and strawberry; however, in other fleshy fruits, the association between lncRNAs and fruit ripening remains largely elusive. Here, we constructed 9 ssRNA-Seq libraries from three different peach (Prunus persica) fruit developmental stages comprising the first and second exponential stages and the fruit-ripening stage. In total, 1500 confident lncRNAs from 887 loci were obtained according to the bioinformatics analysis. The lncRNAs identified in peach fruits showed distinct characteristics compared with protein-coding mRNAs, including lower expression levels, lower complexity of alternative splicing, shorter isoforms and smaller numbers of exons. Expression analysis identified 575 differentially expressed lncRNAs (DELs) classified into 6 clusters, among which members of Clusters 1, 2, 4 and 5 were putatively associated with fruit development and ripening processes. Quantitative real-time PCR revealed that the DELs indeed had stage-specific expression patterns in peach fruits. GO and KEGG enrichment analysis revealed that DELs might be associated with fruit-ripening-related physiological and metabolic changes, such as flavonoid biosynthesis, fruit texture softening, chlorophyll breakdown and aroma compound accumulation. Finally, the similarity analysis of lncRNAs within different plant species indicated the low sequence conservation of lncRNAs. Our study reports a large number of fruit-expressed lncRNAs and identifies fruit development phase-specific expressed lncRNA members, which highlights their potential functions in fruit development and ripening processes and lays the foundations for future functional research.

Integrative analysis of transcriptomic and metabolomic profiles reveals new insights into the molecular foundation of fruit quality formation in Citrullus lanatus (Thunb.) Matsum. & Nakai

In this study, an integrated transcriptome and metabolome analysis was used to explore the molecular foundation of fruit quality in two parent lines of Citrullus lanatus with distinct flesh characteristics, including ‘14-1’ (sweet, red, and soft) and ‘W600’ (bitter, light yellow, and firm), as well as the corresponding F1 population (bitter, light yellow, and firm). Numerous differentially expressed genes (DEGs) were identified in the fruit samples: 3,766 DEGs for ‘14-1’ vs. ‘W600’, 2,767 for ‘14-1’ vs. F1, and 1,178 for F1 vs. ‘W600’ at the transition stage; and 4,221 for ‘14-1’ vs. ‘W600’, 2,447 for ‘14-1’ vs. F1, and 446 for F1 vs. ‘W600’ at the maturity stage. Weighted gene co-expression network analysis (WGCNA) revealed that a gene module including 1,111 DEGs was closely associated with flesh taste and color, and another gene module including 1,575 DEGs contributed significantly to flesh texture. The metabolomic results showed that there were 447 differential metabolites (DMs) for ‘14-1’ vs. ‘W600’ fruits, 394 for ‘14-1’ vs. F1, and 298 for F1 vs. ‘W600’ at the maturity stage. Combining WGNCA and metabolomic results, several DEGs and DMs were further identified as hub players in fruit quality formation: six DEGs with four DMs for flesh sweetness; six DEGs with 13 DMs for bitterness; nine DEGs with 10 DMs for flesh color; and nine DEGs with four DMs for flesh texture. Altogether, these observations not only expand our knowledge of the molecular basis of fruit quality in watermelon, but also provide potential targets for future watermelon improvement.

Transcriptome profiling of genes associated with fruit firmness in the melon variety 'Baogua' (Cucumis melo ssp. agrestis Jeffrey)

Fruit firmness is an important trait of melons due to its effect on fresh fruit consumption, storage, and transport. However, information on the expression of genes influencing the fruit firmness of 'Baogua' (BG) melon (Cucumis melo ssp. agrestis Jeffrey) remains rare. This study aimed to identify the key genes associated with the firmness of BG fruit sampled at 14 and 28 days after pollination (dap) via transcriptome sequencing. A total of 1113 up-regulated and 2224 down-regulated differentially expressed genes (DEGs) were identified. The main Gene Ontology terms assigned to the DEGs were phosphotransferase activity, alcohol group as acceptor, protein phosphorylation, and protein kinase activity. The enriched KEGG pathways involving the DEGs were starch and sucrose metabolism, diterpenoid biosynthesis, plant hormone signal transduction, and MAPK signaling pathway-plant. In addition, qRT-PCR verified that four GAL genes, namely, CmGAL1-4, were differentially expressed at 0, 7, 14, 21, and 28 dap. Our data revealed that CmGAL1 expression was highest at 21 dap. However, the expression levels of CmGAL2-4 were highest at 14 dap. The sequence of CmGAL1 was similar to the sequences of homologs from melon and cucumber. Subcellular localization analysis revealed CmGAL1 was located in the cell membrane and cytoplasm. Our findings implied that fruit development at 14 dap, which is a key time-point, varies considerably from fruit development at 28 dap. Our present study provides new information on the genes associated with BG fruit firmness and help improve the storage and transport of BG fruit prior to processing.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01131-5.

Identification of Long Non-Coding RNAs Associated with Tomato Fruit Expansion and Ripening by Strand-Specific Paired-End RNA Sequencing

As emerging essential regulators in plant development, long non-coding RNAs (lncRNAs) have been extensively investigated in multiple horticultural crops, as well as in different tissues of plants. Tomato fruits are an indispensable part of people’s diet and are consumed as fruits and vegetables. Meanwhile, tomato is widely used as a model to study the ripening mechanism in fleshy fruit. Although increasing evidence shows that lncRNAs are involved in lots of biological processes in tomato plants, the comprehensive identification of lncRNAs in tomato fruit during its expansion and ripening and their functions are partially known. Here, we performed strand-specific paired-end RNA sequencing (ssRNA-seq) of tomato Heinz1706 fruits at five different developmental stages, as well as flowers and leaves. We identified 17,674 putative lncRNAs by referencing the recently released SL4.0 and annotation ITAG4.0 in tomato plants. Many lncRNAs show different expression patterns in fleshy fruit at different developmental stages compared with leaves or flowers. Our results indicate that lncRNAs play an important role in the regulation of tomato fruit expansion and ripening, providing informative lncRNA candidates for further studies in tomato fruits. In addition, we also summarize the recent advanced progress in lncRNAs mediated regulation on horticultural fruits. Hence, our study updates the understanding of lncRNAs in horticultural plants and provides resources for future studies relating to the expansion and ripening of tomato fruits.

Nanopore-Based Comparative Transcriptome Analysis Reveals the Potential Mechanism of High-Temperature Tolerance in Cotton (Gossypium hirsutum L.)

Extreme high temperatures are threatening cotton production around the world due to the intensification of global warming. To cope with high-temperature stress, heat-tolerant cotton cultivars have been bred, but the heat-tolerant mechanism remains unclear. This study selected heat-tolerant (‘Xinluzao36′) and heat-sensitive (‘Che61-72′) cultivars of cotton treated with high-temperature stress as plant materials and performed comparative nanopore sequencing transcriptome analysis to reveal the potential heat-tolerant mechanism of cotton. Results showed that 120,605 nonredundant sequences were generated from the raw reads, and 78,601 genes were annotated. Differentially expressed gene (DEG) analysis showed that a total of 19,600 DEGs were screened; the DEGs involved in the ribosome, heat shock proteins, auxin and ethylene signaling transduction, and photosynthesis pathways may be attributed to the heat tolerance of the heat-tolerant cotton cultivar. This study also predicted a total of 5118 long non-coding RNAs (lncRNAs)and 24,462 corresponding target genes. Analysis of the target genes revealed that the expression of some ribosomal, heat shock, auxin and ethylene signaling transduction-related and photosynthetic proteins may be regulated by lncRNAs and further participate in the heat tolerance of cotton. This study deepens our understandings of the heat tolerance of cotton.

Advances in Non-Coding RNA Sequencing

Non-coding RNAs (ncRNAs) comprise a set of abundant and functionally diverse RNA molecules. Since the discovery of the first ncRNA in the 1960s, ncRNAs have been shown to be involved in nearly all steps of the central dogma of molecular biology. In recent years, the pace of discovery of novel ncRNAs and their cellular roles has been greatly accelerated by high-throughput sequencing. Advances in sequencing technology, library preparation protocols as well as computational biology helped to greatly expand our knowledge of which ncRNAs exist throughout the kingdoms of life. Moreover, RNA sequencing revealed crucial roles of many ncRNAs in human health and disease. In this review, we discuss the most recent methodological advancements in the rapidly evolving field of high-throughput sequencing and how it has greatly expanded our understanding of ncRNA biology across a large number of different organisms.