Transcriptome analysis unravels spatiotemporal modulation of phytohormone-pathway expression underlying gibberellin-induced parthenocarpic fruit set in San Pedro-type fig (Ficus carica L.)

Parthenocarpy, a pollination-independent fruit set mechanism to ensure yield stability

Fruit development is essential for flowering plants' reproduction and a significant food source. Climate change threatens fruit yields due to its impact on pollination and fertilization processes, especially vulnerable to extreme temperatures, insufficient light, and pollinator decline. Parthenocarpy, the development of fruit without fertilization, offers a solution, ensuring yield stability in adverse conditions and enhancing fruit quality. Parthenocarpic fruits not only secure agricultural production but also exhibit improved texture, appearance, and shelf life, making them desirable for food processing and other applications. Recent research unveils the molecular mechanisms behind parthenocarpy, implicating transcription factors (TFs), noncoding RNAs, and phytohormones such as auxin, gibberellin (GA), and cytokinin (CK). Here we review recent findings, construct regulatory models, and identify areas for further research.

Transcriptomic and Hormone Analyses Provide Insight into the Regulation of Axillary Bud Outgrowth of Eucommia ulmoides Oliver

An essential indicator of Eucommia ulmoides Oliver (E. ulmoides) is the axillary bud; the growth and developmental capacity of axillary buds could be used to efficiently determine the structural integrity of branches and plant regeneration. We obtained axillary buds in different positions on the stem, including upper buds (CK), tip buds (T1), and bottom buds (T2), which provided optimal materials for the study of complicated regulatory networks that control bud germination. This study used transcriptomes to analyze the levels of gene expression in three different types of buds, and the results showed that 12,131 differentially expressed genes (DEGs) were discovered via the pairwise comparison of transcriptome data gathered from CK to T2, while the majority of DEGs (44.38%) were mainly found between CK and T1. These DEGs were closely related to plant hormone signal transduction and the amino acid biosynthesis pathway. We also determined changes in endogenous hormone contents during the process of bud germination. Interestingly, except for indole-3-acetic acid (IAA) content, which showed a significant upward trend (p < 0.05) in tip buds on day 4 compared with day 0, the other hormones showed no significant change during the process of germination. Then, the expression patterns of genes involved in IAA biosynthesis and signaling were examined through transcriptome analysis. Furthermore, the expression levels of genes related to IAA biosynthesis and signal transduction were upregulated in tip buds. Particularly, the expression of the IAA degradation gene Gretchen Hagen 3 (GH3.1) was downregulated on day 4, which may support the concept that endogenous IAA promotes bud germination. Based on these data, we propose that IAA synthesis and signal transduction lead to morphological changes in tip buds during the germination process. On this basis, suggestions to improve the efficiency of the production and application of E. ulmoides are put forward to provide guidance for future research.

Genome-wide identification and characterization of parthenocarpic fruit set-related gene homologs in cucumber (Cucumis sativus L.)

Cucumber (Cucumis sativus L.), a major horticultural crop, in the family Cucurbitaceae is grown and consumed globally. Parthenocarpy is an ideal trait for many fruit and vegetables which produces seedless fruit desired by consumers. The seedlessness occurs when fruit develops without fertilization which can be either natural or induced. So far, a limited number of genes regulating parthenocarpic fruit set have been reported in several fruit or vegetable crops, most of which are involved in hormone biosynthesis or signalling. Although parthenocarpic cucumber has been widely used in commercial production for a long time; its genetic basis is not well understood. In this study, we retrieved thirty five parthenocarpy fruit-set related genes (PRGs) from bibliomic data in various plants. Thirty-five PRG homologs were identified in the cucumber genome via homology-based search. An in silico analysis was performed on phylogenetic tree, exon-intron structure, cis-regulatory elements in the promoter region, and conserved domains of their deduced proteins, which provided insights into the genetic make-up of parthenocarpy-related genes in cucumber. Simple sequence repeat (SSR) sequences were mined in these PRGs, and 31 SSR markers were designed. SSR genotyping identified three SSRs in two polymorphic genes. Quantitative real-time PCR of selected genes was conducted in five cucumber lines with varying degrees of parthenocarpic fruit set capacities, which revealed possible association of their expression with parthenocarpy. The results revealed that homologs CsWD40 and CsPIN-4 could be considered potential genes for determination of parthenocarpy as these genes showed parental polymorphism and differential gene expression in case of parthenocarpic and non-parthenocarpic parents.

AP2/ERF genes associated with superfast fig (Ficus carica L.) fruit ripening

Fig fruits have significant health value and are culturally important. Under suitable climatic conditions, fig fruits undergo a superfast ripening process, nearly doubling in size, weight, and sugar content over three days in parallel with a sharp decrease in firmness. In this study, 119 FcAP2/ERF genes were identified in the fig genome, namely 95 ERFs, 20 AP2s, three RAVs, and one soloist. Most of the ERF subfamily members (76) contained no introns, whereas the majority of the AP2 subfamily members had at least two introns each. Three previously published transcriptome datasets were mined to discover expression patterns, encompassing the fruit peel and flesh of the 'Purple Peel' cultivar at six developmental stages; the fruit receptacle and flesh of the 'Brown Turkey' cultivar after ethephon treatment; and the receptacle and flesh of parthenocarpic and pollinated fruits of the 'Brown Turkey' cultivar. Eighty-three FcAP2/ERFs (68 ERFs, 13 AP2s, one RAV, and one soloist) were expressed in the combined transcriptome dataset. Most FcAP2/ERFs were significantly downregulated (|log₂(fold change) | ≥ 1 and p-adjust < 0.05) during both normal fruit development and ethephon-induced accelerated ripening, suggesting a repressive role of these genes in fruit ripening. Five significantly downregulated ERFs also had repression domains in the C-terminal. Seven FcAP2/ERFs were identified as differentially expressed during ripening in all three transcriptome datasets. These genes were strong candidates for future functional genetic studies to elucidate the major FcAP2/ERF regulators of the superfast fig fruit ripening process.

Proteome Analysis of Vacuoles Isolated from Fig (Ficus carica L.) Flesh during Fruit Development

Fruit flesh cell vacuoles play a pivotal role in fruit growth and quality formation. In the present study, intact vacuoles were carefully released and collected from protoplasts isolated from flesh cells at five sampling times along fig fruit development. Label-free quantification and vacuole proteomic analysis identified 1,251 proteins, 1,137 of which were recruited as differentially abundant proteins (DAPs) by fold change ≥ 1.5, P < 0.05. DAPs were assigned to 10 functional categories; among them, 238, 186, 109, 93 and 90 were annotated as metabolism, transport proteins, membrane fusion or vesicle trafficking, protein fate and stress response proteins, respectively. Decreased numbers of DAPs were uncovered along fruit development. The overall changing pattern of DAPs revealed two major proteome landscape conversions in fig flesh cell vacuoles: the first occurred when fruit developed from late-stage I to mid-stage II, and the second occurred when the fruit started ripening. Metabolic proteins related to glycosidase, lipid and extracellular proteins contributing to carbohydrate storage and vacuole expansion, and protein-degrading proteins determining vacuolar lytic function were revealed. Key tonoplast proteins contributing to vacuole expansion, cell growth and fruit quality formation were also identified. The revealed comprehensive changes in the vacuole proteome during flesh development were compared with our previously published vacuole proteome of grape berry. The information expands our knowledge of the vacuolar proteome and the protein basis of vacuole functional evolution during fruit development and quality formation.

Genome-wide analysis of JAZ family genes expression patterns during fig (Ficus carica L.) fruit development and in response to hormone treatment

BACKGROUND

Jasmonate-ZIM domain (JAZ) repressors negatively regulate signal transduction of jasmonates, which regulate plant development and immunity. However, no comprehensive analysis of the JAZ gene family members has been done in the common fig (Ficus carica L.) during fruit development and hormonal treatment.

RESULTS

In this study, 10 non-redundant fig JAZ family genes (FcJAZs) distributed on 7 chromosomes were identified in the fig genome. Phylogenetic and structural analysis showed that FcJAZ genes can be grouped into 5 classes. All the classes contained relatively complete TIFY and Jas domains. Yeast two hybrid (Y2H) results showed that all FcJAZs proteins may interact with the identified transcription factor, FcMYC2. Tissue-specific expression analysis showed that FcJAZs were highly expressed in the female flowers and roots. Expression patterns of FcJAZs during the fruit development were analyzed by RNA-Seq and qRT-PCR. The findings showed that, most FcJAZs were significantly downregulated from stage 3 to 5 in the female flower, whereas downregulation of these genes was observed in the fruit peel from stage 4 to 5. Weighted-gene co-expression network analysis (WGCNA) showed the expression pattern of FcJAZs was correlated with hormone signal transduction and plant-pathogen interaction. Putative cis-elements analysis of FcJAZs and expression patterns of FcJAZs which respond to hormone treatments revealed that FcJAZs may regulate fig fruit development by modulating the effect of ethylene or gibberellin.

CONCLUSIONS

This study provides a comprehensive analysis of the FcJAZ family members and provides information on FcJAZs contributions and their role in regulating the common fig fruit development.

Transcriptome Profiles of Leaves and Roots of Goldenrain Tree (Koelreuteria paniculata Laxm.) in Response to Cadmium Stress

Cadmium (Cd) pollution is a widespread environmental problem. In this study, we explored the transcriptome and biochemical responses of goldenrain tree (Koelreuteria paniculata Laxm.) leaves and roots to Cd stress. Leaf and root growth decreased substantially under Cd stress (50 mg/L CdCl2), but leaf and root antioxidant mechanisms were significantly activated. In RNA-seq analysis, roots treated with 25 mg/L CdCl2 featured enriched GO terms in cellular components related to intracellular ribonucleoprotein complex, ribonucleoprotein complex, and macromolecular complex. In leaves under Cd stress, most differentially expressed genes were enriched in the cellular component terms intrinsic component of membrane and membrane part. Weighted gene co-expression network analysis and analysis of module-trait relations revealed candidate genes associated with superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities and malondialdehyde (MDA). Ten transcription factors responded to Cd stress expression, including those in C2H2, MYB, WRKY, and bZIP families. Transcriptomic analysis of goldenrain tree revealed that Cd stress rapidly induced the intracellular ribonucleoprotein complex in the roots and the intrinsic component of membrane in the leaves. The results also indicate directions for further analyses of molecular mechanisms of Cd tolerance and accumulation in goldenrain tree.

AGAMOUS Gene as a New Sex-Identification Marker in Fig (Ficus carica L.) Is More Efficient Than RAN1

Fig is an ancient gynodioecious fruit tree with females for commercial fruit production and hermaphrodites (males) sometimes used as pollen providers. An early sex-identification method would improve breeding efficiency. Three AGAMOUS (AG) genes were recruited from the Ficus carica genome using AG sequences from Ficus microcarpa and Ficus hispida. FcAG was 5230 bp in length, with 7 exons and 6 introns, and a 744-bp coding sequence. The gene was present in both female and male fig genomes, with a 15-bp deletion in the 7th exon. The other two AG genes (FcAG2-Gall_Stamen and FcAG3-Gall_Stamen) were male-specific, without the 15-bp deletion (759-bp coding sequence), and were only expressed in the gall and stamen of the male fig fruit. Using the deletion as the forward primer (AG-Marker), male plants were very efficiently identified by the presence of a 146-bp PCR product. The previously reported fig male and female polymorphism gene RESPONSIVE-TO-ANTAGONIST1 (RAN1) was also cloned and compared between male and female plants. Fifteen SNPs were found in the 3015-bp protein-coding sequence. Among them, 12 SNPs were identified as having sex-differentiating capacity by checking the sequences of 27 known male and 24 known female cultivars. A RAN1-Marker of 608 bp, including 6 SNPs, was designed, and a PCR and sequencing-based method was verified with 352 fig seedlings from two hybrid populations. Our results confirmed that the newly established AG-Marker is as accurate as the RAN1-Marker, and provide new clues to understanding Ficus sex determination.

Papain-Like Cysteine Protease Gene Family in Fig (Ficus carica L.): Genome-Wide Analysis and Expression Patterns

The papain-like cysteine proteases (PLCPs) are the most abundant family of cysteine proteases in plants, with essential roles in biotic/abiotic stress responses, growth and senescence. Papain, bromelain and ficin are widely used in food, medicine and other industries. In this study, 31 PLCP genes (FcPCLPs) were identified in the fig (Ficus carica L.) genome by HMM search and manual screening, and assigned to one of nine subfamilies based on gene structure and conserved motifs. SAG12 and RD21 were the largest subfamilies with 10 and 7 members, respectively. The FcPCLPs ranged from 1,128 to 5,075 bp in length, containing 1-10 introns, and the coding sequence ranged from 624 to 1,518 bp, encoding 207-505 amino acids. Subcellular localization analysis indicated that 24, 2, and 5 PLCP proteins were targeted to the lysosome/vacuole, cytoplasm and extracellular matrix, respectively. Promoter (2,000 bp upstream) analysis of FcPLCPs revealed a high number of plant hormone and low temperature response elements. RNA-seq revealed differential expression of 17 FcPLCPs in the inflorescence and receptacle, and RD21 subfamily members were the major PLCPs expressed in the fruit; 16 and 5 FcPLCPs responded significantly to ethylene and light, respectively. Proteome analyses revealed 18 and 5 PLCPs in the fruit cell soluble proteome and fruit latex, respectively. Ficins were the major PLCP in fig fruit, with decreased abundance in inflorescences, but increased abundance in receptacles of commercial-ripe fruit. FcRD21B/C and FcALP1 were aligned as the genes encoding the main ficin isoforms. Our study provides valuable multi-omics information on the FcPLCP family and lays the foundation for further functional studies.

Transcriptome analysis reveals the effects of grafting on sweetpotato scions during the full blooming stages

BACKGROUND

Sweetpotato (Ipomoea batatas) is a hexaploid plant and generally most genotypes do not flower at all in sub-tropics. Heterografting was carried out between sweetpotato cultivar 'Xushu 18' and Japanese morning glory (Ipomoea nil). With sweetpotato as 'scion' and I. nil as 'rootstock', sweetpotato was induced flowering in the autumn. However, little is known about the molecular mechanisms underlying sweetpotato responses to grafting, especially during the full blooming stages.

OBJECTIVES

To investigate the poorly understood molecular responses underlying the grafting-induced phenotypic processes in sweetpotato at full anthesis.

METHODS

In this study, to explore the transcriptome diversity and complexity of sweetpotato, PacBio Iso-Seq and Illumina RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression of five tissues: scion flowers (SF), scion leaves (SL), scion stems (SS), own-rooted leaves (OL) and own-rooted stems (OS).

RESULTS

A total of 138,151 transcripts were generated with an average length of 2255 bp, and more than 72% (100,396) of the transcripts were full-length. During full blooming, to examine the difference in gene expression of sweetpotato under grafting and natural growth conditions, 7905, 7795 and 15,707 differentially expressed genes were detected in pairwise comparisons of OS versus SS, OL versus SL and SL versus SF, respectively. Moreover, differential transcription of genes associated with anthocyanin biosynthesis, light pathway and photosynthesis, ethylene signal transduction pathway was observed in scion responses to grafting.

CONCLUSION

Our study is useful in understanding the molecular basis of grafting-induced flowering in grafted sweetpotatoes, and will lay a foundation for further research on sweetpotato breeding in the future.