Comparative Transcriptional Analysis of Loquat Fruit Identifies Major Signal Networks Involved in Fruit Development and Ripening Process

Transcriptomic Analysis Reveals Potential Gene Regulatory Networks Under Cold Stress of Loquat (Eriobotrya japonica Lindl.)

Loquat (Eriobotrya japonica Lindl. ) is one of the most economically important evergreen fruit crops in China, while it often suffered the injury of cold stress in winter and earlier spring, and the annual yield loss of loquat fruits caused by cold or freezing stress was immeasurable. However, knowledge about the physiological response and molecular mechanism under cold stress is still limited. To investigate the potential regulation mechanism pre- and post-cold stress in loquat and the changes in physiological indicators, a comparative transcriptome analysis was performed against a cold-resistant cv. "Huoju" and a cold-sensitive cv. "Ninghaibai". The results of physiological indicators related to cold resistance indicated that rachis was most sensitive to cold stress and was considered as the representative organ to directly evaluate cold resistance of loquat based on subordinate function analysis. Here, we compared the transcriptome profiles of rachis pre- and under cold stress in "Huoju" and "Ninghaibai". A total of 4,347 and 3,513 differentially expressed genes (DEGs) were detected in "Ninghaibai" and "Huoju", among which 223 and 166 were newly identified genes, respectively, most of them were functionally enriched in plant hormone signal transduction (Huoju: 142; Ninghaibai: 200), and there were higher plant hormone content and related DEG expression levels in "Huoju" than that of "Ninghaibai". Moreover, a total of 3,309 differentially expressed transcription factors (DETFs) were identified, and some DEGs and DETFs were screened to be subjected to co-expression network analysis based on the gene expression profile data. Some candidate DEGs, including UDP-glycosyltransferase (UGT), glycosyltransferase (GT), sugar phosphate/phosphate translocator (SPT), sugar transport protein (STP), proline-rich receptor-like protein kinase (PERK), and peroxidise (POD), were significantly affected by cold stress, and the expression level of these genes obtained from real-time quantitative RT-PCR was consistent with the pattern of transcriptome profile, which suggested that these genes might play the vital roles in cold resistance of loquat. Our results provide an invaluable resource for the identification of specific genes and TFs and help to clarify gene transcription during the cold stress response of loquat.

Comparative Transcriptomic Analysis of Differentially Expressed Transcripts Associated with Flowering Time of Loquat (Eriobotya japonica Lindl.)

Flowering is an important phenophase of plant species, however, knowledge about the regulatory mechanism controlling flowering cues in loquat is limited. To identify candidate genes regulating flowering time in loquat, we used RNA-Seq technology to conduct a comparative transcriptome analysis of differentiating apical buds collected from the early-flowering variety ‘Baiyu’ and the late-flowering variety ‘Huoju’. A total of 28,842 differentially expressed transcripts (DETs) were identified. Of these, 42 DETs controlled flowering time while 17 other DETs were associated with the ABA signaling pathway. Compared with those in ‘Huoju’, EjFT, EjFY, EjFLK, and EjCAL1-like were significantly upregulated in ‘Baiyu’. Moreover, transcripts of the ABA 8′-hydroxylases (EjABH2, EjABH4, and EjABH4-like2), the ABA receptors (EjPYL4/8), and the bZIP transcription factor EjABI5-like were upregulated in ‘Baiyu’ compared with ‘Huoju’. Hence, they might regulate loquat flowering time. There was no significant difference between ‘Baiyu’ and ‘Huoju’ in terms of IAA content. However, the ABA content was about ten-fold higher in the apical buds of ‘Baiyu’ than in those of ‘Huoju’. The ABA:IAA ratio sharply rose and attained a peak during bud differentiation. Thus, ABA is vital in regulating floral bud formation in loquat. The results of the present study help clarify gene transcription during loquat flowering.

Comparative transcriptome analysis reveals key genes potentially related to organic acid and sugar accumulation in loquat

Organic acids and sugars are the primary components that determine the quality and flavor of loquat fruits. In the present study, major organic acids, sugar content, enzyme activities, and the expression of related genes were analyzed during fruit development in two loquat cultivars, ’JieFangZhong’ (JFZ) and ’BaiLi’ (BL). Our results showed that the sugar content increased during fruit development in the two cultivars; however, the organic acid content dramatically decreased in the later stages of fruit development. The differences in organic acid and sugar content between the two cultivars primarily occured in the late stage of fruit development and the related enzymes showed dynamic changes in activies during development. Phosphoenolpyruvate carboxylase (PEPC) and mNAD malic dehydrogenase (mNAD-MDH) showed higher activities in JFZ at 95 days after flowering (DAF) than in BL. However, NADP-dependent malic enzyme (NADP-ME) activity was the lowest at 95 DAF in both JFZ and BL with BL showing higher activity compared with JFZ. At 125 DAF, the activity of fructokinase (FRK) was significantly higher in JFZ than in BL. The activity of sucrose synthase (SUSY) in the sucrose cleavage direction (SS-C) was low at early stages of fruit development and increased at 125 DAF. SS-C activity was higher in JFZ than in BL. vAI and sucrose phosphate synthase (SPS) activities were similar in the two both cultivars and increased with fruit development. RNA-sequencing was performed to determine the candidate genes for organic acid and sugar metabolism. Our results showed that the differentially expressed genes (DEGs) with the greated fold changes in the later stages of fruit development between the two cultivars were phosphoenolpyruvate carboxylase 2 (PEPC2), mNAD-malate dehydrogenase (mNAD-MDH), cytosolic NADP-ME (cyNADP-ME2), aluminum-activated malate transporter (ALMT9), subunit A of vacuolar H⁺-ATPase (VHA-A), vacuolar H⁺-PPase (VHP1), NAD-sorbitol dehydrogenase (NAD-SDH), fructokinase (FK), sucrose synthase in sucrose cleavage (SS-C), sucrose-phosphate synthase 1 (SPS1), neutral invertase (NI), and vacuolar acid invertase (vAI). The expression of 12 key DEGs was validated by quantitative reverese transcription PCR (RT-qPCR). Our findings will help understand the molecular mechanism of organic acid and sugar formation in loquat, which will aid in breeding high-quality loquat cultivars.

Dissection of the incidence and severity of purple spot physiological disorder in loquat fruit through a physiological and molecular approach

Loquat (Eriobotrya japonica) fruit marketability is affected by the incidence and severity of purple spot (PS), a pre-harvest physiological disorder showing an evident skin discoloration with depressed surface. Despite its impact in limiting the cultivation and economic potential of loquat, the etiology of this disorder is still poorly understood. To this end, our study aimed to investigate and disclose possible mechanisms underlying PS development. The intensity and severity of PS in three loquat cultivars ('Morphitiki', 'Karantoki' and 'Obusa') was phenotypically monitored during successive on-tree fruit developmental stages. 'Obusa' fruits harvested at commercial maturity stage showed the highest incidence of purple spot (58.6%), while 'Morphitiki' fruits did not show any symptoms. 'Karantoki' fruits demonstrated an intermediate response, with 31.3% of the fruit being affected. Thereafter, fruits with 30-50% PS severity were selected and used for further analysis; peel tissue was removed from both symptomatic and asymptomatic tissue of the same fruit for all examined cultivars. 'Karantoki' fruit with PS were characterized by the highest accumulation of total soluble sugars, sucrose, glucose and fructose contents, while the concentration of these primary metabolites was the lowest in asymptomatic fruit of 'Obusa', exception made for the sucrose. The incidence of PS was also transcriptionally investigated by assessing the mRNA profile of important genes involved in polyphenolic (PAL1, PAL2 and PPO1) and carbohydrate (CWI2, CWI3, SPS1, SPS2, NI2, NI3, SuSy, HXK, FRK and VI) pathway. The enhanced expression levels of CWI3 and VI genes in symptomatic fruit of the highly susceptible cultivar 'Obusa' highlight a cultivar-specific type of response. Notably, SuSy registered significantly suppressed levels in symptomatic tissue of both 'Obusa' and 'Karantoki'. To what extent PPO is associated with PS incidence and whether the etiology of the disorder can be assigned to an oxidative process triggered and coordinated by its action need to be further elucidated. The aforementioned genes are suggested to be further examined as potential markers towards a more sophisticated and informed characterization of purple spot detection in loquat fruit.

Single-molecule real-time sequencing of the full-length transcriptome of loquat under low-temperature stress

In this study, third-generation full-length (FL) transcriptome sequencing was performed of loquat using single-molecule real-time(SMRT) sequencing from the pooled cDNA of embryos of young loquat fruit under different low temperatures (three biological replicates for treatments of 1°C, -1°C, and -3°C, for 12 h or 24 h) and the control group(three biological replicates for treatments of room temperature), Illumina sequencing was used to correct FL transcriptome sequences. A total of 3 PacBio Iso-Seq libraries (1–2 kb, 2–3 kb and 3–6 kb) and 21 Illumina transcriptome libraries were constructed, a total of 13.41 Gb of clean reads were generated, which included 215,636 reads of insert (ROIs) and 121,654 FL, non-chimaric (FLNC) reads. Transcript clustering analysis of the FLNC reads revealed 76,586 consensus isoforms, and a total of 12,520 high-quality transcript sequences corrected with non-FL sequences were used for subsequent analysis. After the redundant reads were removed, 38,435 transcripts were obtained. A total of 27,905 coding DNA sequences (CDSs) were identified, and 407 long non-coding RNAs (lncRNAs) were ultimately predicted. Additionally, 24,832 simple sequence repeats (SSRs) were identified, and a total of 1,295 alternative splicing (AS) events were predicted. Furthermore, 37,993 transcripts were annotated in eight functional databases. This is the first study to perform SMRT sequencing of the FL transcriptome of loquat. The obtained transcriptomic data are conducive for further exploration of the mechanism of loquat freezing injury and thus serve as an important theoretical basis for generating new loquat material and for identifying new ways to improve loquat cold resistance.

An Integrative Analysis of Transcriptome, Proteome and Hormones Reveals Key Differentially Expressed Genes and Metabolic Pathways Involved in Flower Development in Loquat

Flower development is a vital developmental process in the life cycle of woody perennials, especially fruit trees. Herein, we used transcriptomic, proteomic, and hormone analyses to investigate the key candidate genes/proteins in loquat (Eriobotrya japonica) at the stages of flower bud differentiation (FBD), floral bud elongation (FBE), and floral anthesis (FA). Comparative transcriptome analysis showed that differentially expressed genes (DEGs) were mainly enriched in metabolic pathways of hormone signal transduction and starch and sucrose metabolism. Importantly, the DEGs of hormone signal transduction were significantly involved in the signaling pathways of auxin, gibberellins (GAs), cytokinin, ethylene, abscisic acid (ABA), jasmonic acid, and salicylic acid. Meanwhile, key floral integrator genes FLOWERING LOCUS T (FT) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) and floral meristem identity genes SQUAMOSA PROMOTER BINDING LIKE (SPL), LEAFY (LFY), APETALA1 (AP1), and AP2 were significantly upregulated at the FBD stage. However, key floral organ identity genes AGAMOUS (AG), AP3, and PISTILLATA (PI) were significantly upregulated at the stages of FBE and FA. Furthermore, transcription factors (TFs) such as bHLH (basic helix-loop-helix), NAC (no apical meristem (NAM), Arabidopsis transcription activation factor (ATAF1/2) and cup-shaped cotyledon (CUC2)), MYB_related (myeloblastosis_related), ERF (ethylene response factor), and C2H2 (cysteine-2/histidine-2) were also significantly differentially expressed. Accordingly, comparative proteomic analysis of differentially accumulated proteins (DAPs) and combined enrichment of DEGs and DAPs showed that starch and sucrose metabolism was also significantly enriched. Concentrations of GA3 and zeatin were high before the FA stage, but ABA concentration remained high at the FA stage. Our results provide abundant sequence resources for clarifying the underlying mechanisms of the flower development in loquat.

Transcriptome analysis provides insights into the regulation of metabolic processes during postharvest cold storage of loquat (Eriobotrya japonica) fruit

Loquat (Eriobotrya japonica) fruit accumulates lignin during postharvest storage under chilling conditions (0 °C), while low-temperature conditioning (LTC; 5 °C for 6 days followed by transfer to 0 °C) or heat treatment (HT; 40 °C for 4 h followed by transfer to 0 °C) can alleviate lignification. Here we compared transcriptome profiles of loquat fruit samples under LTC or HT to those stored at 0 °C at five time points from day 1 to day 8 after treatment. High-throughput transcriptome sequences were de novo assembled into 53,319 unique transcripts with an N50 length of 1306 bp. A total of 2235 differentially expressed genes were identified in LTC, and 1020 were identified in HT compared to 0 °C. Key genes in the lignin biosynthetic pathway, including EjPAL2, EjCAD1, EjCAD3, 4CL, COMT, and HCT, were responsive to LTC or HT treatment, but they showed different expression patterns during the treatments, indicating that different structural genes could regulate lignification at different treatment stages. Coexpression network analysis showed that these candidate biosynthetic genes were associated with a number of transcription factors, including those belonging to the AP2, MYB, and NAC families. Gene ontology (GO) enrichment analysis of differentially expressed genes indicated that biological processes such as stress responses, cell wall and lignin metabolism, hormone metabolism, and metal ion transport were significantly affected under LTC or HT treatment when compared to 0 °C. Our analyses provide insights into transcriptome responses to postharvest treatments in loquat fruit.

Metabolic Dynamics During Loquat Fruit Ripening and Postharvest Technologies

Loquat is an important fruit widely cultivated worldwide with high commercial value. During loquat fruit development, ripening, and storage, many important metabolites undergo dramatic changes, resulting in accumulation of a diverse mixture of nutrients. Given the value of loquat fruit, significant progresses have been achieved in understanding the metabolic changes during fruit ripening and storage, as well as postharvest technologies applied in loquat fruit in recent years. The objective of the present review is to summarize currently available knowledge and provide new references for improving loquat fruit quality.

De novo Transcriptome Profiling of Flowers, Flower Pedicels and Pods of Lupinus luteus (Yellow Lupine) Reveals Complex Expression Changes during Organ Abscission

Yellow lupine (Lupinus luteus L., Taper c.), a member of the legume family (Fabaceae L.), has an enormous practical importance. Its excessive flower and pod abscission represents an economic drawback, as proper flower and seed formation and development is crucial for the plant's productivity. Generative organ detachment takes place at the basis of the pedicels, within a specialized group of cells collectively known as the abscission zone (AZ). During plant growth these cells become competent to respond to specific signals that trigger separation and lead to the abolition of cell wall adhesion. Little is known about the molecular network controlling the yellow lupine organ abscission. The aim of our study was to establish the divergences and similarities in transcriptional networks in the pods, flowers and flower pedicels abscised or maintained on the plant, and to identify genes playing key roles in generative organ abscission in yellow lupine. Based on de novo transcriptome assembly, we identified 166,473 unigenes representing 219,514 assembled unique transcripts from flowers, flower pedicels and pods undergoing abscission and from control organs. Comparison of the cDNA libraries from dropped and control organs helped in identifying 1,343, 2,933 and 1,491 differentially expressed genes (DEGs) in the flowers, flower pedicels and pods, respectively. In DEG analyses, we focused on genes involved in phytohormonal regulation, cell wall functioning and metabolic pathways. Our results indicate that auxin, ethylene and gibberellins are some of the main factors engaged in generative organ abscission. Identified 28 DEGs common for all library comparisons are involved in cell wall functioning, protein metabolism, water homeostasis and stress response. Interestingly, among the common DEGs we also found an miR169 precursor, which is the first evidence of micro RNA engaged in abscission. A KEGG pathway enrichment analysis revealed that the identified DEGs were predominantly involved in carbohydrate and amino acid metabolism, but some other pathways were also targeted. This study represents the first comprehensive transcriptome-based characterization of organ abscission in L. luteus and provides a valuable data source not only for understanding the abscission signaling pathway in yellow lupine, but also for further research aimed at improving crop yields.