REVIEW article

Front. Plant Sci., 16 January 2026

Sec. Plant Metabolism and Chemodiversity

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1750317

Transcriptomics in solanaceous crop improvement: advances and opportunities

  • 1. School of Horticulture, Anhui Agricultural University, Hefei, China

  • 2. Department of Plant Pathology, CCS Haryana Agricultural University, Hisar, India

  • 3. Department of Horticulture, Virginia Polytechnic Institute and State University, Blacksburg, VA, United States

  • 4. Horticultural Sciences Department, University of Florida, Gainesville, FL, United States

  • 5. Department of Plant and Soil Sciences, Mississippi State University, Starkville, MS, United States

  • 6. School of Life Science, Huaibei Normal University, Huaibei, Anhui, China

Article metrics

View details

795

Views

100

Downloads

Abstract

The solanaceae family includes several economically important crops. Tomato (Solanum lycopersicum), potato (S. tuberosum), eggplant (S. melongena), and pepper (Capsicum annuum) are among the most consumed vegetables worldwide. Over the past few decades, consumers have often expressed dissatisfaction with declining fruit quality, especially flavor, while growers’ profits remain tied primarily to yield and stress resistance. Thus, reconciling these divergent demands provides a clear roadmap for improving solanaceous crops. To completely understand molecular mechanisms underlying biological processes associated with fruit quality, yield, and stresses, there is a need to focus on gene expression at the RNA level. With the availability of next-generation sequencing (NGS) technologies, transcriptomics platforms have helped to improve our knowledge about RNA-based gene regulatory networks. The current review discusses the recent literature on transcriptomics techniques and their progress status in major solanaceae crops across important agronomic traits. In addition to the current improvements, we have also discussed future perspectives of transcriptomics in genetic engineering.

1 Introduction

Increasing temperatures and shrinking water availability in major fruit and vegetable production regions threatens the availability of nutritious food globally (Cammarano et al., 2022; Wang et al., 2017). Among horticultural crops, the Solanaceae family represents a strategic target for genetic improvement due to its nutritional importance and genetic diversity (Kaur et al., 2023). With the global population increasing to an estimated 9.7 billion by 2050, it is necessary to accelerate genetic gain in highly nutritious Solanaceae crops (Keating et al., 2014; Sandhu et al., 2022; Kantar et al., 2025).

The Solanaceae family encompasses ~100 genera and 2,500 species (Olmstead et al., 2008) among which tomato (S. lycopersicum), pepper (C. annuum L.), potato (S. tuberosum), and eggplant (S. melongena) are significant contributors to nutrition availability. Furthermore, secondary metabolites from these crops are used in pharmaceuticals, flavoring agents, and fragrances (Fiesel et al., 2022). However, these solanaceous crops face a critical trade-off between consumer preferences and commercial production objectives. While consumers prioritize fruit quality attributes like flavor, freshness, and appearance (Zhang et al., 2023), enhanced tolerance to heat, drought, and biotic stresses is of utmost importance to farmers. This diversity in trait priorities necessitates a multi-trait breeding strategy that simultaneously addresses yield, stress resilience and fruit quality traits (Kaur et al., 2023; Singh et al., 2025).

To bridge the gap between consumer demands and production constraints while adapting to climate change, solanaceous crop breeding can exploit the knowledge derived from molecular and genetic mechanisms of tolerance to heat, drought, biotic stresses and flavor. Next-generation sequencing (NGS) and transcriptomics have enhanced our ability to characterize genome-wide expression patterns to construct functional gene regulatory networks for these complex traits. Transcriptomic studies have successfully identified candidate genes conferring resistance to biotic stresses (Meng et al., 2021) abiotic stresses (Altaf et al., 2022), enhanced flavor compounds (Gao et al., 2019), improved post-harvest quality (Cruz-Mendívil et al., 2015; Popovsky-Sarid et al., 2017; Wahrenburg et al., 2021) and nutritional content (Kim et al., 2014) (Figure 1). In solanaceous crops, transcriptomic approaches are particularly valuable for dissecting complex, polygenic traits where phenotypic selection alone proves insufficient. Developmental processes and stress responses are governed by intricate signaling cascades that operate across temporal and spatial scales, with coordinated gene expression patterns ultimately determining phenotypic outcomes. However, existing reviews have focused on individual species or single trait categories, while a cross-species, multi-trait synthesis examining transcriptomic insights across stress tolerance, yield, and quality traits remains absent in solanaceous crops. This review paper addresses the aforementioned knowledge gap and provides a snapshot of recent advances in transcriptomics of solanaceous crop improvement emphasizing translational insights for multi-trait breeding strategies.

Figure 1

Diagram showing how transcriptomics is used to analyze abiotic stress, biotic stress, yield, and quality traits. Central arrow points to quantifying differential gene expression responses. Below, it highlights gene ontology, pathway analysis, heatmaps, global profiling, quantification, comparative transcriptomic analysis, alignment, clustering, and integration with other “omics.”

An overview of the utilization of transcriptomics techniques in economically important traits in solanaceous crops.

2 Timeline of the development of transcriptomics technology in Solanaceae crops

The integration of transcriptomics with genomic resources has transformed candidate gene identification in Solanaceae breeding. By integrating gene expression profiles with genomic regions identified through genome-wide association studies and quantitative trait loci (QTL) mapping, breeders can prioritize causal genes from among multiple candidates and develop functional markers that capture allelic variation affecting expression levels. This expression quantitative trait loci (eQTL) approach enables more precise marker-assisted selection and accelerates the validation of genes controlling complex traits. Multiple technological platforms have emerged to interrogate transcriptomes, each offering distinct advantages for specific applications. These range from hybridization-based microarrays to high-throughput RNA sequencing (RNA-seq), single-cell transcriptomics, and specialized methods for profiling small regulatory RNAs. We review these approaches sequentially (Figure 2), emphasizing their applications in solanaceous crop improvement and the insights each technology has provided for breeding programs.

Figure 2

Timeline illustrating significant milestones in genome sequencing: 1965 marks the sequencing of the first tRNA, followed by the discovery of reverse transcriptase in 1970, Northern and Sanger sequencing in 1977, the polymerase chain reaction in 1984, and automation of Sanger sequencing in 1987. Later milestones include sequencing the Arabidopsis genome in 2000, potato in 2011, tomato in 2012, pepper in 2013, and eggplant in 2014.

Timeline of the major advancements contributing to transcriptomics emergence as a tool in solanaceous crop improvement.

Early transcriptomic approaches were constrained by their single-gene focus. Northern blotting enabled detection and quantification of specific transcripts through gel separation and probe hybridization (Alwine et al., 1977), while reverse transcription PCR and cDNA library construction allowed targeted gene expression analysis (Mullis et al., 1986). However, these methods could not capture the genome-wide expression dynamics necessary for understanding complex regulatory networks underlying multi-trait phenotypes.

The need for simultaneous profiling of thousands of genes drove the development of microarray platforms in the late 1990s. Microarrays operate by hybridizing fluorescently labeled RNA samples to gene-specific probes arrayed at known positions on solid support. Early platforms utilized cDNA probes derived from expressed sequence tags (ESTs), while later iterations employed synthesized oligonucleotides; either short probes with multiple sequences per gene (25-mers on Affymetrix chips) or longer probes with single or few representatives (>60-mers) (Yamada et al., 2003). In Solanaceae, microarrays enabled the first genome-scale expression studies, including identification of ripening-related transcripts in tomato (Karlova et al., 2011). Despite these advances, microarrays suffered from substantial limitations: 40-60% of genes showed no detectable expression due to high background noise, expression distributions were highly skewed toward low-abundance transcripts, and probe design required prior sequence knowledge (Klebanov and Yakovlev, 2007; Zhao et al., 2014).

To address microarray limitations while maintaining high throughput, tag-based sequencing methods emerged, including serial analysis of gene expression (Matsumura et al., 1999), cap analysis of gene expression (Agarwal et al., 2014), and massively parallel signature sequencing (Meyers et al., 2004). These approaches sequenced short tags derived from defined transcript positions, eliminating hybridization bias and enabling discovery of novel transcripts without prior sequence knowledge. However, short tag lengths (typically 14–21 nucleotides) prevented discrimination of alternative splice isoforms and paralogs with high sequence similarity highlighting critical limitations for studying gene families in polyploid crops like potato, hence ultimately restricting widespread adoption (Wang et al., 2009).

RNA-seq fundamentally transformed transcriptome profiling by replacing hybridization-based detection with direct sequencing and digital quantification of cDNA molecules derived from RNA samples (Wang et al., 2009). This approach eliminates probe design requirements and hybridization biases, enabling unbiased detection of all transcripts including protein-coding mRNAs and diverse non-coding RNAs (lncRNA, miRNA, circRNA, rRNA, and tRNA) without prior sequence knowledge. By counting sequence reads aligned to genomic or transcript references, RNA-seq provides absolute quantification across a dynamic range exceeding five orders of magnitude, detecting both highly abundant transcripts and rare, lowly expressed genes that evade microarray detection. This sensitivity has proven particularly valuable in Solanaceae crops for capturing expression of biosynthetic pathway genes (often expressed at low levels) controlling flavor compound production in tomato and pepper (Gao et al., 2019).

Beyond improved sensitivity, RNA-seq enables three critical capabilities for dissecting complex traits: (1) Discovery of novel transcripts and unannotated genes through de novo assembly, expanding the functional repertoire of crop genomes (Liu et al., 2015); (2) Characterization of alternative splicing events and isoform diversity, revealing how single genes generate multiple protein products with distinct functions (Keller et al., 2017; Yang et al., 2020c); and (3) Detection of post-transcriptional modifications including RNA editing events that alter protein sequence (Zouine et al., 2014; Amiryousefi et al., 2018). In Solanaceae, these capabilities have been leveraged to identify candidate genes controlling fruit ripening and quality traits in tomato (Zhu et al., 2015), characterize the transcriptional architecture of heat and drought responses across multiple genotypes in pepper (Kim et al., 2024) and profile pathogen-responsive transcriptional networks in eggplant (Chen et al., 2018). The introduction of short-read sequencing platforms particularly Illumina’s sequencing-by-synthesis (Egan et al., 2012) and Ion Torrent’s semiconductor-based sequencing (Merriman et al., 2012) dramatically reduced per-sample costs from thousands to hundreds of dollars, democratizing access and enabling large-scale expression profiling across diverse genotypes, environments, and developmental stages.

Two technological advances have further expanded RNA-seq’s applications in crop genomics. First, long-read sequencing platforms including Pacific Biosciences (PacBio) isoform sequencing (Iso-Seq) and Oxford Nanopore direct RNA sequencing sequence full-length transcripts in single reads (>1–10 kb), eliminating the need for computational isoform reconstruction and enabling definitive characterization of alternative splicing, fusion transcripts, and allele-specific expression in heterozygous and polyploid genomes like potato. Second, single-cell RNA-seq (scRNA-seq) resolves transcriptomes of individual cells within heterogeneous tissues, revealing cell-type-specific expression programs that are masked in bulk tissue analyses (Stuart and Satija, 2019; Aldridge and Teichmann, 2020). In plants, scRNA-seq has begun elucidating cell-type-specific responses to developmental cues and environmental stresses (Seyfferth et al., 2021), though technical challenges remain: low transcript capture efficiency (5-20%), high per-cell costs limiting sample sizes, ambient RNA contamination from lysed cells, and critically for plant tissues with thick cell walls difficulties in generating high-quality single-cell suspensions. Emerging spatial transcriptomics technologies address this last limitation by preserving tissue architecture while profiling spatially resolved expression, offering promise for understanding fruit development and ripening processes where spatial gradients in gene expression determine quality traits.

3 Transcriptomics of fruit quality traits and regulatory networks

Solanaceous fruits are an essential source of phytochemicals and secondary metabolites that influence plant defense, fruit quality, UV protection, and signaling pathways (Carvalho Lemos et al., 2019). Transcriptomic approaches have proven particularly powerful for dissecting these complex biosynthetic pathways, where multi-enzyme cascades, tissue-specific expression, and coordinated regulation cannot be inferred from genome sequences alone. The quality traits of several solanaceous crops that have been examined using transcriptomics are listed in Table 1.

Table 1

Crop Technique Quality trait Characteristics References
eggplant qRT-PCR anthocyanin and CGA transcript abundance in different tissues (Docimo et al., 2016)
eggplant RNA-seq anthocyanin screening for DEGs from transcriptomic data of eggplant peel during the transition from darkness to varying durations of light exposure (Li et al., 2018a)
eggplant RNA-seq browning potential DEGs associated with browning between the two eggplant cultivars were screened (Liu et al., 2021b)
pepper qRT-PCR AsA gene expression of biosynthetic genes is inversely related to AsA (Alós et al., 2013)
pepper RNA-seq capsaicinoids unigenes (transcripts) were assembled by RNA-seq for the mixture of placenta and pericarp of pungent pepper, and new candidate genes involved in capsaicinoid biosynthesis were predicted (Liu et al., 2013)
pepper RNA-seq capsaicinoids and AsA RNA-Seq was used to obtain transcriptomes of whole Serrano-type chili pepper fruits, and genes expression related to capsaicinoid and ascorbic acid biosynthesis were analyzed (Martínez-López et al., 2014)
pepper RNA-seq quercetin and derivatives biosynthetic pathway of quercetin and derivatives from sweet pepper fruits were analysis, and expression profile of genes was analysed (Guevara et al., 2021)
pepper RNA-seq flavonoid, phenylalanine, tyrosine, tryptophan and phenylpropanoid combined transcriptome and metabolome analysis showed that the genes and metabolites that responded to UV-C were mainly involved in flavonoid, phenylalanine, tyrosine, tryptophan and phenylpropanoid biosynthesis (Ma et al., 2021a)
pepper RNA-seq fruit development and ripening(shape, size, ethylene and secondary metabolites biosynthesis) transcriptomic analysis on immature stage and mature stage of a wild and a cultivated pepper variety, DEGs related to fruit ripening were futher screened (Razo-Mendivil et al., 2021)
pepper RNA-seq ascorbate peroxidase expression analysis of CaAPX genes by RNA-Seq at different stages of pepper fruit ripening (González-Gordo et al., 2022)
pepper multi-omics flavonoids and carotenoids DEGs were obtained and used for gene co-expression network analysis (WGCNA) (Liu et al., 2020)
pepper multi-omics methyl jasmonate model of methyl jasmonate-regulated postharvest chilling injury in green pepper fruits was studied by integrating transcriptomics, proteomics, and metabolomics (Fu et al., 2022)
potato microarray carotenoids differentially expressed microarray probes present within the major QTL for total tuber carotenoid content (Campbell et al., 2014)
potato multi-omics soluble sugars combined transcriptomic and proteomic analysis was conducted on potato tubers to investigate the mechanism of cold responses during postharvest storage (Lin et al., 2019)
potato multi-omics anthocyanins genomic data and the transcriptome data was to identify the crucial major genes that affect anthocyanin biosynthesis (Tengkun et al., 2019)
potato multi-omics browning postcutting browning were evaluated with transcriptomics and untargeted metabolomics to investigate the molecular mechanism (Qiao et al., 2022)
tobacco microarray carotenoids microarray was adapted for transcriptome comparison of tobacco leaves derived from three cultivated regions (Lei et al., 2013)
tobacco transcriptome cadmium(cd) DEGs obtained in the root and leaf samples of two Nicotiana species under the same Cd stress concentration (Zhang et al., 2021c)
tomato RNA-seq chilling DEGs responsed to chilling (Zhang et al., 2016)
tomato RNA-seq malate accumulation DEGs determined by RNA-seq in overexpression line (Ye et al., 2017)
tomato RNA-seq ethylene synthesis, signal transduction and carotenoid biosynthesis RNA-Seq was used to identify genes that were differentially expressed in 1-MCP-treated tomato fruits (Bobokalonov et al., 2018)
tomato RNA-seq biosynthesis of terpenoids, alcohols and esters characterize the transcriptomic profiles of cherry tomato fruit, harvested during postharvest storage under different temperatures and after 1-MCP treatment (Zou et al., 2018)
tomato multi-omics appearance, flavor, texture different tomato genotypes follow distinct transcriptomic, metabolomic (D’Esposito et al., 2017)
tomato multi-omics flavor, size, and production RNA-seq were capture genes associated structural variants (Alonge et al., 2020)
tomato multi-omics steroidal glycoalkaloids and pericarp metabolites genes with major effect on fruit metabolome (Zhu et al., 2018)
tomato multi-omics fruit flavor identified new genes and a rare allele for fruit flavor (Gao et al., 2019)
tomato multi-omics carotenoids identified DEGs related to ripening and carotenoid pathways (Zuo et al., 2020)
tomato multi-omics fruit quality (nitrate, soluble sugar, etc.) identification of DEGs in tomato fruits under substrate-Based and soil cultivation conditions (Guo et al., 2022)

Previous studies to understand quality traits in the solanaceae family using transcriptomics techniques.

3.1 Regulation of ascorbic acid metabolism

Ascorbic acid (AsA) is the most abundant antioxidant in plant cells. Pepper is among the vegetables with the highest AsA content (240 mg/100 g fresh weight) (Marín et al., 2004). In pepper, RNA-seq and targeted expression profiling have illuminated the molecular basis of nutritional quality and flavor compounds that define consumer preferences and commercial value. The expression levels of genes involved in AsA biosynthesis, recycling, and degradation have been analyzed in sweet pepper using quantitative real-time PCR, where an inverse correlation between AsA concentrations and the expression of the biosynthetic genes was observed (Alós et al., 2013).

3.2 Capsaicinoid biosynthesis and pungency determination

Capsaicinoids confer the characteristic pungency of chili peppers while providing protection against herbivores, insects, and microorganisms (Tewksbury and Nabhan, 2001). RNA-seq transcriptome profiling of pungent chili pepper ‘Xiaomila’ (C. frutescens) identified three structural genes encoding key enzymes in capsaicinoid biosynthesis: dihydroxy acid dehydratase (DHAD), threonine deaminase (TD), and prephenate aminotransferase (PAT) (Liu et al., 2013). Developmental time-course analysis in Serrano ‘Tampiqueño’ chili pepper revealed distinct temporal dynamics for these biosynthetic genes: DHAD expression remained stable during early fruit development (10–40 days post-anthesis, DPA) before declining at maturity (40–60 DPA), while TD expression decreased progressively throughout all developmental stages (Martínez-López et al., 2014). Comparative transcriptomics between pungent and non-pungent pepper placental tissue across development, alongside orthologous gene expression in tomato, identified the evolutionary innovations underlying capsaicinoid biosynthesis (Kim et al., 2014). This study revealed that activation of branched-chain amino acid aminotransferase ketoacyl-ACP synthase, and the acyltransferasegene during evolution within Solanaceae enabled the specialized biosynthetic capacity of chili peppers. These genes showed elevated expression specifically in pungent fruit placental tissue, confirming the tissue-specific nature of capsaicinoid biosynthesis and providing candidate markers for pungency breeding programs.

3.3 Ethylene-regulated ripening and flavor compound biosynthesis

Tomato has served as the archetypal model for understanding climacteric fruit ripening, where ethylene orchestrates dramatic metabolic transitions affecting color, flavor, and texture (Cara and Giovannoni, 2008; Mathieu et al., 2009). RNA-seq analysis of fruits treated with 1-methylcyclopropene (1-MCP), an ethylene inhibitor that extends shelf-life, identified 5,683 significantly differentially expressed genes (DEGs), revealing the extensive transcriptional network under ethylene control (Bobokalonov et al., 2018). 1-MCP treatment downregulated genes governing ethylene synthesis and signal transduction alongside carotenoid biosynthetic genes, linking ethylene perception to pigment accumulation that defines tomato fruit color. Integration of RNA-seq with DNA methylation profiling further demonstrated that ripening-associated genes, including ethylene biosynthesis enzymes (ACO2), signaling components (CTR1), ethylene response factors (ERF2, ERF5, ERF17, ERF114), and carotenoid pathway genes (PSY1, NSY, CrtR-b2, LCY1, NCED), are regulated through coordinated transcriptional and epigenetic mechanisms (Zuo et al., 2020).Storage temperature interactions with ethylene signaling profoundly affect flavor volatile production. Comparative transcriptomics across storage temperatures (25°C, 10°C, and 4°C) revealed that cold storage (4°C) suppressed expression of five genes responsible for biosynthesis of terpenoids, alcohols, and ester; key contributors to tomato flavor compared to moderate refrigeration (10°C) (Zou et al., 2018). These findings provide molecular explanations for the well-known consumer complaint that refrigerated tomatoes lack flavor, demonstrating that transcriptional suppression of volatile biosynthesis rather than metabolite degradation underlies this quality loss.

3.4 Phenylalanine and fatty acid-derived flavor networks

Beyond ethylene-regulated processes, tomato flavor depends heavily on phenylpropanoid and fatty acid-derived volatiles. Phenylalanine serves as a precursor for numerous flavor compounds, including methyl salicylate, a key phenylpropanoid contributor to tomato aroma (Tieman et al., 2010; Dal Cin et al., 2011; Kaur et al., 2023). Paradoxically, RNA-seq analysis revealed that organic nutrient supplementation increased phenylalanine accumulation but decreased phenylalanine-derived volatiles, including methyl salicylate, due to altered expression of LeAADC, LePAR, and SlSAMT genes (Li et al., 2022b). This counterintuitive result demonstrates that precursor availability alone does not determine volatile production and transcriptional control of biosynthetic and regulatory genes governs metabolic flux, with implications for strategies aimed at improving flavor. Enzymatic browning in fresh-cut potato products, caused by polyphenol oxidase activity on phenolic substrates, represents another critical quality defect (Cabezas-Serrano et al., 2009). Comparative transcriptomics between browning-resistant (‘Kexin 13’) and browning-susceptible (‘Yunshu 505’) cultivars revealed that DEGs were enriched in lipid metabolism and amino acid metabolism pathways, with phenolic substrates produced through phenylpropanoid biosynthesis closely associated with browning development (Qiao et al., 2022; Darshan et al., 2020). These findings provide both candidate genes for marker development and physiological targets for breeding low-browning potato varieties suitable for fresh-cut processing markets.

Fatty acid-derived volatiles, encompassing both long-chain and short-chain compounds, constitute another major flavor class. Comparative transcriptomics between processing tomato M82 and an introgression line identified SlLIP1, encoding a lipase responsible for lipid degradation to fatty acids; SlLIP1 knockdown lines showed reduced long-chain fatty acid-derived volatiles (Garbowicz et al., 2018). Similarly, RNA-seq combined with volatile profiling identified SlLIP8 as governing short-chain fatty acid-derived volatile production; expression levels correlated strongly with Z-3-hexen-1-ol and hexyl alcohol content, and SlLIP8 mutants exhibited substantially reduced volatile accumulation (Li et al., 2020). These lipase genes represent promising targets for metabolic engineering or marker-assisted selection to enhance specific volatile profiles, enabling breeding for improved flavor intensity or novel aroma characteristics.

3.5 Anthocyanin and chlorogenic acid regulation in eggplant

In eggplant, anthocyanins and chlorogenic acid (CGA) are the predominant secondary metabolites, contributing to fruit pigmentation, antioxidant capacity, and nutritional value (Jiang et al., 2016). The MYB transcription factor SmMyb1 emerged as a master regulator controlling both anthocyanin and CGA biosynthesis across eggplant tissues (Docimo et al., 2016). Quantitative PCR profiling revealed that SmMyb1 expression was relatively low across most tissues except stems, while functional analysis through transient overexpression in N. benthamiana demonstrated that C-terminal domain deletions specifically disrupted anthocyanin biosynthesis and CGA accumulation, indicating this region confers regulatory specificity. Furthermore, time-course RNA-seq analysis following light exposure (0, 0.5, 4, and 8 hours after bag removal) identified 1,956 DEGs, including 24 structural genes in anthocyanin biosynthesis pathways and 102 transcription factors exhibiting dynamic light-responsive expression changes (Li et al., 2018a). Notably, three photoreceptors: UV Resistance Locus 8 (UVR8), Cryptochrome 3 (CRY3), and (UVR3) were differentially expressed, along with light signal transduction elements COP1 and two SPA proteins, revealing the molecular cascade linking light perception to anthocyanin biosynthesis. This regulatory network provides targets for manipulating pigmentation patterns in eggplant breeding, with potential applications for developing cultivars with enhanced anthocyanin content for premium markets or specialized colors for consumer appeal.

Across Solanaceae crops, transcriptomic studies reveal several unifying principles governing secondary metabolite accumulation and quality trait expression. First, major regulatory hubs, including ethylene signaling in tomato, MYB transcription factors in eggplant, and temperature-responsive networks in potato coordinate expression of multiple biosynthetic pathways simultaneously, enabling pleiotropic effects from single genetic interventions. Second, environmental factors (temperature, light, nutrient availability) exert profound transcriptional effects that can override or modulate developmental programs, necessitating genotype-by-environment considerations in breeding. Third, counterintuitive relationships between precursor availability and end-product accumulation (exemplified by phenylalanine and AsA studies) highlight that transcriptional regulation of metabolic flux, rather than substrate concentration, often limits metabolite production. Finally, tissue-specific and temporal expression patterns indicate that effective breeding strategies must consider not just which genes to manipulate, but when and where they are expressed during development and post-harvest handling. These insights collectively enable more sophisticated breeding approaches that leverage expression-based markers, target regulatory genes for greater pleiotropic benefits, and account for environmental effects on quality trait expression.

4 Transcriptomics of solanaceous crops for biotic stress responses

Biotic stresses, including fungi, bacteria, viruses, nematodes, and insect pests represent the primary constraints limiting solanaceous crop productivity, with plant diseases alone causing substantial global yield losses (Kaur et al., 2019). While disease and pest management strategies include pathogen-free seed, chemical control, and cultural practices, deployment of resistant cultivars remains the most economically effective and environmentally sustainable approach (Mahmoud, 2021). However, resistance breeding requires precise identification and introgression of resistance genes (R-genes) alongside understanding of defense activation mechanisms (Lucas, 2011; Li et al., 2022a). Transcriptomic profiling, initially through microarrays and increasingly through RNA-seq has revolutionized this process by enabling genome-wide characterization of defense responses, identification of novel R-genes, and elucidation of regulatory networks governing pathogen recognition and downstream immunity.

4.1 Fungal pathogen defense mechanisms

Fungal pathogens constitute the predominant biotic threat to solanaceous crops, with late blight disease caused by Phytophthora infestans representing the most devastating example, inflicting up to 70% yield losses in both tomato and potato (Singh et al., 2021a). Comparative transcriptome analysis between P. infestans-infected resistant and susceptible tomato cultivars identified 1,037 DEGs and 688 long non-coding RNAs (lncRNAs), revealing previously unrecognized regulatory layers in disease resistance (Cui et al., 2017). Notably, lncRNA16397 was identified as an antisense transcript regulating expression of SlGRX22, a glutaredoxin gene involved in redox homeostasis; overexpression of both lncRNA16397 and SlGRX22 enhanced P. infestans resistance, demonstrating that non-coding RNAs represent untapped genetic resources for resistance breeding beyond traditional protein-coding R-genes.

The tomato-Cladosporium fulvum pathosystem, a classical model for studying gene-for-gene resistance interactions, exemplifies how transcriptomics elucidates R-gene function. RNA-seq analysis of C. fulvum-infected tomato cultivars carrying the Cf-19 resistance gene identified 418 DEGs mediating resistance responses, providing a comprehensive catalog of downstream genes activated by R-gene recognition (Zhao et al., 2019). This systems-level view reveals that single R-genes trigger coordinated expression changes across multiple defense pathways, informing expectations for resistance durability and potential trade-offs with growth or yield. In potato, transcriptomic studies have dissected defense responses to multiple necrotrophic and hemi biotrophic fungal pathogens. RNA-seq profiling of Alternaria solani (early blight) and Rhizoctonia solani (black scurf) infections revealed upregulation of biotic stress tolerance genes and pathogen defense transcripts, illuminating the distinct transcriptional signatures associated with necrotrophic pathogen lifestyles (Zrenner et al., 2020; Brouwer et al., 2021). Time-course transcriptomics of potato-Fusarium sulphureum (dry rot) interactions identified 1,025, 1,334, and 1,394 upregulated DEGs at 1-, 3-, and 7-days post-inoculation (dpi), respectively, with lignin biosynthesis genes, including PAL, 4CL, CAD, and POD showing progressively elevated expression (Fan et al., 2021). This temporal profiling demonstrates that structural defense responses intensify over time, suggesting that early resistance assessment may underestimate cultivar differences and that breeding programs should evaluate defense responses across multiple infection stages.

Eggplant transcriptomics has leveraged wild relatives to identify resistance mechanisms absent in cultivated germplasm. While Solanum melongena is highly susceptible to the soil-borne fungus Verticillium dahliae (Verticillium wilt), the wild relative S. aculeatissimum exhibits robust resistance. Transcriptome analysis of V. dahliae-infected S. aculeatissimum roots revealed 11,696 upregulated and 5,949 downregulated genes, with enrichment in R-genes, pathogenesis-related (PR) proteins, and phenylpropanoid pathway genes (Zhou et al., 2016). This massive transcriptional reprogramming encompassing nearly 30% of the expressed genome indicates that resistance involves coordinated activation of multiple defense layers rather than single-gene effects, providing numerous targets for introgression into cultivated eggplant through marker-assisted breeding or genomic selection approaches.

4.2 Nematode resistance networks

Nematodes, particularly root-knot nematodes (RKNs) and cyst nematodes, inflict severe damage on vegetable crops through root system colonization. In potato, comparative transcriptomics between potato cyst nematode (PCN)-resistant (‘Kufri Swarna’) and susceptible (‘Kufri Jyoti’) cultivars identified 438 upregulated and 353 downregulated genes in the resistant line, with upregulated genes enriched for disease resistance genes and transcription factors including WRKY, HMG, and MYB families (Chandrasekar et al., 2022). These transcription factors likely serve as master regulators for coordinating expression of downstream defense genes, representing high-value targets for breeding interventions that could confer broad-spectrum resistance. Tobacco transcriptome profiling following Meloidogyne incognita (RKN) invasion revealed 4,354 differentially expressed genes in resistant genotypes, encompassing genes involved in cell wall modification, oxidative burst, signal transduction, and transcriptional regulation by ERF, MYB, and WRKY factors (Li et al., 2018c). The prominence of cell wall modification genes suggests that physical barriers preventing nematode penetration or migration constitute a key resistance mechanism. Similarly, transcriptome analysis of S. torvum, a wild eggplant relative used as nematode-resistant rootstock identified 5,360 DEGs following RKN infection, enriched in protein phosphorylation, defense signal transduction, and plant-pathogen interaction pathways (Zhang et al., 2021b). The convergence on signaling and phosphorylation pathways across multiple nematode-host systems indicates that rapid signal perception and amplification are critical for effective nematode resistance, suggesting that breeding for enhanced early-response signaling components could improve resistance across diverse nematode species.

4.3 Viral disease resistance mechanisms

Unlike fungal and bacterial diseases, viral infections in solanaceous crops lack effective chemical control options, rendering genetic resistance the sole viable management strategy. In pepper, chili leaf curl virus (ChiLCV) threatens production across South Asia. Suppression subtractive hybridization (SSH) followed by transcriptomic profiling of ChiLCV-infected resistant cultivar ‘Punjab Lal’ identified 231 unique ESTs related to cellular and physiological processes governing resistance (Kushwaha et al., 2015). Subsequent microarray analysis comparing susceptible and resistant Capsicum species identified 319 DEGs, with 234 unique genes upregulated more than 2-fold in resistant lines, revealing a broader resistance network than initially captured by SSH (Rai et al., 2016). The progression from SSH to microarray demonstrates how technological advances enable increasingly comprehensive characterization of resistance mechanisms.

In tomato, comprehensive transcriptomic analysis of tomato spotted wilt virus (TSWV)-infected plants carrying the Sw-7 resistance gene identified 1,244 DEGs involved in callose deposition, lignin accumulation, RNA silencing, and transcriptional regulation (Padmanabhan et al., 2019). Functional characterization revealed that PR-5 (pathogenesis-related protein 5) directly participates in Sw-7-mediated resistance, providing a specific molecular target for diagnostic marker development or engineering enhanced TSWV resistance. This study exemplifies how transcriptomics bridges the gap between R-gene identification and mechanistic understanding of resistance execution.

Tobacco transcriptomics has illuminated complex virus-betasatellite interactions. RNA-seq comparison between tobacco plants infected with Tobacco curly shoot virus (TbCSV) alone versus co-infection with Tobacco curly shoot betasatellite (TbCSB) revealed 3,196 and 4,081 DEGs, respectively, with 7 and 13 genes showing remarkable effects on brassinosteroid (BR) and jasmonic acid (JA) biosynthesis and signaling (Li et al., 2018b). The symptom shift from upward to downward leaf curling with betasatellite co-infection correlates with distinct transcriptional reprogramming of hormone pathways, demonstrating that virus-satellite interactions reprogram host physiology beyond the effects of either agent alone. These findings have implications for understanding disease complexes in solanaceous crops where multiple pathogens co-occur.

4.4 Insect pest defense and interspecific variation

Comparative transcriptomics between tomato and eggplant infested by leaf miner (Tuta absoluta) revealed higher expression levels of plant immunity-related transcripts in eggplant, suggesting differential constitutive or induced defense capacity between these solanaceous species (Chen et al., 2021). This interspecific variation indicates that genes underlying eggplant’s enhanced immunity could be introgressed into tomato through wide crosses or identified as targets for engineering improved insect resistance. Understanding species-specific defense architectures enables strategic exploitation of genetic diversity within Solanaceae for pest resistance breeding.

Despite pathogen diversity spanning fungi, nematodes, viruses, and insects, transcriptomic studies across solanaceous crops reveal convergent defense mechanisms and regulatory architectures. First, transcription factor families particularly WRKY, MYB, and ERF emerge repeatedly as master regulators coordinating pathogen-specific responses, representing high-leverage targets for breeding interventions that could provide broad-spectrum resistance. Second, phenylpropanoid pathway activation constitutes a universal defense response across pathogen types, producing antimicrobial compounds, lignin for physical barriers, and signaling molecules for systemic resistance. Third, pathogenesis-related (PR) proteins, including PR-5, chitinases, and β-1,3-glucanases, are consistently upregulated across resistance responses, suggesting that their expression levels could serve as quantitative predictors of resistance capacity. Furthermore, the discovery of lncRNAs regulating defense genes (lncRNA16397/SlGRX22) reveals previously unrecognized regulatory layers, expanding the pool of genetic variation accessible for resistance improvement beyond traditional R-gene pyramiding.

The massive transcriptional reprogramming observed in resistant responses up to 30% of the transcriptome in V. dahliae-resistant eggplant suggests that effective resistance often involves coordinated multi-gene networks rather than single-gene effects. This systems-level perspective supports genomic selection approaches that capture small-effect loci collectively contributing to quantitative resistance, complementing traditional R-gene deployment. Finally, wild relatives (S. aculeatissimum, S. torvum) harbor resistance mechanisms absent in cultivated germplasm (Zhang et al., 2021b), and transcriptomics provides the molecular roadmap for identifying and introgressing these valuable alleles into elite breeding lines, accelerating the development of durable, broad-spectrum resistance across solanaceous crops (Table 2).

Table 2

Crop Technique Biotic stress (organism) Trait References
tomato RNA-seq Tomato yellow leaf curl virus (TYLCV) DEGs involved in cell wall reorganization, transcriptional regulation, defense response, ubiquitination, metabolite synthesis (Chen et al., 2013)
tomato microarray Botrytis cinerea profile the expression of genes for the biosynthesis, modification and signal transduction (Blanco-Ulate et al., 2013)
tomato RNA-seq C. globosum deciphered the defense signaling pathways associated with systemic resistance induced by Chaetomium globosum (Singh et al., 2022)
tomato RNA-seq tomato spotted wilt virus (TSWV), cucumber mosaic virus (CMV), and potato virus Y (PVY) grafting alters tomato transcriptome that increases tolerance to an airborne virus infection (Spanò et al., 2020)
tomato RNA-seq Verticillium dahliae molecular mechanisms of tomato plant susceptibility (Tan et al., 2015)
tomato microarray Ralstonia. solanacearum elucidate the molecular mechanisms in the early resistance response (Ishihara et al., 2012)
tomato microarray Ralstoniasolanacearum transcript level of biochemical and immunohistochemical work on silicon induced resistance (Ghareeb et al., 2011)
tomato RNA-seq TSWV DEGs were identified throughout a disease progression process involving networks of host resistance genes, RNA silencing/antiviral defense genes, and crucial transcriptional and translational regulators. (Padmanabhan et al., 2019)
tomato RNA-seq Phytophthora infestans identification of lncRNA16397 conferring resistance by co-expressing glutaredoxin (Cui et al., 2017)
tomato RNA-seq Phytophthora infestans SpWRKY3 was screened from transcriptome analysis of tomatoes inoculated with and without pathogens (Cui et al., 2018)
tomato RNA-seq Cladosporium fulvum analyse the differences between the response of resistant plants (carrying the Cf-19 gene) and susceptible plants (Moneymaker (MM), carrying the Cf-0 gene) at 0, 7 and 20 days after inoculation (Zhao et al., 2019)
tomato RNA-seq Clavibacter michiganensis subsp. michiganensis transcriptome response of resistant and susceptible tomato lines the infection (Basim et al., 2021)
tomato RNA-seq Meloidogyne incognita resistance against root-knot nematode (Shukla et al., 2018)
tomato RNA-seq TYLCV determine the distribution of the viral transcripts accumulated during the infection (Piedra-Aguilera et al., 2019)
tomato microarray/RNA-seq Potato spindle tuber viroid transcriptomic profile changes in tomato roots systemically infected with mild or severe PSTVd variants (Góra-Sochacka et al., 2019)
potato RNA-seq Potato Virus A (PVA) identify the genes involved in the response to PVA infection and potential candidate genes for improving resistance cultivars (Li et al., 2017)
potato RNA-seq Phytophthora infestans (Pi) exogenous ethylene activates immune and defense responses in potato against late blight (Yang et al., 2020a)
potato RNA-seq P. infestans genes associated with late blight resistance in potato (He et al., 2021)
potato microarray PVY transcript and protein abundance of potato following the infection with PVY (Stare et al., 2017)
potato RNA-seq Alternaria solani identified the differential expression of several potato and A. solani transcripts that present a group of valuable candidates into the roles in immunity or disease development. (Brouwer et al., 2021)
potato RNA-seq Ralstonia solani necrotrophic lifestyle of R. solani AG3-PT during early interaction with its host (Zrenner et al., 2020)
potato RNA-seq Fusarium sulphureum. involvement of lignin synthesis-related pathway and MAPK signaling pathway-plant in resistance (Fan et al., 2021)
potato RNA-seq Synchytrium endobioticum transcriptome analysis of resistance mechanism to potato wart disease (Li et al., 2021a)
potato RNA-seq Potato cyst nematodes identification of disease resistance genes and molecular markers for PCN infestation (Chandrasekar et al., 2022)
potato DeepSAGE P. infestans DeepSAGE transcriptome analysis uncovered novel candidate genes for plant host pathogen interactions (Gyetvai et al., 2012)
pepper RNA-seq Xanthomonas campestris pv. vesicatoria DEGs from transcriptomic analysis of a susceptible cultivar and a resistant cultivar were identified (Gao et al., 2021)
pepper RNA-seq Colletotrichum capsici resistance mechanism and identification of markers for anthracnose resistance in chili (Kumar et al., 2022)
pepper RNA-seq Chilli leaf curl virus(Chi LCV) protein homeostasis and defense in resistant chili plants against Chi LCV infection (Kushwaha et al., 2015)
eggplant RNA-seq V. dahliae transcriptome of the wild eggplant species S. aculeatissimum in response to V. dahliae (Zhou et al., 2016)
eggplant RNA-seq Tuta absoluta transcriptome analysis showed higher transcript changes in infested eggplant than tomato (Chen et al., 2021)
eggplant RNA-seq Meloidogyne incognita (root-knot nematode - RKN) transcriptomic profiles of eggplant and Solanum torvum were compared. Response-related DEGs were identified as being involved in stimulus response, protein phosphorylation, hormone signaling, and plant-pathogen interactions. (Zhang et al., 2021b)
tobacco RNA-seq Tobacco curly shoot virus (TbCSV) and Tobacco curly shoot betasatellite (TbCSB) transcriptome analysis N. benthamiana infected by TbCSV (Li et al., 2018b)
tobacco RNA-seq Cucumber mosaic virus (CMV) mechanisms of CMV tolerance in plants (Liu et al., 2019)
tobacco RNA-seq CMV analysis of N. tabacum infected by CMV during systemic symptom development (Lu et al., 2012)
tobacco RNA-seq Phytophthora nicotianae candidate resistant and susceptible genes against P. nicotianae infections (Meng et al., 2021)
tobacco RNA-seq M. incognita defense-related genes against M. incognita invasion in tobacco (Li et al., 2018c)

Description of several transcriptomics techniques used to study biotic stresses in solanaceous crops.

5 Transcriptomics for solanaceous crops for abiotic stress responses

Plants in natural settings are confronted by multiple abiotic stresses including drought, salinity, temperature extremes, and nutrient limitations that cause significant crop yield losses (Kaur et al., 2022). These stresses ultimately generate reactive oxygen species (ROS) that interfere with protein-protein interactions, disrupt ionic and osmotic homeostasis, and cause tissue and cell death (Singh et al., 2021b), severely hampering physiological growth, yield, and productivity. In response, plants have evolved adaptive mechanisms operating at molecular, cellular, and physiological levels, with transcriptional reprogramming representing the first line of defense (Gong et al., 2020). Central to these responses are transcription factors (TFs); trans-acting regulatory proteins that bind promoter regions to control expression of target genes governing stress tolerance. To date, sixty TF families have been identified in solanaceous crops, with WRKY, DOF, MYB, bZIP, ARF, ERF, HSF, and NAC families most prominently implicated in abiotic stress defense (Tolosa and Zhang, 2020). Transcriptome profiling across tissue types and developmental stages has enabled systematic identification of functional genes and regulatory networks underlying stress tolerance, providing the molecular foundation for breeding climate-resilient cultivars capable of maintaining productivity under increasingly variable environmental conditions.

5.1 Drought and salinity stress responses in tomato

Tomato, sensitive to heat, salinity, cold, and drought stress, has served as the primary model for dissecting abiotic stress transcriptional networks in Solanaceae. Comparative transcriptomics under drought and salinity stress identified the bZIP transcription factors SlAREB1 and SlAREB2 as master regulators, with expression induced in both leaves and root tissues (Orellana et al., 2010). Microarray and cDNA-amplified fragment length polymorphism (cDNA-AFLP) analysis of SlAREB1-overexpressing lines revealed production of oxidative stress-response proteins, defense-related proteins, protease inhibitors, and additional transcription factors, demonstrating that single TF manipulations can trigger cascading activation of entire stress tolerance networks. Key candidate genes associated with drought tolerance included AHG2 (ABA signaling), PRXIIF and PRXQ (peroxidases for ROS detoxification), SAP5 (stress-associated protein), and CCD1 (carotenoid cleavage dioxygenase), collectively spanning signaling perception, ROS management, and metabolic adjustment.

For salinity tolerance specifically, critical candidate genes mapped to HKT1 (high-affinity potassium transporter) on chromosome 7 and LeNHX3 (Na+/H+ exchanger) on chromosome 1, which regulate ion homeostasis by controlling sodium exclusion and compartmentalization. RNA-seq analysis of tomatoes growing under seawater salinity stress showed upregulation of DEGs responsible for growth maintenance, signal transduction, and protein-binding pathways (Usharani et al., 2022), indicating that salinity tolerance requires coordinated transcriptional activation spanning multiple cellular processes beyond ion transport alone. Functional analysis under combined drought and salinity stress revealed elevated levels of hydrogen peroxide (signaling molecule), melatonin (antioxidant), and zeatin riboside (cytokinin) during stress conditions (Zhou et al., 2019) (Table 3), demonstrating that plants deploy overlapping regulatory mechanisms for multiple concurrent stresses; a critical consideration for breeding cultivars adapted to real-world field conditions where multiple stresses co-occur.

Table 3

Crop Technique Abiotic stress Characteristics References
eggplant RNA-seq low temperature plant hormone signal transduction related genes were involved in eggplant response to low temperature (Yang et al., 2020b)
tomato RNA-seq drought In-depth analysis of DEGs revealed that ABA metabolism and polyamine biosynthesis elevated to provide tolerance (Nie et al., 2022)
potato RNA-seq heat genes involved in amino acid biosynthesis and secondary metabolism were mostly induced during heat exposure (Liu et al., 2021a)
potato RNA-seq heat DEGs were clustered into 49 different GO types, reflecting the functional diversity of the heat stress response genes (Tang et al., 2020)
potato RNA-seq cold the expression of 15 heat shock proteins and the accumulation of soluble sugars (mediated by enzymes such as invertase inhibitor and fructokinase) to prevent cellular damage (Lin et al., 2019)
potato RNA-seq drought differential gene expression trend reflected the upregulation of genes involved in carbohydrate metabolism, flavonoid biosynthesis, lipid biosynthesis/transfer, heat shock proteins and secondary metabolites like phenolics and lignins. (Aliche et al., 2021)
potato RNA-seq nitrogen deficiency DEGs like glutamate dehydrogenase, glutamine synthetase and carbonic anhydrasehelped to combat nitrogen stress (Zhang et al., 2020)
pepper RNA-seq heat, cold, salinity, and osmotic stress Analyses of the transcriptome data in this study will provide useful information for basic studies of various stimuli to facilitate the development of stress-resistant pepper cultivars (Kang et al., 2020)
pepper RNA-seq heat DEGs revealed the most enriched pathways were metabolic processes under stress and photosynthesis and light harvesting during stress (Wang et al., 2021)
pepper microarray salt stress signaling genes played important role in heat stress (López−Serrano et al., 2021)
pepper RNA-seq drought calcineurin B-like proteins and CBL-interacting protein kinases are important TFs during multiple stress conditions (Ma et al., 2021b)
tobacco RNA-seq chilling stress/low temperature DEGs were primarily enriched in the following functions: regulating cellular osmotic balance, promoting photosystem II repair, modulating signal transduction pathways such as ABA/GA, and enhancing lipid metabolism and lignin synthesis to improve membrane stability and cell wall mechanical strength (Zhou et al., 2020)
tobacco RNA-seq salt stress the up-regulated genes in salt stress were mainly concentrated in transcription factor WRKY family and PAR1 resistance gene family (Wu et al., 2021)
tobacco RNA-seq drought upregulation of genes involved in arginine, proline, and starch metabolism (Khan et al., 2019)
tobacco RNA-seq salinity acetylcholine plays an important role in downregulating ROS and hydrogen peroxide level (Qin et al., 2021)
tobacco microarray environmental changes DEGs among three cultivated regions were associated with the light reaction of photosystem II, response to stimuli, and secondary metabolism. (Lei et al., 2013)

Description of different TFs and regulation of genes during abiotic stress conditions in plants.

5.2 Heat shock response and temperature stress adaptation

Heat stress transcriptional responses have been extensively characterized across multiple solanaceous crops, revealing conserved heat shock factor (HSF) regulatory networks. In tomato, RNA-seq analysis showed that spermidine, a polyamine plays an important regulatory role in heat stress tolerance, while endoplasmic reticulum (ER) stress responses constitute a critical factor contributing to varietal heat tolerance differences (Usharani et al., 2022). The connection between ER stress and heat tolerance suggests that protein folding capacity and quality control mechanisms limit thermotolerance, providing targets for breeding interventions.

In potato, candidate genes encoding heat shock transcription factors including StHsf004, StHsf007, StHsf009, StHsf014, and StHsf019 showed constitutive expression under non-stress conditions, providing basal thermotolerance capacity, while additional HSFs were induced specifically under heat stress (Tang et al., 2016). Comprehensive transcriptome analysis of heat-sensitive potato variety ‘Agria’ identified 2,949 differentially expressed genes responding to heat stress, resulting in accumulation of heat shock proteins that serve as molecular chaperones preventing protein denaturation and aggregation (Usharani et al., 2022). In pepper, RNA-seq comparison between heat-tolerant and heat-sensitive varieties revealed that the glutathione metabolic pathway plays a critical role in heat tolerance (Wang et al., 2020), providing antioxidant capacity to neutralize heat-induced ROS production. These findings collectively indicate that heat tolerance involves both protective protein expression (chaperones) and ROS management systems, suggesting that breeding strategies should target both mechanisms simultaneously.

Combined drought and heat stress transcriptomics in potato revealed differential regulation of at least 19,000 transcripts, including genes encoding type-1 metacaspase (DMT400032693, programmed cell death), ABI5-binding protein 2 (DMT400036544, ABA signaling), cell wall invertase (DMT400023092, sugar metabolism), tryptophan synthase (DMT400029363, amino acid biosynthesis), and lipid biosynthesis enzymes (DMT400033492) (Demirel et al., 2020). Critically, tolerant varieties encoded fewer transcripts involved in programmed cell death pathways, suggesting that stress tolerance involves not just activating protective mechanisms but also suppressing cell death signaling, a balance that determines whether plants survive or succumb to stress combinations.

5.3 Cold stress and low-temperature acclimation

Cold stress responses involve distinct transcriptional programs from heat stress, with phenylpropanoid pathway activation emerging as a common defense strategy. In potato, functional analysis identified microRNAs stu-novel-miR28211 and stu-novel-miR43095 as regulators of cold tolerance, modulating peroxidase levels in the phenylpropanoid pathway (Liao et al., 2021). The phenylpropanoid pathway produces defensive metabolites, lignin for cellular reinforcement, and antioxidant compounds representing a multi-layered cold protection strategy.

In eggplant, cold stress induced expression of WRKY transcription factors SmWRKY26 and SmWRKY32 (Yang et al., 2020b), while paradoxically downregulating anthocyanin biosynthesis genes, including CHI, 3GT, F3′5′H, DFR2, ANS, and associated transcription factors resulting in lower anthocyanin content at high temperatures (Lv et al., 2019). This temperature-dependent metabolic shift suggests that anthocyanin biosynthesis is optimized for moderate temperatures, with resources redirected toward other protective mechanisms under temperature extremes.

Potato tuber quality during storage presents distinct challenges related to cold-induced sweetening and enzymatic browning. Combined transcriptomic and proteomic analysis of tubers stored at 15°C, 4°C, and 0°C revealed that cold stress triggered soluble sugar accumulation through coordinated regulation of granule-bound starch synthase 1, beta-amylase, invertase inhibitor, and fructokinase (Lin et al., 2019). This cold-induced sweetening reduces processing quality by promoting undesirable browning during frying, a major economic concern for the potato industry. Understanding the transcriptional basis of sugar accumulation enables breeding for cold-sweetening resistant cultivars or identification of storage conditions that minimize quality deterioration.

In tobacco, transcriptome analysis following chilling stress revealed upregulation of 1,170 genes and downregulation of 505 genes at 12- and 24-hours post-treatment (Zhou et al., 2020). These genes regulate cell osmotic potential and soluble sugar levels while promoting PSII (photosystem II) repair, lipid metabolism adjustments, and lignin synthesis to enhance cell wall mechanical strength. The emphasis on PSII repair highlights that photosynthetic apparatus damage represents a primary mechanism of cold injury, and tolerance requires active repair systems rather than passive protection alone.

5.4 Salinity tolerance mechanisms across species

Beyond tomato, salinity transcriptomics in other solanaceous crops has revealed both conserved and species-specific tolerance mechanisms. In salt-tolerant eggplant variety ‘ST118’, DEGs included transcription factors from C2C2-CO-like, WRKY, and MYB families, alongside ion transport genes AKT1 and KAT1 (potassium channels) and SOS1 (salt overly sensitive 1, Na+/H+ antiporter) (Li et al., 2018a). The co-regulation of transcription factors and ion transporters indicates that effective salinity tolerance requires both sensing/signaling components and effector mechanisms for ion homeostasis.

In tobacco, RNA-seq analysis identified 1,654 DEGs involved in salt stress responses, with upregulation of WRKY family and PAR1 family transcription factors providing salinity tolerance (Wu et al., 2021). This study also revealed the critical role of abscisic acid (ABA) in mediating salt resistance, with ABA-responsive genes constituting a major component of the salinity tolerance network. Transgenic tobacco seedlings overexpressing ArSQE (squalene epoxidase) along with bZIP transcription factors showed tolerance to multiple abiotic stresses (Upadhyay et al., 2020), demonstrating that manipulating key metabolic enzymes alongside regulatory factors can confer broad stress tolerance, a promising strategy for developing multi-stress resilient cultivars.

5.5 Drought tolerance networks and signaling cascade

Despite diversity in stress types and solanaceous species, transcriptomic studies reveal several convergent molecular principles governing abiotic stress tolerance. First, a core set of transcription factor families particularly WRKY, MYB, bZIP, ERF, and NAC function as master regulators across all stress types and species, indicating that these TF families possess inherent stress-responsive regulatory capacity. Manipulating expression or DNA-binding specificity of these master regulators represents a high-leverage breeding strategy potentially conferring tolerance to multiple stresses simultaneously.

Second, ROS management through antioxidant systems (glutathione metabolism, peroxidases, carotenoid-derived signals) constitutes a universal stress response mechanism, as ROS generation is the common downstream consequence of most abiotic stresses. Breeding for enhanced antioxidant capacity through selection on glutathione pathway genes or peroxidase expression levels could provide baseline stress tolerance across diverse environmental challenges.

Third, programmed cell death pathways emerge as critical determinants of stress outcomes: tolerant genotypes suppress cell death signaling while maintaining protective responses, whereas susceptible genotypes undergo premature tissue death. This suggests that stress tolerance breeding should focus not only on activating protective mechanisms but also on dampening cell death pathways, a dual strategy rarely considered in conventional breeding programs.

Fourth, the massive transcriptional reprogramming under combined stresses (19,000 transcripts in potato drought+heat stress) indicates that tolerance to stress combinations cannot be predicted from single-stress responses. This necessitates phenotyping and transcriptomic profiling under realistic multi-stress conditions that mirror field environments, rather than controlled single-stress treatments that may poorly predict field performance.

Fifth, tissue-specific and temporal dynamics of stress responses (e.g., SlAREB expression in both leaves and roots; time-course changes in tobacco chilling response) highlight that effective stress tolerance requires coordinated responses across organs and sustained transcriptional commitment over time. Breeding programs should therefore evaluate stress tolerance across multiple tissues and time points to capture genotypes with robust, sustained responses rather than transient stress reactions.

Finally, the discovery of regulatory microRNAs (stu-novel-miR28211, stu-novel-miR43095) and hormone signaling components (ABA, melatonin, cytokinins) reveals additional regulatory layers beyond protein-coding genes. Expanding breeding targets to include regulatory RNAs and hormone pathway components could unlock genetic variation inaccessible through traditional gene-focused approaches. These transcriptomic insights collectively provide a molecular roadmap for developing climate-resilient solanaceous cultivars capable of maintaining productivity under the increasingly variable and extreme environmental conditions projected for agricultural regions worldwide (Table 3).

6 Transcriptomic approaches to dissecting yield architecture and productivity

Yield represents the ultimate breeding target, yet its highly polygenic nature and sensitivity to environmental conditions have historically made it challenging to dissect genetically. Transcriptomics provides a powerful approach for understanding yield determination by capturing gene expression dynamics across tissues, developmental stages, and environmental conditions revealing regulatory networks that translate genetic variation into phenotypic outcomes (Lowe et al., 2017). Unlike single-timepoint genomic analyses, transcriptome profiling captures the active molecular processes governing resource allocation, organ development, and stress integration that collectively determine productivity. In solanaceous crops, yield improvement requires simultaneous optimization of fruit or tuber morphology, nutrient use efficiency, photosynthetic capacity, and reproductive development—all traits amenable to transcriptomic dissection. Moreover, as emphasized in the introduction, achieving genetic gain under climate change necessitates understanding how yield-determining processes respond to heat, drought, and other stresses. This section examines how transcriptomics has illuminated the molecular basis of yield components, providing actionable targets for breeding higher-yielding, climate-resilient cultivars.

6.1 Fruit and tuber morphology: developmental programs determining harvestable yield

Fruit size and shape constitute primary yield determinants in solanaceous crops, with morphological variation arising from early developmental programs that establish organ identity and growth trajectories. In tomato, genetic and transcriptomic analyses identified ENO (enolase) as a regulator of fruit size, with RNA-seq revealing that a promoter mutation selected during domestication increased ENO expression in cultivated tomatoes, driving larger fruit development (Yuste-Lisbona et al., 2020). This finding demonstrates that cis-regulatory variation, rather than protein-coding changes, can generate major morphological shifts, suggesting that breeding programs should screen regulatory regions, not just coding sequences, to unlock yield-enhancing alleles.

In potato, tuber shape determination involves complex developmental transitions from stolon to tuber. Combined genetic and transcriptomic analysis using RNA-seq across tissue sections identified five candidate genes governing tuber morphology: Soltu.DM.10G018510 (nuclear shuttle interacting protein), Soltu.DM.10G018550 (aminophospholipid ATPase, membrane trafficking), Soltu.DM.10G018580 (seven-transmembrane receptor), Soltu.DM.10G018620 (HSI2-like transcriptional repressor), and Soltu.DM.10G018660 (lipid transfer protein) (Fan et al., 2022). Critically, transcriptome sequencing revealed that developmental divergence between round and elongated tuber morphologies initiates as early as the stolon hook stage, much earlier than previously recognized. This temporal resolution demonstrates that yield-determining developmental decisions occur during vegetative rather than reproductive phases, indicating that breeding for tuber shape requires screening at early developmental stages and that selection based on mature tuber phenotypes may miss critical genetic variation affecting early developmental commitment.

In eggplant, comparative transcriptomics between four wild accessions (Solanum insanum) and 16 domesticated lines identified genomic regions targeted during domestication that could enhance yield improvement (Page and Chapman, 2021). Analysis of fixation index (FST) and nucleotide diversity (π) revealed low diversity in domesticated populations, indicating selection sweeps in yield-related regions. These domestication-targeted regions likely harbor alleles balancing fruit size against plant vigor or stress tolerance, trade-offs that could be re-optimized through genomics-informed breeding that reintroduces favorable wild alleles into elite backgrounds.

6.2 Integrating morphology with stress tolerance: multi-trait transcriptional networks

An innovative study demonstrated that SlMX1, encoding a MIXTA-like MYB transcription factor that regulates trichome formation, simultaneously influences yield, stress tolerance, and secondary metabolism (Ewas et al., 2022). Comparative transcriptomics between SlMX1-overexpressing transgenic tomatoes and wild-type plants revealed that SlMX1 upregulated and downregulated genes governing cell growth, primary metabolism, and secondary metabolite biosynthesis, ultimately improving fruit yield while enhancing resistance to both biotic and abiotic stresses. This pleiotropic TF represents an ideal breeding target for addressing the consumer-producer trade-off. The identification of such master regulators through transcriptomics enables breeding strategies that avoid negative genetic correlations by selecting regulatory variants affecting multiple beneficial traits concurrently.

Similarly, an unexpected finding demonstrated that RNA N6-methyladenosine (m6A) modification, an epitranscriptomic layer regulating mRNA stability and translation increases potato yield and biomass in field trials (Yu et al., 2021). This reveals that yield improvement can be achieved not just through altering which genes are transcribed, but through modulating post-transcriptional RNA processing, expanding the pool of genetic variation accessible for breeding beyond traditional expression quantitative trait loci (eQTL).

6.3 Nutrient use efficiency: metabolic networks linking resource acquisition to productivity

Nitrogen use efficiency (NUE), the capacity to maintain productivity under limiting nitrogen availability represents a critical breeding target for sustainable agriculture and yield stability under variable nutrient conditions. In potato, RNA-seq analysis identified superoxide dismutase (ROS detoxification), probable protein phosphatase 2C (signaling), and high-affinity nitrate transporters as major candidate genes governing NUE. The involvement of ROS detoxification enzymes indicates that nitrogen limitation generates oxidative stress, and that NUE requires coordinated upregulation of both nitrogen acquisition and antioxidant systems, a systems-level understanding unattainable from genomic analysis alone.

In tobacco, RNA-seq following nitrogen deficiency identified 428 upregulated and 213 downregulated genes with distinct pathway enrichment patterns (Li et al., 2021b). Upregulated genes were enriched in MAPK signaling (stress perception), sesquiterpenoid and triterpenoid biosynthesis (secondary metabolism), and arginine and proline metabolism (compatible solutes), while downregulated genes were enriched in photosynthesis, nitrogen metabolism, and amino acid biosynthesis. This transcriptional reprogramming reveals an apparent paradox: nitrogen-deficient plants downregulate nitrogen assimilation and amino acid biosynthesis, the very processes needed to utilize available nitrogen. This suggests that severe nitrogen limitation triggers a resource conservation strategy prioritizing survival (stress signaling, compatible solutes) over growth, with implications for fertilization timing: applications must precede severe deficiency to maintain productive metabolism.

In eggplant, comparative transcriptomics between nitrogen-use efficient and inefficient genotypes revealed that genes involved in light reactions, particularly light-harvesting complexes (LHCs) and ferredoxin-NADP reductase (FNR) were upregulated in efficient genotypes (Mauceri et al., 2021). This counterintuitive finding indicates that NUE depends critically on photosynthetic efficiency: efficient genotypes maintain higher photosynthetic capacity under nitrogen limitation, generating sufficient carbon skeletons and ATP to maximize utilization of available nitrogen. This demonstrates that breeding for NUE requires selecting for maintained photosynthetic function under nutrient stress, not just enhanced nitrogen uptake, a physiological target that could be achieved through marker-assisted selection on LHC and FNR expression levels.

6.4 Inflorescence architecture and reproductive development

Inflorescence architecture directly determines potential fruit number per plant, representing a major yield component. Meristem transcriptomic profiling in tomato revealed that allelic dosage at specific loci controlling inflorescence branching improves architecture and yield (Soyk et al., 2017), demonstrating that yield increases can be achieved through optimizing allelic combinations rather than identifying novel genes, a principle applicable to genomic selection strategies that capture small-effect quantitative trait loci collectively contributing to complex traits.

Stigma position, which affects pollination efficiency and fertilization success, represents another critical reproductive trait. Comparative transcriptomics between exserted-stigma (improved pollination) and inserted-stigma genotypes identified 801 DEGs—566 upregulated, 235 downregulated (Riccini et al., 2021). Upregulated genes in exserted-stigma lines included those governing cell wall structure and metabolism, particularly sucrose metabolism genes (std1, sps, SuSy), the auxin efflux facilitator PIN-FORMED3, and two xyloglucan endotransglucosylase-hydrolases (XTH1, XTH7) mediating cell wall expansion. Hormone-related genes were prominently affected, including three Auxin Response Factors (ARFs), two Auxin/Indole-3-Acetic Acid proteins (Aux/IAA), two Arabidopsis Response Regulators (ARRs), and two GIBBERELLIN-STIMULATED TRANSCRIPT 1 (GAST1)-like genes. The convergence on auxin and gibberellin pathways indicates that stigma elongation results from coordinated hormone signaling driving localized cell expansion, providing candidate genes for breeding improved pollination efficiency under heat stress conditions that disrupt normal stigma development.

Parthenocarpy, fruit development without fertilization offers yield stability under pollination-limiting conditions (heat stress, pollinator scarcity). Microarray analysis of a parthenocarpic tomato mutant (pat) with pleiotropic floral traits identified novel genes regulating fruit set (Ruiu et al., 2015). Differentially expressed genes segregated into regulatory function (RF) and pollination-dependent (PD) groups. RF genes including LeT6/TKn2 (KNOTTED-like homeobox), SlBPEp, SlCRCb, and SlDEF showed downregulation in the pat mutant, while PD genes enriched in auxin perception and response including ARF and Aux/IAA family members (SlARF9, SlIAA14) exhibited altered expression patterns. This regulatory network dissection provides molecular handles for engineering facultative parthenocarpy that maintains normal fruit development under optimal conditions while enabling fruit set under stress, a breeding objective directly addressing climate resilience.

During fruit ripening, oligonucleotide microarray analysis revealed that most plastid genes are transcriptionally downregulated (Kahlau and Bock, 2008), with transcriptional, post-transcriptional, and translational controls all contributing to developmental regulation of plastid gene expression. This multilayered regulation of photosynthetic capacity during ripening affects carbon partitioning between growth and fruit filling, with implications for breeding optimized source-sink relationships.

6.5 Photosynthetic optimization and light response networks

Photosynthetic efficiency directly determines carbon fixation capacity and ultimately yield potential. Leaf morphology affects photosynthetic efficiency through altered light interception and physiological parameters, impacting fruit sugar accumulation and yield (Nicotra et al., 2011). RNA-seq analysis investigating associations between photosynthetic gene expression and leaf morphology in tomato identified transcriptional programs linking leaf shape variation to photosynthetic capacity (Chitwood et al., 2013), providing candidate genes for breeding leaf architectures optimized for specific canopy management or planting densities.

A dominant locus for continuous light tolerance on chromosome 7, identified through RNA-seq, increased yield by 20% in modern tomato hybrid lines carrying the introgressed allele (Velez-Ramirez et al., 2014). This substantial yield gain from manipulating photoperiod response demonstrates that breeding for altered circadian regulation or photoperiod sensitivity could unlock major productivity increases, particularly for protected cultivation systems where light duration can be controlled.

In pepper, RNA-seq analysis of CO2 enrichment responses identified 149 DEGs promoting photosynthesis and yield (Zhang et al., 2021a). Two genes, Capana00g004478 (downregulated) and Capana08g001515 (upregulated) both homologous to AtBGLU45 encoding β-glucosidase involved in carbohydrate and glutathione metabolism, showed opposite expression changes, suggesting antagonistic regulation of carbon partitioning and ROS management under elevated CO2. These findings have implications for breeding cultivars optimized for protected cultivation with CO2 supplementation or for adaptation to rising atmospheric CO2 levels.

Light quality also profoundly affects yield through altered photosynthetic gene expression. RNA-seq analysis of red/blue light treatments in pepper identified 96 genes involved in light reaction pathways, with expression of CPRF2, Paggis, HLIPS, GIGANTEA, LSH1, and FTSH varying by light quality (Tang et al., 2019). Gene Ontology enrichment revealed significant enrichment in xyloglucan:xyloglucosyl transferase activity and carbohydrate metabolic processes closely associated with light responses. The involvement of cell wall modification genes indicates that light quality affects not just photosynthesis but also growth through altered cell expansion, suggesting that LED-based protected cultivation systems could be optimized for both light spectrum and intensity to maximize yield.

6.6 Environmental effects on yield: transcriptional basis of genotype-by-environment interactions

Yield expression depends critically on environmental conditions, with transcriptomics revealing molecular mechanisms underlying genotype-by-environment (G×E) interactions. RNA-seq analysis of soilless substrate-based cultivation versus soil cultivation in tomato detected 476 DEGs, 321 upregulated, 155 downregulated, associated with 10.1% yield increase in soilless systems (Guo et al., 2022). This demonstrates that cultivation practices alter transcriptional programs governing growth and fruit quality, providing molecular explanations for agronomic observations and enabling prediction of which genotypes will respond favorably to specific management practices.

Under salinity stress, Affymetrix microarray analysis identified eight genes, JAZ8 (jasmonate signaling), Polcalcin jun (calcium binding), ER5, ATPase 3 (ion transport), XTH (cell wall modification), WRKY (transcription factor), and TSW12 showing significant positive correlations with total yield (Alsadon et al., 2015). These genes were recommended for incorporation into salinity-tolerant breeding programs, exemplifying how transcriptomics identifies expression-based markers predicting yield maintenance under stress a critical capability for breeding climate-resilient cultivars.

Fruit soluble solids and yield, both commercially important traits, result largely from fruit primary carbon metabolism (Causse et al., 2004). Microarray analysis of six introgression lines measuring Brix and yield revealed that known introgressed alleles altered expression of carbon metabolism enzymes at 20 days after anthesis (Baxter et al., 2005), demonstrating that yield-determining metabolic processes operate during specific developmental windows. This temporal specificity indicates that breeding for altered sugar content or yield requires considering not just which genes are present but when they are expressed during fruit development.

In tobacco, overexpression of GhKTI12 (a Kunitz trypsin inhibitor from cotton) enhanced seed yield and biomass. RNA-seq analysis of transgenic lines identified 518 upregulated and 406 downregulated genes common across overexpression lines (Myat et al., 2022). Most downregulated genes were MADS-box genes regulating flowering time—including MADS6, AP1, AP3, AGL8, AGL6, SEP1, and SEP2, while most upregulated genes governed cell division and differentiation, including RD21, TET8, KTN80, AOX1, AOX2, CP1, and KIC. This coordinated repression of flowering genes with simultaneous activation of cell division programs suggests that yield improvement resulted from extended vegetative growth and enhanced meristematic activity, demonstrating that engineering developmental timing can substantially increase biomass production.

Across diverse yield components and solanaceous species, transcriptomic studies reveal several actionable principles for breeding higher-yielding, climate-adapted cultivars. First, yield determination involves coordinated transcriptional networks spanning morphological development, metabolic regulation, hormone signaling, and stress responses not isolated single-gene effects. Master transcription factors like SlMX1 that regulate multiple yield-determining processes simultaneously represent high-leverage breeding targets enabling multi-trait improvement without negative genetic correlations.

Second, temporal dynamics are critical, developmental decisions determining final yield occur surprisingly early (potato tuber shape at stolon hook stage; carbon metabolism at 20 DAA in tomato fruit), necessitating phenotyping and transcriptome profiling across multiple developmental time points rather than focusing solely on mature organs. This temporal complexity also suggests that breeding for altered developmental timing, as demonstrated by GhKTI12 overexpression could substantially increase productivity by extending growth phases or optimizing source-sink relationships.

Third, yield under stress requires maintaining key physiological processes rather than just activating stress responses, nitrogen-use efficient eggplant maintains photosynthetic gene expression, salt-tolerant tomato maintains ion homeostasis genes, and parthenocarpic lines maintain auxin signaling. This indicates that breeding for yield stability should prioritize genotypes showing minimal transcriptional disruption of core metabolic processes under stress, rather than maximal stress-responsive gene activation.

Fourth, environmental factors (substrate type, light quality, CO2 levels, salinity) substantially alter yield-determining transcriptional programs, with cultivars showing differential transcriptional plasticity. Understanding these G×E interactions through multi-environment transcriptomics enables prediction of which genotypes will perform optimally under specific management practices or environmental conditions, essential for breeding cultivars adapted to specific production systems (protected cultivation, organic systems, saline soils) or future climate scenarios.

Fifth, epitranscriptomic regulation (RNA methylation) and cis-regulatory variation (ENO promoter) both contribute to yield variation, expanding breeding targets beyond traditional protein-coding genes. Genomic selection models can incorporate eQTL and regulatory variants alongside structural variants to capture the full spectrum of functional variation affecting yield.

Finally, the identification of expression-based markers showing consistent correlations with yield across environments (JAZ8, XTH, WRKY in tomato salinity study) provides practical tools for marker-assisted selection or genomic prediction, enabling early-generation selection on molecular profiles before yield itself can be accurately phenotyped. These transcriptomic insights collectively provide a molecular foundation for breeding solanaceous cultivars that achieve the genetic gain acceleration necessary to meet projected food demand while maintaining productivity under the increasingly variable environmental conditions imposed by climate change.

7 Conclusions and future prospects

7.1 Major achievements: from gene discovery to regulatory network elucidation

Transcriptomics has fundamentally transformed solanaceous crop improvement by revealing the molecular architecture underlying complex agronomic traits that have historically resisted genetic dissection. This review has synthesized progress across four major domains: (1) phytochemical biosynthesis controlling nutritional quality and flavor, (2) biotic stress resistance mechanisms spanning diverse pathogen types, (3) abiotic stress tolerance networks enabling climate adaptation, and (4) yield component determination integrating developmental, metabolic, and environmental factors. Several unifying principles emerge from this synthesis that have direct implications for accelerating genetic gain while addressing the consumer-producer trade-offs.

First, transcriptomic studies have consistently revealed that economically important traits whether capsaicinoid pungency, disease resistance, drought tolerance, or fruit size are governed not by single genes but by coordinated regulatory networks orchestrated by master transcription factors. The repeated identification of WRKY, MYB, bZIP, ERF, HSF, and NAC transcription factor families across all trait categories indicates these represent high-leverage breeding targets: manipulating single regulatory genes can trigger cascading expression changes across entire biosynthetic or defense pathways, potentially enabling multi-trait improvement without the negative genetic correlations that plague conventional phenotypic selection. For instance, SlMX1 in tomato simultaneously improves yield, trichome-mediated pest resistance, and stress tolerance, directly addressing the challenge of balancing production imperatives with consumer quality demands.

Second, temporal and spatial dynamics of gene expression are critical for trait manifestation but are frequently overlooked in breeding programs focused on mature phenotypes. Developmental decisions determining potato tuber shape occur at the stolon hook stage, capsaicinoid biosynthesis is restricted to placental tissue, and carbon metabolism during specific fruit development windows (20 days after anthesis in tomato) substantially influences final yield and quality. This temporal-spatial specificity implies that breeding strategies must incorporate stage-specific and tissue-specific transcriptome profiling to capture the full spectrum of functional variation, and that phenotypic selection at inappropriate developmental stages may miss causal genetic variation entirely.

Third, the discovery of regulatory RNAs (lncRNA16397 regulating SlGRX22 in P. infestans resistance; microRNAs modulating cold tolerance in potato), epigenetic modifications (DNA methylation affecting tomato ripening genes), epitranscriptomic regulation (RNA m6A modification improving potato yield), and cis-regulatory variation (ENO promoter mutation increasing tomato fruit size) reveals that breeding programs focused exclusively on protein-coding variation are accessing only a fraction of functional genetic diversity. Expression quantitative trait loci (eQTL) mapping and regulatory variant analysis should be systematically integrated into genomic selection models to capture this hidden genetic variation.

Fourth, transcriptomics has illuminated the molecular basis of genotype-by-environment interactions that confound breeding progress. Storage temperature profoundly affects tomato flavor volatile biosynthesis, nitrogen availability triggers distinct metabolic reprogramming strategies in efficient versus inefficient genotypes, and combined stresses elicit transcriptional responses unpredictable from single-stress treatments. This G×E understanding enables breeding strategies tailored to specific production environments—whether protected cultivation with controlled CO2 and light quality, organic systems with limited nitrogen inputs, or climate-stressed regions with co-occurring heat and drought.

7.2 Current limitations and knowledge gaps

Despite substantial progress, significant gaps constrain the translation of transcriptomic insights into breeding outcomes. First, the vast majority of studies employ bulk RNA-seq from heterogeneous tissues, masking cell-type-specific expression programs that may be functionally critical. For example, understanding why capsaicinoid biosynthesis is restricted to placental cells, or how vascular tissue mediates systemic stress signaling, requires single-cell or spatial transcriptomics approaches that are only beginning to be applied in solanaceous crops. The technical challenges of generating high-quality protoplasts from plants with thick cell walls have limited single-cell RNA-seq adoption, though emerging nuclei isolation protocols may circumvent this barrier.

Second, most transcriptomic studies are descriptive rather than functional, identifying differentially expressed genes without validating causal roles through reverse genetics. While candidate gene lists proliferate in the literature, relatively few have been functionally characterized through knockout, overexpression, or natural allelic variation studies. This validation bottleneck prevents confident prioritization of breeding targets and risks pursuing expression changes that are correlated with but not causal for desirable phenotypes. Integration of transcriptomics with genome editing technologies, particularly multiplex CRISPR approaches enabling simultaneous knockout of multiple candidates could accelerate functional validation.

Third, temporal resolution remains insufficient in most studies. Single-timepoint or coarse time-course analyses miss transient expression changes, oscillatory dynamics, and circadian-regulated processes that may be functionally important. High-resolution time-series experiments, while expensive, are essential for capturing the dynamic transcriptional responses underlying developmental transitions, stress acclimation, and fruit ripening all processes where timing critically determines outcomes.

Fourth, integration across omics layers, transcriptome, proteome, metabolome, epigenome remains limited despite universal recognition of its importance. Transcript abundance often correlates poorly with protein levels due to post-transcriptional and translational regulation, and metabolite concentrations reflect both biosynthesis and turnover rates that are not directly predictable from enzyme expression alone. The AsA inverse correlation between biosynthetic gene expression and accumulation, and the phenylalanine paradox where substrate accumulation doesn’t predict volatile production, exemplify why multi-omics integration is essential for accurate trait prediction and engineering.

Fifth, wild relatives remain underexploited despite harboring resistance and quality alleles absent from cultivated germplasm (Kaur et al., 2023; Malik et al., 2022). While several studies have examined wild eggplant (S. aculeatissimum, S. torvum) and identified massive transcriptional differences underlying superior stress tolerance or pest resistance, systematic transcriptomic characterization across wild Solanum diversity is lacking. Comprehensive wild germplasm transcriptome profiling could reveal novel biosynthetic pathways, regulatory mechanisms, and stress response networks unavailable in domesticated gene pools.

7.3 Emerging technologies and future directions

Several technological advances promise to address current limitations and accelerate trait discovery. Long-read sequencing platforms (PacBio Iso-Seq, Oxford Nanopore direct RNA sequencing) enable full-length transcript characterization, definitively resolving alternative splicing patterns, fusion transcripts, and allele-specific expression in heterozygous and polyploid genomes. Alternative splicing has been shown to generate protein isoforms with distinct or even antagonistic functions, yet remains poorly characterized in solanaceous crops. Comprehensive isoform catalogs for major cultivars and trait-relevant tissues would reveal previously invisible functional variation and enable isoform-specific marker development.

Spatial transcriptomics technologies, including Visium, Slide-seq, and emerging in situ sequencing methods preserve tissue architecture while profiling spatially resolved gene expression, enabling reconstruction of cell-type-specific expression programs without physical cell isolation. For solanaceous crops, spatial transcriptomics could illuminate: (1) gradients of gene expression during fruit ripening that determine quality trait distribution, (2) cell-type-specific pathogen responses revealing why certain tissues resist infection, (3) vascular-specific expression of nutrient transporters governing use efficiency, and (4) localized stress signaling networks determining whether damage remains contained or triggers whole-plant responses.

Multi-omics data integration through machine learning and network inference approaches can identify key regulatory nodes controlling complex traits. Co-expression network analysis has already identified hub genes regulating tomato flavor (Zhu et al., 2022), and extension to multi-layer networks integrating transcriptome-proteome-metabolome-phenome data would enable prediction of which genetic perturbations yield desired phenotypic outcomes. Such predictive models could guide gene editing targets, predict transgressive segregants in breeding populations, and identify optimal allelic combinations for genomic selection.

Pan-genome and structural variation (SV) analyses integrated with transcriptomics represent another frontier. A tomato pan-genome encompassing 238,490 SVs across 100 diverse lines demonstrated that structural variants, including presence/absence variations, copy number changes, and inversions, underlie fruit flavor, size, and yield variation (Alonge et al., 2020). Critically, an SV at the glycosyltransferase NSGT1 was shown to control phenylpropanoid volatile levels affecting flavor. SVs may frequently affect trait variation through altering gene expression (creating or disrupting regulatory elements) rather than protein function, making SV-eQTL mapping essential for comprehensive trait dissection. Moreover, SVs are often missed by short-read resequencing but readily detected by long-read sequencing, suggesting that current genomic selection models ignoring structural variation are capturing incomplete genetic architectures.

7.4 Translational applications: from transcriptome to improved cultivars

Realizing the breeding impact of transcriptomic discoveries requires systematic implementation in breeding programs, which remains limited. Expression-based markers, genes whose expression levels correlate consistently with trait values across environments could enable early-generation selection before phenotypes fully manifest. For example, identifying seedling-stage expression signatures predicting mature plant drought tolerance or fruit quality would dramatically accelerate breeding cycles. Developing robust expression-based prediction models requires large-scale, multi-environment transcriptome profiling of breeding populations, a substantial investment but potentially transformative for genetic gain rates.

Genomic selection models incorporating expression data alongside DNA sequence variation improve prediction accuracy for complex traits, particularly under novel environmental conditions not represented in training populations. eQTL mapping identifies variants affecting gene expression, which can be weighted in genomic prediction models based on biological priors (e.g., known pathway involvement). This functionally informed genomic selection could improve prediction of trait performance in untested environments or under climate scenarios not yet experienced, addressing the critical need for breeding cultivars adapted to future rather than historical climates.

Gene editing technologies, particularly base editing and prime editing enabling precise nucleotide changes coupled with transcriptomic identification of causal regulatory variants, enable direct engineering of desired expression patterns without creating transgenic plants. The ENO promoter mutation increasing tomato fruit size, or the SlAREB1 expression-enhancing variants conferring drought tolerance, could be engineered into elite backgrounds through non-transgenic editing, circumventing regulatory barriers to GMO deployment while leveraging transcriptomic insights.

Speed breeding protocols combined with transcriptome-guided selection could dramatically accelerate genetic gain. By profiling expression signatures at early growth stages under controlled stress treatments, breeders could select superior genotypes within months rather than years. This is particularly valuable for traits like multi-stress tolerance where field phenotyping requires multiple seasons and locations but transcriptional signatures might predict field performance from controlled environment tests.

7.5 Addressing the multi-trait breeding challenge

The introduction highlighted the critical challenge of simultaneously improving yield, stress resilience, and fruit quality traits, often showing negative genetic correlations when selected phenotypically. Transcriptomics has revealed why these correlations exist and how they might be broken. Stress tolerance often involves activating resource-expensive defense pathways (massive transcriptional reprogramming in V. dahliae-resistant eggplant) that divert resources from growth and fruit development. However, identification of genotypes maintaining yield under stress through minimal transcriptional disruption of core metabolic processes, rather than maximal stress response activation, suggests an alternative breeding strategy: selecting for stress resilience (maintaining function) rather than stress resistance (activating defenses).

Similarly, the inverse relationship between tomato refrigeration-induced flavor loss and shelf-life reflects transcriptional suppression of volatile biosynthesis at cold storage temperatures. Understanding this mechanistic basis enables targeted breeding: identifying or engineering cultivars that maintain volatile biosynthesis at low temperatures would simultaneously improve flavor and enable cold storage, breaking the apparent trade-off. Transcriptomics thus provides the mechanistic understanding necessary to design breeding strategies that overcome negative correlations through targeting underlying regulatory mechanisms rather than treating traits as independent.

Climate change imposes an urgent imperative to accelerate genetic gain in solanaceous crops while simultaneously improving stress resilience and maintaining quality attributes demanded by consumers and processors. Transcriptomics has transitioned from a gene discovery tool to a systems-level technology illuminating regulatory networks, developmental dynamics, environmental plasticity, and multi-trait integration. The path forward requires: (1) systematic application of emerging technologies (spatial transcriptomics, long-read sequencing) to address current knowledge gaps, (2) functional validation of candidate genes through genome editing, (3) multi-omics integration for predictive trait modeling, (4) large-scale expression profiling of breeding populations for genomic selection implementation, and (5) exploration of wild germplasm to access novel regulatory mechanisms.

Critically, realizing the breeding impact of transcriptomic insights demands closer collaboration between molecular biologists generating transcriptome data and breeders making selection decisions. Transcriptomic discoveries must be translated into actionable selection tools expression-based markers, eQTL-informed genomic selection, and guide RNA designs for precision editing. The substantial investment in transcriptome profiling made by the research community over the past two decades positions solanaceous crop breeding to meet future challenges, but only if these molecular insights are systematically integrated into breeding practice. The tools, knowledge, and opportunity exist; what remains is implementation to deliver the climate-resilient, high-quality, productive cultivars that global food security demands.

Statements

Author contributions

XQ: Conceptualization, Writing – original draft. BJ: Conceptualization, Investigation, Writing – original draft. JQ: Conceptualization, Investigation, Writing – original draft. JMS: Conceptualization, Writing – original draft. IK: Conceptualization, Writing – original draft. RH: Writing – original draft. JDS: Writing – original draft. LS: Writing – original draft. HZ: Writing – review & editing. XY: Writing – review & editing. XL: Funding acquisition, Writing – review & editing. GK: Conceptualization, Writing – original draft. HL: Conceptualization, Funding acquisition, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This project was funded by grants from the Postdoctoral Research Foundation of Anhui Agricultural University (03087003), the National Natural Science Foundation of China (32402592), Anhui Agricultural University High-level Talent Introduction Funds (rc322302) and Anhui Overseas Talent Recruitment Program.

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    Agarwal P. Parida S. K. Mahto A. Das S. Mathew I. E. Malik N. et al . (2014). Expanding frontiers in plant transcriptomics in aid of functional genomics and molecular breeding. Biotechnol. J.9, 14801492. doi: 10.1002/biot.201400063

  • 2

    Aldridge S. Teichmann S. A. (2020). Single cell transcriptomics comes of age. Nat. Commun.11, 4307. doi: 10.1038/s41467-020-18158-5

  • 3

    Aliche E. B. Gengler T. Hoendervangers I. Oortwijn M. Bachem C. W. B. Borm T. et al . (2021). Transcriptomic responses of potato to drought stress. Potato Res. 65, 289305. doi: 10.1007/s11540-021-09527-8

  • 4

    Alonge M. Wang X. Benoit M. Soyk S. Pereira L. Zhang L. et al . (2020). Major impacts of widespread structural variation on gene expression and crop improvement in tomato. Cell182, 145161.e23. doi: 10.1016/j.cell.2020.05.021

  • 5

    Alós E. Rodrigo M. J. Zacarías L. (2013). Transcriptomic analysis of genes involved in the biosynthesis, recycling and degradation of L-ascorbic acid in pepper fruits (Capsicum annuum L.). Plant Sci.207, 211. doi: 10.1016/j.plantsci.2013.02.007

  • 6

    Alsadon A. A. Ibrahim A. A. Wahb-Allah M. A. Ali A. A. M. Sadder M. T. (2015). Tomato under salinity stress: correlation between growth and yield components and responsive genes. Acta Hortic.1081, 111119. doi: 10.17660/ActaHortic.2015.1081.11

  • 7

    Altaf M. A. Mandal S. Behera B. Mangal V. Naz S. Kumar R. et al . (2022). Salinity stress tolerance in solanaceous crops: current understanding and its prospects in genome editing. J. Plant Growth Regul. 42, 40204036. doi: 10.1007/s00344-022-10890-0

  • 8

    Alwine J. C. Kemp D. J. Stark G. R. (1977). Method for detection of specific RNAs in agarose gels by transfer to diazobenzyloxymethyl-paper and hybridization with DNA probes. Proc. Natl. Acad. Sci. U.S.A.74, 53505354. doi: 10.1073/pnas.74.12.5350

  • 9

    Amiryousefi A. Hyvönen J. Poczai P. (2018). The chloroplast genome sequence of bittersweet (Solanum dulcamara): Plastid genome structure evolution in Solanaceae. PLoS One13, e0196069. doi: 10.1371/journal.pone.0196069

  • 10

    Basim H. Basim E. Tombuloglu H. Unver T. (2021). Comparative transcriptome analysis of resistant and cultivated tomato lines in response to Clavibacter michiganensis subsp. michiganensis. Genomics113, 24552467. doi: 10.1016/j.ygeno.2021.05.033

  • 11

    Baxter C. J. Sabar M. Quick W. P. Sweetlove L. J. (2005). Comparison of changes in fruit gene expression in tomato introgression lines provides evidence of genome-wide transcriptional changes and reveals links to mapped QTLs and described traits. J. Exp. Bot.56, 15911604. doi: 10.1093/jxb/eri154

  • 12

    Blanco-Ulate B. Vincenti E. Powell A. L. T. Cantu D. (2013). Tomato transcriptome and mutant analyses suggest a role for plant stress hormones in the interaction between fruit and Botrytis cinerea. Front. Plant Sci.4. doi: 10.3389/fpls.2013.00142

  • 13

    Bobokalonov J. Liu Y. Shahrin T. Liu L. (2018). Transcriptomic analysis on the regulation of tomato ripening by the ethylene inhibitor 1-methylcyclopropene. JPS7, 49. doi: 10.5539/jps.v7n2p49

  • 14

    Brouwer S. M. Brus-Szkalej M. Saripella G. V. Liang D. Liljeroth E. Grenville-Briggs L. J. (2021). Transcriptome analysis of potato infected with the necrotrophic pathogen Alternaria solani. Plants10, 2212. doi: 10.3390/plants10102212

  • 15

    Cabezas-Serrano A. B. Amodio M. L. Cornacchia R. Rinaldi R. Colelli G. (2009). Suitability of five different potato cultivars (Solanum tuberosum L.) to be processed as fresh-cut products. Postharvest Biol. Technol.53, 138144. doi: 10.1016/j.postharvbio.2009.03.009

  • 16

    Cammarano D. Jamshidi S. Hoogenboom G. Ruane A. C. Niyogi D. Ronga D. (2022). Processing tomato production is expected to decrease by 2050 due to the projected increase in temperature. Nat. Food3, 437444. doi: 10.1038/s43016-022-00521-y

  • 17

    Campbell R. Pont S. D. A. Morris J. A. McKenzie G. Sharma S. K. Hedley P. E. et al . (2014). Genome-wide QTL and bulked transcriptomic analysis reveals new candidate genes for the control of tuber carotenoid content in potato (Solanum tuberosum L.). Theor. Appl. Genet.127, 19171933. doi: 10.1007/s00122-014-2349-0

  • 18

    Cara B. Giovannoni J. J. (2008). Molecular biology of ethylene during tomato fruit development and maturation. Plant Sci.175, 106113. doi: 10.1016/j.plantsci.2008.03.021

  • 19

    Carvalho Lemos V. Reimer J. J. Wormit A. (2019). Color for life: Biosynthesis and distribution of phenolic compounds in pepper (Capsicum annuum). Agriculture9, 81. doi: 10.3390/agriculture9040081

  • 20

    Causse M. Duffe P. Gomez M. C. Buret M. Damidaux R. Zamir D. et al . (2004). A genetic map of candidate genes and QTLs involved in tomato fruit size and composition. J. Exp. Bot.55, 16711685. doi: 10.1093/jxb/erh207

  • 21

    Chandrasekar S. Natarajan P. Mhatre P. H. Mahajan M. Nivitha S. Palanisamy V. E. et al . (2022). RNA-Seq of cyst nematode infestation of potato (Solanum tuberosum L.): A comparative transcriptome analysis of resistant and susceptible cultivars. Plants11, 1008. doi: 10.3390/plants11081008

  • 22

    Chen L.-M. Li X.-W. He T.-J. Li P.-J. Liu Y. Zhou S.-X. et al . (2021). Comparative biochemical and transcriptome analyses in tomato and eggplant reveal their differential responses to Tuta absoluta infestation. Genomics113, 21082121. doi: 10.1016/j.ygeno.2021.05.002

  • 23

    Chen T. Lv Y. Zhao T. Li N. Yang Y. Yu W. et al . (2013). Comparative transcriptome profiling of a resistant vs. susceptible tomato (Solanum lycopersicum) cultivar in response to infection by tomato yellow leaf curl virus. PloS One8, e80816. doi: 10.1371/journal.pone.0080816

  • 24

    Chen N. Yu B. Dong R. Lei J. Chen C. Cao B. (2018). RNA-Seq-derived identification of differential transcription in the eggplant (Solanum melongena) following inoculation with bacterial wilt. Gene644, 137147. doi: 10.1016/j.gene.2017.11.003

  • 25

    Chitwood D. H. Kumar R. Headland L. R. Ranjan A. Covington M. F. Ichihashi Y. et al . (2013). A quantitative genetic basis for leaf morphology in a set of precisely defined tomato introgression lines. Plant Cell25, 24652481. doi: 10.1105/tpc.113.112391

  • 26

    Cruz-Mendívil A. López-Valenzuela J. A. Calderón-Vázquez C. L. Vega-García M. O. Reyes-Moreno C. Valdez-Ortiz A. (2015). Transcriptional changes associated with chilling tolerance and susceptibility in ‘Micro-Tom’ tomato fruit using RNA-Seq. Postharvest Biol. Technol.99, 141151. doi: 10.1016/j.postharvbio.2014.08.009

  • 27

    Cui J. Luan Y. Jiang N. Bao H. Meng J. (2017). Comparative transcriptome analysis between resistant and susceptible tomato allows the identification of lncRNA16397 conferring resistance to Phytophthora infestans by co-expressing glutaredoxin. Plant J.89, 577589. doi: 10.1111/tpj.13408

  • 28

    Cui J. Xu P. Meng J. Li J. Jiang N. Luan Y. (2018). Transcriptome signatures of tomato leaf induced by Phytophthora infestans and functional identification of transcription factor SpWRKY3. Theor. Appl. Genet.131, 787800. doi: 10.1007/s00122-017-3035-9

  • 29

    D’Esposito D. Ferriello F. Molin A. D. Diretto G. Sacco A. Minio A. et al . (2017). Unraveling the complexity of transcriptomic, metabolomic and quality environmental response of tomato fruit. BMC Plant Biol.17, 66. doi: 10.1186/s12870-017-1008-4

  • 30

    Dal Cin V. Tieman D. M. Tohge T. McQuinn R. de Vos R. C. H. Osorio S. et al . (2011). Identification of genes in the phenylalanine metabolic pathway by ectopic expression of a MYB transcription factor in tomato fruit. Plant Cell23, 27382753. doi: 10.1105/tpc.111.086975

  • 31

    Darshan K. Aggarwal R. Bashyal B. M. Singh J. Shanmugam V. Gurjar M. S. et al . (2020). Transcriptome profiling provides insights into potential antagonistic mechanisms involved in Chaetomium globosum against Bipolaris sorokiniana. Front. Microbiol.11. doi: 10.3389/fmicb.2020.578115

  • 32

    Demirel U. Morris W. L. Ducreux L. J. M. Yavuz C. Asim A. Tindas I. et al . (2020). Physiological, biochemical, and transcriptional responses to single and combined abiotic stress in stress-tolerant and stress-sensitive potato genotypes. Front. Plant Sci.11. doi: 10.3389/fpls.2020.00169

  • 33

    Docimo T. Francese G. Ruggiero A. Batelli G. De Palma M. Bassolino L. et al . (2016). Phenylpropanoids accumulation in eggplant fruit: Characterization of biosynthetic genes and regulation by a MYB transcription factor. Front. Plant Sci.6. doi: 10.3389/fpls.2015.01233

  • 34

    Egan A. N. Schlueter J. Spooner D. M. (2012). Applications of next-generation sequencing in plant biology. Am. J. Bot.99, 175185. doi: 10.3732/ajb.1200020

  • 35

    Ewas M. Harlina P. W. Shahzad R. Khames E. Ali F. Nishawy E. et al . (2022). Constitutive expression of SlMX1 gene improves fruit yield and quality, health-promoting compounds, fungal resistance and delays ripening in transgenic tomato plants. J. @ Plant Interact.17, 517536. doi: 10.1080/17429145.2022.2066730

  • 36

    Fan G. Wang Q. Xu J. Chen N. Zhu W. Duan S. et al . (2022). Fine mapping and candidate gene prediction of tuber shape controlling ro locus based on integrating genetic and transcriptomic analyses in potato. Int. J. Mol. Sci.23, 1470. doi: 10.3390/ijms23031470

  • 37

    Fan Y. Zhang W. Kang Y. Shi M. Yang X. Yu H. et al . (2021). Physiological and dynamic transcriptome analysis of two potato varieties reveal response of lignin and MAPK signal to dry rot caused by Fusarium sulphureum. Sci. Hortic.289, 110470. doi: 10.1016/j.scienta.2021.110470

  • 38

    Fiesel P. D. Parks H. M. Last R. L. Barry C. S. (2022). Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae. Nat. Prod. Rep.39, 14381464. doi: 10.1039/d2np00003b

  • 39

    Fu A. Zheng Y. Lv Y. Watkins C. B. Bai C. Ma L. et al . (2022). Multi-omics analysis reveals specific modifications associated with reduced chilling injury in bell pepper fruit by methyl jamonate. Postharvest Biol. Technol.185, 111799. doi: 10.1016/j.postharvbio.2021.111799

  • 40

    Gao L. Gonda I. Sun H. Ma Q. Bao K. Tieman D. M. et al . (2019). The tomato pan-genome uncovers new genes and a rare allele regulating fruit flavor. Nat. Genet.51, 10441051. doi: 10.1038/s41588-019-0410-2

  • 41

    Gao S. Wang F. Niran J. Li N. Yin Y. Yu C. et al . (2021). Transcriptome analysis reveals defense-related genes and pathways against Xanthomonas campestris pv. vesicatoria in pepper (Capsicum annuum L.). PLoS One16, e0240279. doi: 10.1371/journal.pone.0240279

  • 42

    Garbowicz K. Liu Z. Alseekh S. Tieman D. Taylor M. Kuhalskaya A. et al . (2018). Quantitative trait loci analysis identifies a prominent gene involved in the production of fatty acid-derived flavor volatiles in tomato. Mol. Plant11, 11471165. doi: 10.1016/j.molp.2018.06.003

  • 43

    Ghareeb H. Bozsó Z. Ott P. G. Repenning C. Stahl F. Wydra K. (2011). Transcriptome of silicon-induced resistance against Ralstonia solanacearum in the silicon non-accumulator tomato implicates priming effect. Physiol. Mol. Plant Pathol.75, 8389. doi: 10.1016/j.pmpp.2010.11.004

  • 44

    Gong Z. Xiong L. Shi H. Yang S. Herrera-Estrella L. R. Xu G. et al . (2020). Plant abiotic stress response and nutrient use efficiency. Sci. China Life Sci.63, 635674. doi: 10.1007/s11427-020-1683-x

  • 45

    González-Gordo S. Rodríguez-Ruiz M. López-Jaramillo J. Muñoz-Vargas M. A. Palma J. M. Corpas F. J. (2022). Nitric oxide (NO) differentially modulates the ascorbate peroxidase (APX) isozymes of sweet pepper (Capsicum annuum L.) fruits. Antioxidants (Basel)11, 765. doi: 10.3390/antiox11040765

  • 46

    Góra-Sochacka A. Więsyk A. Fogtman A. Lirski M. Zagórski-Ostoja W. (2019). Root transcriptomic analysis reveals global changes induced by systemic infection of Solanum lycopersicum with mild and severe variants of potato spindle tuber viroid. Viruses11, 992. doi: 10.3390/v11110992

  • 47

    Guevara L. Domínguez-Anaya M.Á. Ortigosa A. González-Gordo S. Díaz C. Vicente F. et al . (2021). Identification of compounds with potential therapeutic uses from sweet pepper (Capsicum annuum L.) fruits and their modulation by nitric oxide (NO). Int. J. Mol. Sci.22, 4476. doi: 10.3390/ijms22094476

  • 48

    Guo J. Dong L. Kandel S. L. Jiao Y. Shi L. Yang Y. et al . (2022). Transcriptomic and Metabolomic Analysis Provides Insights into the Fruit Quality and Yield Improvement in Tomato under Soilless Substrate-Based Cultivation. Agronomy12, 923. doi: 10.3390/agronomy12040923

  • 49

    Gyetvai G. Sønderkær M. Göbel U. Basekow R. Ballvora A. Imhoff M. et al . (2012). The transcriptome of compatible and incompatible interactions of potato (Solanum tuberosum) with Phytophthora infestans revealed by DeepSAGE analysis. PLoS One7, e31526. doi: 10.1371/journal.pone.0031526

  • 50

    He M. Zhou Y. Ye G. Zheng J. Meng Y. Wang J. et al . (2021). Serial transcriptome analysis reveals genes associated with late blight resistance in potato cultivar Qingshu 9. Agronomy11, 1919. doi: 10.3390/agronomy11101919

  • 51

    Ishihara T. Mitsuhara I. Takahashi H. Nakaho K. (2012). Transcriptome analysis of quantitative resistance-specific response upon Ralstonia solanacearum infection in tomato. PLoS One7, e46763. doi: 10.1371/journal.pone.0046763

  • 52

    Jiang M. Liu Y. Ren L. Lian H. Chen H. (2016). Molecular cloning and characterization of anthocyanin biosynthesis genes in eggplant (Solanum melongena L.). Acta Physiol. Plant38, 163. doi: 10.1007/s11738-016-2172-0

  • 53

    Kahlau S. Bock R. (2008). Plastid transcriptomics and translatomics of tomato fruit development and chloroplast-to-chromoplast differentiation: chromoplast gene expression largely serves the production of a single protein. Plant Cell20, 856874. doi: 10.1105/tpc.107.055202

  • 54

    Kang W.-H. Sim Y. M. Koo N. Nam J.-Y. Lee J. Kim N. et al . (2020). Transcriptome profiling of abiotic responses to heat, cold, salt, and osmotic stress of Capsicum annuum L. Sci. Data7, 17. doi: 10.1038/s41597-020-0352-7

  • 55

    Kantar M. Ewing P. Bančič J. Fumia N. Garba I. Jamshidi S. et al . (2025). Computational design for more engaged, impactful, and dynamic agricultural research. Crop Sci.65, e70034. doi: 10.1002/csc2.70034

  • 56

    Karlova R. Rosin F. M. Busscher-Lange J. Parapunova V. Do P. T. Fernie A. R. et al . (2011). Transcriptome and metabolite profiling show that APETALA2a is a major regulator of tomato fruit ripening. Plant Cell23, 923941. doi: 10.1105/tpc.110.081273

  • 57

    Kaur G. Abugu M. Tieman D. (2023). The dissection of tomato flavor: biochemistry, genetics, and omics. Front. Plant Sci.14. doi: 10.3389/fpls.2023.1144113

  • 58

    Kaur G. Lujan P. Sanogo S. Steiner R. Puppala N. (2019). Assessing in vitro efficacy of certain fungicides to control Sclerotinia sclerotiorum in peanut. Arch. Of Phytopathol. And Plant Prot.52, 184199. doi: 10.1080/03235408.2019.1603350

  • 59

    Kaur I. Watts S. Raya C. Raya J. Kariyat R. (2022). “ Surface warfare: plant structural defenses challenge caterpillar feeding,” in Caterpillars in the middle: tritrophic interactions in a changing world Fascinating Life Sciences. Eds. MarquisR. J.KopturS. ( Springer International Publishing, Cham), 6592. doi: 10.1007/978-3-030-86688-4_3

  • 60

    Keating B. A. Herrero M. Carberry P. S. Gardner J. Cole M. B. (2014). Food wedges: Framing the global food demand and supply challenge towards 2050. Glob. Food Sec.3, 125132. doi: 10.1016/j.gfs.2014.08.004

  • 61

    Keller M. Hu Y. Mesihovic A. Fragkostefanakis S. Schleiff E. Simm S. (2017). Alternative splicing in tomato pollen in response to heat stress. DNA Res.24, 205217. doi: 10.1093/dnares/dsw051

  • 62

    Khan R. Zhou P. Ma X. Zhou L. Wu Y. Ullah Z. et al . (2019). Transcriptome Profiling, Biochemical and Physiological Analyses Provide New Insights towards Drought Tolerance in Nicotiana tabacum L. Genes (Basel)10, 1041. doi: 10.3390/genes10121041

  • 63

    Kim N. Lee J. Yeom S. I. Kang N. J. Kang W. H. (2024). The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper. Sci. Data11, 381. doi: 10.1038/s41597-024-03239-7

  • 64

    Kim S. Park M. Yeom S.-I. Kim Y.-M. Lee J. M. Lee H.-A. et al . (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet.46, 270278. doi: 10.1038/ng.2877

  • 65

    Klebanov L. Yakovlev A. (2007). How high is the level of technical noise in microarray data? Biol. Direct2, 9. doi: 10.1186/1745-6150-2-9

  • 66

    Kumar R. Rai A. Rai A. C. Singh V. K. Singh M. Singh P. M. et al . (2022). De novo assembly, differential gene expression and pathway analyses for anthracnose resistance in chilli (Capsicum annuum L.). J. Plant Biochem. Biotechnol.31, 124138. doi: 10.1007/s13562-021-00668-y

  • 67

    Kushwaha N. Sahu P. P. Prasad M. Chakraborty S. (2015). Chilli leaf curl virus infection highlights the differential expression of genes involved in protein homeostasis and defense in resistant chilli plants. Appl. Microbiol. Biotechnol.99, 47574770. doi: 10.1007/s00253-015-6415-6

  • 68

    Lei B. Zhao X.-H. Zhang K. Zhang J. Ren W. Ren Z. et al . (2013). Comparative transcriptome analysis of tobacco (Nicotiana tabacum) leaves to identify aroma compound-related genes expressed in different cultivated regions. Mol. Biol. Rep.40, 345357. doi: 10.1007/s11033-012-2067-0

  • 69

    Li Y. Chang D. Zhang X. Shi H. Yang H. (2021b). RNA-Seq, physiological, and biochemical analysis of burley tobacco response to nitrogen deficiency. Sci. Rep.11, 16802. doi: 10.1038/s41598-021-93363-w

  • 70

    Li J. Chitwood-Brown J. Kaur G. Labate J. A. Vallad G. E. Lee T. G. et al . (2022a). Novel Sources of Resistance to Fusarium oxysporum f. sp. lycopersici Race 3 Among Solanum pennellii Accessions. J. Amer. Soc Hortic. Sci.147, 3544. doi: 10.21273/JASHS05080-21

  • 71

    Li P. Fan R. Peng Z. Qing Y. Fang Z. (2021a). Transcriptome analysis of resistance mechanism to potato wart disease. Open Life Sci.16, 475481. doi: 10.1515/biol-2021-0045

  • 72

    Li J. He Y.-J. Zhou L. Liu Y. Jiang M. Ren L. et al . (2018a). Transcriptome profiling of genes related to light-induced anthocyanin biosynthesis in eggplant (Solanum melongena L.) before purple color becomes evident. BMC Genomics19, 201. doi: 10.1186/s12864-018-4587-z

  • 73

    Li Y. Hu X. Chen J. Wang W. Xiong X. He C. (2017). Integrated mRNA and microRNA transcriptome analysis reveals miRNA regulation in response to PVA in potato. Sci. Rep.7, 16925. doi: 10.1038/s41598-017-17059-w

  • 74

    Li W. Lu X. Li J. (2022b). The effect of organic nutrient solution on flavor in ripe cherry tomato fruit—Transcriptome and metabolomic analyses. Environ. Exp. Bot.194, 104721. doi: 10.1016/j.envexpbot.2021.104721

  • 75

    Li X. Tieman D. Liu Z. Chen K. Klee H. J. (2020). Identification of a lipase gene with a role in tomato fruit short-chain fatty acid-derived flavor volatiles by genome-wide association. Plant J.104, 631644. doi: 10.1111/tpj.14951

  • 76

    Li K. Wu G. Li M. Ma M. Du J. Sun M. et al . (2018b). Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus. Virol. J.15, 138. doi: 10.1186/s12985-018-1044-1

  • 77

    Li X. Xing X. Tian P. Zhang M. Huo Z. Zhao K. et al . (2018c). Comparative transcriptome profiling reveals defense-related genes against Meloidogyne incognita invasion in tobacco. Molecules23, 2081. doi: 10.3390/molecules23082081

  • 78

    Liao H. Wang Q. Zhang N. Fu Y. Wu G. Ren X. et al . (2021). High-Throughput MicroRNA and mRNA Sequencing Reveals that MicroRNAs may be Involved in Peroxidase-Mediated Cold Tolerance in Potato. Plant Mol. Biol. Rep.39, 577594. doi: 10.1007/s11105-020-01272-5

  • 79

    Lin Q. Xie Y. Guan W. Duan Y. Wang Z. Sun C. (2019). Combined transcriptomic and proteomic analysis of cold stress induced sugar accumulation and heat shock proteins expression during postharvest potato tuber storage. Food Chem.297, 124991. doi: 10.1016/j.foodchem.2019.124991

  • 80

    Liu B. Kong L. Zhang Y. Liao Y. (2021a). Gene and metabolite integration analysis through transcriptome and metabolome brings new insight into heat stress tolerance in potato (Solanum tuberosum L.). Plants10, 103. doi: 10.3390/plants10010103

  • 81

    Liu S. Li W. Wu Y. Chen C. Lei J. (2013). De novo transcriptome assembly in chili pepper (Capsicum frutescens) to identify genes involved in the biosynthesis of capsaicinoids. PLoS One8, e48156. doi: 10.1371/journal.pone.0048156

  • 82

    Liu Y. Lin-Wang K. Deng C. Warran B. Wang L. Yu B. et al . (2015). Comparative transcriptome analysis of white and purple potato to identify genes involved in anthocyanin biosynthesis. PLoS One10, e0129148. doi: 10.1371/journal.pone.0129148

  • 83

    Liu Y. Lv J. Liu Z. Wang J. Yang B. Chen W. et al . (2020). Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in pepper fruit (Capsicum annuum L.). Food Chem.306, 125629. doi: 10.1016/j.foodchem.2019.125629

  • 84

    Liu X. Zhang A. Zhao J. Shang J. Zhu Z. Wu X. et al . (2021b). Transcriptome profiling reveals potential genes involved in browning of fresh-cut eggplant (Solanum melongena L.). Sci. Rep.11, 16081. doi: 10.1038/s41598-021-94831-z

  • 85

    Liu D. Zhao Q. Cheng Y. Li D. Jiang C. Cheng L. et al . (2019). Transcriptome analysis of two cultivars of tobacco in response to Cucumber mosaic virus infection. Sci. Rep.9, 3124. doi: 10.1038/s41598-019-39734-w

  • 86

    López-Serrano L. Calatayud Á. López-Galarza S. Serrano R. Bueso E. (2021). Uncovering salt tolerance mechanisms in pepper plants: a physiological and transcriptomic approach. BMC Plant Biol.21, 169. doi: 10.1186/s12870-021-02938-2

  • 87

    Lowe R. Shirley N. Bleackley M. Dolan S. Shafee T. (2017). Transcriptomics technologies. PLoS Comput. Biol.13, e1005457. doi: 10.1371/journal.pcbi.1005457

  • 88

    Lu J. Du Z.-X. Kong J. Chen L.-N. Qiu Y.-H. Li G.-F. et al . (2012). Transcriptome analysis of Nicotiana tabacum infected by Cucumber mosaic virus during systemic symptom development. PLoS One7, e43447. doi: 10.1371/journal.pone.0043447

  • 89

    Lucas J. A. (2011). Advances in plant disease and pest management. J. Agric. Sci.149, 91114. doi: 10.1017/S0021859610000997

  • 90

    Lv L. L. Feng X. F. Li W. Li K. (2019). High temperature reduces peel color in eggplant (Solanum melongena) as revealed by RNA-seq analysis. Genome62, 503512. doi: 10.1139/gen-2019-0021

  • 91

    Ma X. Li Y. Gai W.-X. Li C. Gong Z.-H. (2021b). The CaCIPK3 gene positively regulates drought tolerance in pepper. Hortic. Res.8, 216. doi: 10.1038/s41438-021-00651-7

  • 92

    Ma L. Wang Q. Li L. Grierson D. Yuan S. Zheng S. et al . (2021a). UV-C irradiation delays the physiological changes of bell pepper fruit during storage. Postharvest Biol. Technol.180, 111506. doi: 10.1016/j.postharvbio.2021.111506

  • 93

    Mahmoud G. A.-E. (2021). “ Biotic stress to legumes: fungal diseases as major biotic stress factor,” in Sustainable agriculture reviews 51: legume agriculture and biotechnology vol 2 Sustainable Agriculture Reviews. Eds. GuleriP.KumarV.LichtfouseE. ( Springer International Publishing, Cham), 181212. doi: 10.1007/978-3-030-68828-8_7

  • 94

    Malik P. Huang M. Neelam K. Bhatia D. Kaur R. Yadav B. et al . (2022). Genotyping-by-sequencing based investigation of population structure and genome wide association studies for seven agronomically important traits in a set of 346 Oryza rufipogon accessions. Rice15, 37. doi: 10.1186/s12284-022-00582-4

  • 95

    Marín A. Ferreres F. Tomás-Barberán F. A. Gil M. I. (2004). Characterization and quantitation of antioxidant constituents of sweet pepper (Capsicum annuum L.). J. Agric. Food Chem.52, 38613869. doi: 10.1021/jf0497915

  • 96

    Martínez-López L. A. Ochoa-Alejo N. Martínez O. (2014). Dynamics of the chili pepper transcriptome during fruit development. BMC Genomics15, 143. doi: 10.1186/1471-2164-15-143

  • 97

    Mathieu S. Cin V. D. Fei Z. Li H. Bliss P. Taylor M. G. et al . (2009). Flavour compounds in tomato fruits: identification of loci and potential pathways affecting volatile composition. J. Exp. Bot.60, 325337. doi: 10.1093/jxb/ern294

  • 98

    Matsumura H. Nirasawa S. Terauchi R. (1999). Technical advance: transcript profiling in rice (Oryza sativa L.) seedlings using serial analysis of gene expression (SAGE). Plant J.20, 719726. doi: 10.1046/j.1365-313x.1999.00640.x

  • 99

    Mauceri A. Abenavoli M. R. Toppino L. Panda S. Mercati F. Aci M. M. et al . (2021). Transcriptomics reveal new insights into molecular regulation of nitrogen use efficiency in Solanum melongena. J. Exp. Bot.72, 42374253. doi: 10.1093/jxb/erab121

  • 100

    Meng H. Sun M. Jiang Z. Liu Y. Sun Y. Liu D. et al . (2021). Comparative transcriptome analysis reveals resistant and susceptible genes in tobacco cultivars in response to infection by Phytophthora nicotianae. Sci. Rep.11, 809. doi: 10.1038/s41598-020-80280-7

  • 101

    Merriman B. Ion Torrent R&D Team Rothberg J. M. (2012). Progress in ion torrent semiconductor chip based sequencing. Electrophoresis33, 33973417. doi: 10.1002/elps.201200424

  • 102

    Meyers B. C. Vu T. H. Tej S. S. Ghazal H. Matvienko M. Agrawal V. et al . (2004). Analysis of the transcriptional complexity of Arabidopsis thaliana by massively parallel signature sequencing. Nat. Biotechnol.22, 10061011. doi: 10.1038/nbt992

  • 103

    Mullis K. Faloona F. Scharf S. Saiki R. Horn G. Erlich H. (1986). Specific enzymatic amplification of DNA in vitro: the polymerase chain reaction. Cold Spring Harb. Symp. Quant. Biol.51 Pt 1, 263273. doi: 10.1101/SQB.1986.051.01.032

  • 104

    Myat A. A. Zhou Y. Gao Y. Zhao X. Liang C. Abid M. A. et al . (2022). Overexpression of ghkti12 enhances seed yield and biomass production in nicotiana tabacum. Genes (Basel)13, 426. doi: 10.3390/genes13030426

  • 105

    Nicotra A. B. Leigh A. Boyce C. K. Jones C. S. Niklas K. J. Royer D. L. et al . (2011). The evolution and functional significance of leaf shape in the angiosperms. Funct. Plant Biol.38, 535. doi: 10.1071/FP11057

  • 106

    Nie S. Yang Z. Xiao C. Yang R. Wang D. (2022). Physiological and transcriptomic analyses of the effects of slbri1 expression levels on the drought tolerance of tomato seedlings. J. Plant Growth Regul. 42, 795809. doi: 10.1007/s00344-022-10587-4

  • 107

    Olmstead R. G. Bohs L. Migid H. A. Santiago-Valentin E. Garcia V. F. Collier S. M. (2008). A molecular phylogeny of the Solanaceae. Taxon57, 11591181. doi: 10.1002/tax.574010

  • 108

    Orellana S. Yañez M. Espinoza A. Verdugo I. González E. Ruiz-Lara S. et al . (2010). The transcription factor SlAREB1 confers drought, salt stress tolerance and regulates biotic and abiotic stress-related genes in tomato. Plant Cell Environ.33, 21912208. doi: 10.1111/j.1365-3040.2010.02220.x

  • 109

    Padmanabhan C. Ma Q. Shekasteband R. Stewart K. S. Hutton S. F. Scott J. W. et al . (2019). Comprehensive transcriptome analysis and functional characterization of PR-5 for its involvement in tomato Sw-7 resistance to tomato spotted wilt tospovirus. Sci. Rep.9, 7673. doi: 10.1038/s41598-019-44100-x

  • 110

    Page A. M. L. Chapman M. A. (2021). Identifying genomic regions targeted during eggplant domestication using transcriptome data. J. Hered.112, 519525. doi: 10.1093/jhered/esab035

  • 111

    Piedra-Aguilera Á. Jiao C. Luna A. P. Villanueva F. Dabad M. Esteve-Codina A. et al . (2019). Integrated single-base resolution maps of transcriptome, sRNAome and methylome of Tomato yellow leaf curl virus (TYLCV) in tomato. Sci. Rep.9, 2863. doi: 10.1038/s41598-019-39239-6

  • 112

    Popovsky-Sarid S. Borovsky Y. Faigenboim A. Parsons E. P. Lohrey G. T. Alkalai-Tuvia S. et al . (2017). Genetic and biochemical analysis reveals linked QTLs determining natural variation for fruit post-harvest water loss in pepper (Capsicum). Theor. Appl. Genet.130, 445459. doi: 10.1007/s00122-016-2825-9

  • 113

    Qiao L. Gao M. Wang Y. Tian X. Lu L. Liu X. (2022). Integrated transcriptomic and metabolomic analysis of cultivar differences provides insights into the browning mechanism of fresh-cut potato tubers. Postharvest Biol. Technol.188, 111905. doi: 10.1016/j.postharvbio.2022.111905

  • 114

    Qin C. Ahanger M. A. Lin B. Huang Z. Zhou J. Ahmed N. et al . (2021). Comparative transcriptome analysis reveals the regulatory effects of acetylcholine on salt tolerance of Nicotiana benthamiana. Phytochemistry181, 112582. doi: 10.1016/j.phytochem.2020.112582

  • 115

    Rai V. P. Rai A. Kumar R. Kumar S. Kumar S. Singh M. et al . (2016). Microarray analyses for identifying genes conferring resistance to pepper leaf curl virus in chilli pepper (Capsicum spp.). Genom. Data9, 140142. doi: 10.1016/j.gdata.2016.08.002

  • 116

    Razo-Mendivil F. G. Hernandez-Godínez F. Hayano-Kanashiro C. Martínez O. (2021). Transcriptomic analysis of a wild and a cultivated varieties of Capsicum annuum over fruit development and ripening. PLoS One16, e0256319. doi: 10.1371/journal.pone.0256319

  • 117

    Riccini A. Picarella M. E. De Angelis F. Mazzucato A. (2021). Bulk RNA-Seq analysis to dissect the regulation of stigma position in tomato. Plant Mol. Biol.105, 263285. doi: 10.1007/s11103-020-01086-9

  • 118

    Ruiu F. Picarella M. E. Imanishi S. Mazzucato A. (2015). A transcriptomic approach to identify regulatory genes involved in fruit set of wild-type and parthenocarpic tomato genotypes. Plant Mol. Biol.89, 263278. doi: 10.1007/s11103-015-0367-1

  • 119

    Sandhu K. S. Shiv A. Kaur G. Meena M. R. Raja A. K. Vengavasi K. et al . (2022). Integrated approach in genomic selection to accelerate genetic gain in sugarcane. Plants11, 2139. doi: 10.3390/plants11162139

  • 120

    Seyfferth C. Renema J. Wendrich J. R. Eekhout T. Seurinck R. Vandamme N. et al . (2021). Advances and opportunities in single-cell transcriptomics for plant research. Annu. Rev. Plant Biol. 72, 847866. doi: 10.1146/annurev-arplant-081720-010120

  • 121

    Shukla N. Yadav R. Kaur P. Rasmussen S. Goel S. Agarwal M. et al . (2018). Transcriptome analysis of root-knot nematode (Meloidogyne incognita)-infected tomato (Solanum lycopersicum) roots reveals complex gene expression profiles and metabolic networks of both host and nematode during susceptible and resistance responses. Mol. Plant Pathol.19, 615633. doi: 10.1111/mpp.12547

  • 122

    Singh J. Aggarwal R. Bashyal B. M. Darshan K. Meena B. R. Yadav J. et al . (2022). Temporal transcriptome of tomato elucidates the signaling pathways of induced systemic resistance and systemic acquired resistance activated by Chaetomium globosum. Front. Genet.13. doi: 10.3389/fgene.2022.1048578

  • 123

    Singh J. Aggarwal R. Bashyal B. M. Darshan K. Parmar P. Saharan M. S. et al . (2021a). Transcriptome reprogramming of tomato orchestrate the hormone signaling network of systemic resistance induced by Chaetomium globosum. Front. Plant Sci.12. doi: 10.3389/fpls.2021.721193

  • 124

    Singh S. Kaur I. Kariyat R. (2021b). The multifunctional roles of polyphenols in plant-herbivore interactions. Int. J. Mol. Sci.22, 1442. doi: 10.3390/ijms22031442

  • 125

    Singh J. Merchant A. Mayor L. Mbaye M. Gho C. Cooper M. et al . (2025). Understanding rates of genetic gain in sorghum [Sorghum bicolor (L.) Moench] in the United States. Plant Genome18, e70122. doi: 10.1002/tpg2.70122

  • 126

    Soyk S. Lemmon Z. H. Oved M. Fisher J. Liberatore K. L. Park S. J. et al . (2017). Bypassing negative epistasis on yield in tomato imposed by a domestication gene. Cell169, 11421155.e12. doi: 10.1016/j.cell.2017.04.032

  • 127

    Spanò R. Ferrara M. Gallitelli D. Mascia T. (2020). The role of grafting in the resistance of tomato to viruses. Plants9, 1042. doi: 10.3390/plants9081042

  • 128

    Stare T. Stare K. Weckwerth W. Wienkoop S. Gruden K. (2017). Comparison between proteome and transcriptome response in potato (Solanum tuberosum L.) leaves following potato virus Y (PVY) infection. Proteomes5, 14. doi: 10.3390/proteomes5030014

  • 129

    Stuart T. Satija R. (2019). Integrative single-cell analysis. Nat. Rev. Genet.20, 257272. doi: 10.1038/s41576-019-0093-7

  • 130

    Tan G. Liu K. Kang J. Xu K. Zhang Y. Hu L. et al . (2015). Transcriptome analysis of the compatible interaction of tomato with Verticillium dahliae using RNA-sequencing. Front. Plant Sci.6. doi: 10.3389/fpls.2015.00428

  • 131

    Tang R. Gupta S. K. Niu S. Li X.-Q. Yang Q. Chen G. et al . (2020). Transcriptome analysis of heat stress response genes in potato leaves. Mol. Biol. Rep.47, 43114321. doi: 10.1007/s11033-020-05485-5

  • 132

    Tang Z. Yu J. Xie J. Lyu J. Feng Z. Dawuda M. et al . (2019). Physiological and Growth Response of Pepper (Capsicum annum L.) Seedlings to Supplementary Red/Blue Light Revealed through Transcriptomic Analysis. Agronomy9, 139. doi: 10.3390/agronomy9030139

  • 133

    Tang R. Zhu W. Song X. Lin X. Cai J. Wang M. et al . (2016). Genome-wide identification and function analyses of heat shock transcription factors in potato. Front. Plant Sci.7. doi: 10.3389/fpls.2016.00490

  • 134

    Tengkun N. Dongdong W. Xiaohui M. Yue C. Qin C. (2019). Analysis of key genes involved in potato anthocyanin biosynthesis based on genomics and transcriptomics data. Front. Plant Sci.10. doi: 10.3389/fpls.2019.00603

  • 135

    Tewksbury J. J. Nabhan G. P. (2001). Seed dispersal. Directed deterrence by capsaicin in chilies. Nature412, 403404. doi: 10.1038/35086653

  • 136

    Tieman D. Zeigler M. Schmelz E. Taylor M. G. Rushing S. Jones J. B. et al . (2010). Functional analysis of a tomato salicylic acid methyl transferase and its role in synthesis of the flavor volatile methyl salicylate. Plant J.62, 113123. doi: 10.1111/j.1365-313X.2010.04128.x

  • 137

    Tolosa L. N. Zhang Z. (2020). The role of major transcription factors in solanaceous food crops under different stress conditions: current and future perspectives. Plants9, 56. doi: 10.3390/plants9010056

  • 138

    Upadhyay S. Jeena G. S. Kumar S. Shukla R. K. (2020). Asparagus racemosus bZIP transcription factor-regulated squalene epoxidase (ArSQE) promotes germination and abiotic stress tolerance in transgenic tobacco. Plant Sci.290, 110291. doi: 10.1016/j.plantsci.2019.110291

  • 139

    Usharani T. R. Ganachari M. Sowmya H. D. (2022). “ Omics in vegetable crops under the family Solanaceae,” in Omics in horticultural crops (London: Academic Press), 303332. doi: 10.1016/B978-0-323-89905-5.00022-7

  • 140

    Velez-Ramirez A. I. van Ieperen W. Vreugdenhil D. van Poppel P. M. J. A. Heuvelink E. Millenaar F. F. (2014). A single locus confers tolerance to continuous light and allows substantial yield increase in tomato. Nat. Commun.5, 4549. doi: 10.1038/ncomms5549

  • 141

    Wahrenburg Z. Benesch E. Lowe C. Jimenez J. Vulavala V. K. R. S. et al . (2021). Transcriptional regulation of wound suberin deposition in potato cultivars with differential wound healing capacity. Plant J.107, 7799. doi: 10.1111/tpj.15275

  • 142

    Wang X. Gao L. Jiao C. Stravoravdis S. Hosmani P. S. Saha S. et al . (2020). Genome of Solanum pimpinellifolium provides insights into structural variants during tomato breeding. Nat. Commun.11, 5817. doi: 10.1038/s41467-020-19682-0

  • 143

    Wang Z. Gerstein M. Snyder M. (2009). RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet.10, 5763. doi: 10.1038/nrg2484

  • 144

    Wang C. Hawthorne D. Qin Y. Pan X. Li Z. Zhu S. (2017). Impact of climate and host availability on future distribution of Colorado potato beetle. Sci. Rep.7, 4489. doi: 10.1038/s41598-017-04607-7

  • 145

    Wang F. Yin Y. Yu C. Li N. Shen S. Liu Y. et al . (2021). Transcriptomics analysis of heat stress-induced genes in pepper (Capsicum annuum L.) seedlings. Horticulturae7, 339. doi: 10.3390/horticulturae7100339

  • 146

    Wu H. Li H. Zhang W. Tang H. Yang L. (2021). Transcriptional regulation and functional analysis of Nicotiana tabacum under salt and ABA stress. Biochem. Biophys. Res. Commun.570, 110116. doi: 10.1016/j.bbrc.2021.07.011

  • 147

    Yamada K. Lim J. Dale J. M. Chen H. Shinn P. Palm C. J. et al . (2003). Empirical analysis of transcriptional activity in the Arabidopsis genome. Science302, 842846. doi: 10.1126/science.1088305

  • 148

    Yang X. Chen L. Yang Y. Guo X. Chen G. Xiong X. et al . (2020a). Transcriptome analysis reveals that exogenous ethylene activates immune and defense responses in a high late blight resistant potato genotype. Sci. Rep.10, 21294. doi: 10.1038/s41598-020-78027-5

  • 149

    Yang Y. Liu J. Zhou X. Liu S. Zhuang Y. (2020b). Identification of WRKY gene family and characterization of cold stress-responsive WRKY genes in eggplant. PeerJ8, e8777. doi: 10.7717/peerj.8777

  • 150

    Yang Z. Yang Z. Yang C. Wang Z. Chen D. Xie Y. et al . (2020c). Identification and genetic analysis of alternative splicing of long non-coding RNAs in tomato initial flowering stage. Genomics112, 897907. doi: 10.1016/j.ygeno.2019.06.005

  • 151

    Ye J. Wang X. Hu T. Zhang F. Wang B. Li C. et al . (2017). An InDel in the promoter of Al-ACTIVATED MALATE TRANSPORTER9 selected during tomato domestication determines fruit malate contents and aluminum tolerance. Plant Cell29, 22492268. doi: 10.1105/tpc.17.00211

  • 152

    Yu Q. Liu S. Yu L. Xiao Y. Zhang S. Wang X. et al . (2021). RNA demethylation increases the yield and biomass of rice and potato plants in field trials. Nat. Biotechnol.39, 15811588. doi: 10.1038/s41587-021-00982-9

  • 153

    Yuste-Lisbona F. J. Fernández-Lozano A. Pineda B. Bretones S. Ortíz-Atienza A. García-Sogo B. et al . (2020). ENO regulates tomato fruit size through the floral meristem development network. Proc. Natl. Acad. Sci. U.S.A.117, 81878195. doi: 10.1073/pnas.1913688117

  • 154

    Zhang J. Bai R. Shang M. Xu X. Song H. Zheng S. et al . (2021a). Identification of pepper genes involved in the response to CO2 enrichment using RNA-Seq analysis. Hortic. Environ. Biotechnol.62, 113123. doi: 10.1007/s13580-020-00300-3

  • 155

    Zhang Y. Chao J. Li X. Zhang C. Khan R. Du S. et al . (2021c). Comparative transcriptome combined with biochemical and physiological analyses provide new insights toward cadmium accumulation with two contrasting Nicotiana species. Physiol. Plant173, 369383. doi: 10.1111/ppl.13431

  • 156

    Zhang J. Liu S. Zhu X. Chang Y. Wang C. Ma N. et al . (2023). A comprehensive evaluation of tomato fruit quality and identification of volatile compounds. Plants12, 2947. doi: 10.3390/plants12162947

  • 157

    Zhang B. Tieman D. M. Jiao C. Xu Y. Chen K. Fei Z. et al . (2016). Chilling-induced tomato flavor loss is associated with altered volatile synthesis and transient changes in DNA methylation. Proc. Natl. Acad. Sci. U.S.A.113, 1258012585. doi: 10.1073/pnas.1613910113

  • 158

    Zhang J. Wang Y. Zhao Y. Zhang Y. Zhang J. Ma H. et al . (2020). Transcriptome analysis reveals Nitrogen deficiency induced alterations in leaf and root of three cultivars of potato (Solanum tuberosum L.). PLoS One15, e0240662. doi: 10.1371/journal.pone.0240662

  • 159

    Zhang M. Zhang H. Tan J. Huang S. Chen X. Jiang D. et al . (2021b). Transcriptome analysis of eggplant root in response to root-knot nematode infection. Pathogens10, 470. doi: 10.3390/pathogens10040470

  • 160

    Zhao S. Fung-Leung W.-P. Bittner A. Ngo K. Liu X. (2014). Comparison of RNA-Seq and microarray in transcriptome profiling of activated T cells. PloS One9, e78644. doi: 10.1371/journal.pone.0078644

  • 161

    Zhao T. Liu W. Zhao Z. Yang H. Bao Y. Zhang D. et al . (2019). Transcriptome profiling reveals the response process of tomato carrying Cf-19 and Cladosporium fulvum interaction. BMC Plant Biol.19, 572. doi: 10.1186/s12870-019-2150-y

  • 162

    Zhou X. Bao S. Liu J. Zhuang Y. (2016). De Novo sequencing and analysis of the transcriptome of the wild eggplant species Solanum Aculeatissimum in response to Verticillium dahliae. Plant Mol. Biol. Rep.34, 11931203. doi: 10.1007/s11105-016-0998-7

  • 163

    Zhou P. Khan R. Li Q. Liu G. Xu N. Yang Y. et al . (2020). Transcriptomic analyses of chilling stress responsiveness in leaves of tobacco (Nicotiana tabacum) Seedlings. Plant Mol. Biol. Rep.38, 113. doi: 10.1007/s11105-019-01167-0

  • 164

    Zhou R. Kong L. Yu X. Ottosen C.-O. Zhao T. Jiang F. et al . (2019). Oxidative damage and antioxidant mechanism in tomatoes responding to drought and heat stress. Acta Physiol. Plant41, 20. doi: 10.1007/s11738-019-2805-1

  • 165

    Zhu F. Jadhav S. S. Tohge T. Salem M. A. Li J. M. Giovannoni J. J. et al . (2022). A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator. Plant Physiol. 190, 250266. doi: 10.1093/plphys/kiac200

  • 166

    Zhu G. Wang S. Huang Z. Zhang S. Liao Q. Zhang C. et al . (2018). Rewiring of the fruit metabolome in tomato breeding. Cell172, 249261.e12. doi: 10.1016/j.cell.2017.12.019

  • 167

    Zhu B. Yang Y. Li R. Fu D. Wen L. Luo Y. et al . (2015). RNA sequencing and functional analysis implicate the regulatory role of long non-coding RNAs in tomato fruit ripening. J. Exp. Bot.66, 44834495. doi: 10.1093/jxb/erv203

  • 168

    Zou J. Chen J. Tang N. Gao Y. Hong M. Wei W. et al . (2018). Transcriptome analysis of aroma volatile metabolism change in tomato (Solanum lycopersicum) fruit under different storage temperatures and 1-MCP treatment. Postharvest Biol. Technol.135, 5767. doi: 10.1016/j.postharvbio.2017.08.017

  • 169

    Zouine M. Fu Y. Chateigner-Boutin A.-L. Mila I. Frasse P. Wang H. et al . (2014). Characterization of the tomato ARF gene family uncovers a multi-levels post-transcriptional regulation including alternative splicing. PLoS One9, e84203. doi: 10.1371/journal.pone.0084203

  • 170

    Zrenner R. Genzel F. Verwaaijen B. Wibberg D. Grosch R. (2020). Necrotrophic lifestyle of Rhizoctonia solani AG3-PT during interaction with its host plant potato as revealed by transcriptome analysis. Sci. Rep.10, 12574. doi: 10.1038/s41598-020-68728-2

  • 171

    Zuo J. Grierson D. Courtney L. T. Wang Y. Gao L. Zhao X. et al . (2020). Relationships between genome methylation, levels of non-coding RNAs, mRNAs and metabolites in ripening tomato fruit. Plant J.103, 980994. doi: 10.1111/tpj.14778

Summary

Keywords

crop quality, gene regulation, microarray, RNA-seq, solanaceous, transcriptomics

Citation

Qiao X, Jian B, Qiu J, Singh J, Kaur I, Hao R, Singh J, Singh L, Zhang H, Yin X, Li X, Kaur G and Li H (2026) Transcriptomics in solanaceous crop improvement: advances and opportunities. Front. Plant Sci. 16:1750317. doi: 10.3389/fpls.2025.1750317

Received

20 November 2025

Revised

11 December 2025

Accepted

22 December 2025

Published

16 January 2026

Volume

16 - 2025

Edited by

José Pinela, Instituto Nacional de Investigação Agrária e Veterinária (INIAV), Portugal

Reviewed by

Xiangjia Min, Youngstown State University, United States

Hemasundar Alavilli, Sejong University, Republic of Korea

Shan Li, Chinese Academy of Sciences (CAS), China

Updates

Copyright

*Correspondence: Xiang Li, ; Gurleen Kaur, ; Haisu Li,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics