Gene co-expression analysis of tomato seed maturation reveals tissue-specific regulatory networks and hubs associated with the acquisition of desiccation tolerance and seed vigour

20 years of the Bio-Analytic Resource for Plant Biology

The Bio-Analytic Resource for Plant Biology ('the BAR', at https://bar.utoronto.ca) is celebrating its 20th year in operation in 2025. The BAR encompasses and provides visualization tools for large 'omics data sets from plants. The BAR covers data from Arabidopsis, tomato, wheat, barley and 29 other plant species (with data for 2 others to be released soon). These data include nucleotide and protein sequence data, gene expression data, protein-protein and protein-DNA interactions, protein structures, subcellular localizations, and polymorphisms. The data are stored in more than 200 relational databases holding 186 GB of data and are presented to the researchers via web apps. These web apps provide data analysis and visualization tools. Some of the most popular tools are eFP ('electronic fluorescent pictograph') Browsers, ePlants and ThaleMine (an Arabidopsis-specific instance of InterMine). The BAR was designated a Global Core Biodata Resource in 2023. Like other GCBRs, the BAR has excellent operational stability, provides access without login requirement, and provides an API for researchers to be able to access BAR data programmatically. We present in this update a new overarching search tool called Gaia that permits easy access to all BAR data, powered by machine learning and artificial intelligence.

Effect of berry maturity stages on the germination and protein constituents of African nightshade (Solanum scabrum) seeds

African nightshade (Solanum scabrum) is a vegetable of great importance in several African countries. Production by seed is constrained by limited access to high quality seed, leaving farmers unable to meet the growing demand. The aim of this study was to investigate effects of berry maturity stages (mature green and purple) on germination and protein components of African nightshade seeds. Nine accessions were used for the germination assays and three selected for the proteomics. The seeds harvested at the purple ripe stage showed a higher germination percentage than seeds harvested at the mature green stage. Seeds of three accessions, which showed different germination responses depending on maturity, were subjected to gel-based proteome comparison followed by mass spectrometry. Analysis revealed a total of 563 spots, of which 19-23% were differentially abundant between the two harvest times. This is the first study to utilise proteomics in African nightshade. Four types of seed storage proteins were identified. Identified proteins were associated with seed maturation, indicating more advanced maturation stages in accessions with higher germination capacity. In addition, a reference proteome map of S. scabrum was created using GelMap platform ( www.gelmap.de ). Colour change of African nightshade berries seems to be a suitable indicator for seed maturity.

Sleeping but not defenceless: seed dormancy and protection

To ensure their vital role in disseminating the species, dormant seeds have developed adaptive strategies to protect themselves against pathogens and predators. This is orchestrated through the synthesis of an array of constitutive defences that are put in place in a developmentally regulated manner, which are the focus of this review. We summarize the defence activity and the nature of the molecules coming from the exudate of imbibing seeds that leak into their vicinity, also referred to as the spermosphere. As a second layer of protection, the dual role of the seed coat will be discussed; as a physical barrier and a multi-layered reservoir of defence compounds that are synthesized during seed development. Since imbibed dormant seeds can persist in the soil for extensive periods, we address the question of whether during this time a constitutively regulated defence programme is switched on to provide further protection, via the well-defined pathogenesis-related (PR) protein family. In addition, we review the hormonal and signalling pathways that might be involved in the interplay between dormancy and defence and point out questions that need further attention.

Domestication has altered gene expression and secondary metabolites in pea seed coat

The mature seed in legumes consists of an embryo and seed coat. In contrast to knowledge about the embryo, we know relatively little about the seed coat. We analyzed the gene expression during seed development using a panel of cultivated and wild pea genotypes. Gene co-expression analysis identified gene modules related to seed development, dormancy, and domestication. Oxidoreductase genes were found to be important components of developmental and domestication processes. Proteomic and metabolomic analysis revealed that domestication favored proteins involved in photosynthesis and protein metabolism at the expense of seed defense. Seed coats of wild peas were rich in cell wall-bound metabolites and the protective compounds predominated in their seed coats. Altogether, we have shown that domestication altered pea seed development and modified (mostly reduced) the transcripts along with the protein and metabolite composition of the seed coat, especially the content of the compounds involved in defense. We investigated dynamic profiles of selected identified phenolic and flavonoid metabolites across seed development. These compounds usually deteriorated the palatability and processing of the seeds. Our findings further provide resources to study secondary metabolism and strategies for improving the quality of legume seeds which comprise an important part of the human protein diet.

Exploring the gene expression network involved in the heat stress response of a thermotolerant tomato genotype

BACKGROUND

The increase in temperatures due to the current climate change dramatically affects crop cultivation, resulting in yield losses and altered fruit quality. Tomato is one of the most extensively grown and consumed horticultural products, and although it can withstand a wide range of climatic conditions, heat stress can affect plant growth and development specially on the reproductive stage, severely influencing the final yield. In the present work, the heat stress response mechanisms of one thermotolerant genotype (E42) were investigated by exploring its regulatory gene network. This was achieved through a promoter analysis based on the identification of the heat stress elements (HSEs) mapping in the promoters, combined with a gene co-expression network analysis aimed at identifying interactions among heat-related genes.

RESULTS

Results highlighted 82 genes presenting HSEs in the promoter and belonging to one of the 52 gene networks obtained by the GCN analysis; 61 of these also interact with heat shock factors (Hsfs). Finally, a list of 13 candidate genes including two Hsfs, nine heat shock proteins (Hsps) and two GDSL esterase/lipase (GELPs) were retrieved by focusing on those E42 genes exhibiting HSEs in the promoters, interacting with Hsfs and showing variants, compared to Heinz reference genome, with HIGH and/or MODERATE impact on the translated protein. Among these, the Gene Ontology annotation analysis evidenced that only LeHsp100 (Solyc02g088610) belongs to a network specifically involved in the response to heat stress.

CONCLUSIONS

As a whole, the combination of bioinformatic analyses carried out on genomic and trascriptomic data available for tomato, together with polymorphisms detected in HS-related genes of the thermotolerant E42 allowed to determine a subset of candidate genes involved in the HS response in tomato. This study provides a novel approach in the investigation of abiotic stress response mechanisms and further studies will be conducted to validate the role of the highlighted genes.

Arabidopsis HSFA9 Acts as a Regulator of Heat Response Gene Expression and the Acquisition of Thermotolerance and Seed Longevity

Heat-shock transcription factors (HSFs) are crucial for regulating plant responses to heat and various stresses, as well as for maintaining normal cellular functions and plant development. HSFA9 and HSFA2 are two of the Arabidopsis class A HSFs and their expressions are dramatically induced in response to heat shock (HS) stress among all 21 Arabidopsis HSFs. However, the detailed biological roles of their cooperation have not been fully characterized. In this study, we employed an integrated approach that combined bioinformatics, molecular genetics and computational analysis to identify and validate the molecular mechanism that controls seed longevity and thermotolerance in Arabidopsis. The acquisition of tolerance to deterioration was accompanied by a significant transcriptional switch that involved the induction of primary metabolism, reactive oxygen species and unfolded protein response, as well as the regulation of genes involved in response to dehydration, heat and hypoxia. In addition, the cis-regulatory motif analysis in normal stored and controlled deterioration treatment (CDT) seeds confirmed the CDT-repressed genes with heat-shock element (HSE) in their promoters. Using a yeast two-hybrid and molecular dynamic interaction assay, it is shown that HSFA9 acted as a potential regulator that can interact with HSFA2. Moreover, the knock-out mutants of both HSFA9 and HSFA2 displayed a significant reduction in seed longevity. These novel findings link HSF transcription factors with seed deterioration tolerance and longevity.

OsCSD2 and OsCSD3 Enhance Seed Storability by Modulating Antioxidant Enzymes and Abscisic Acid in Rice

Seed deterioration during storage poses a significant challenge to rice production, leading to a drastic decline in both edible quality and viability, thereby impacting overall crop yield. This study aimed to address this issue by further investigating candidate genes associated with two previously identified QTLs for seed storability through genome association analysis. Among the screened genes, two superoxide dismutase (SOD) genes, OsCSD2 (Copper/zinc Superoxide Dismutase 2) and OsCSD3, were selected for further study. The generation of overexpression and CRISPR/Cas9 mutant transgenic lines revealed that OsCSD2 and OsCSD3 play a positive regulatory role in enhancing rice seed storability. Subsequent exploration of the physiological mechanisms demonstrated that overexpression lines exhibited lower relative electrical conductivity, indicative of reduced cell membrane damage, while knockout lines displayed the opposite trend. Furthermore, the overexpression lines of OsCSD2 and OsCSD3 showed significant increases not only in SOD but also in CAT and POD activities, highlighting an augmented antioxidant system in the transgenic seeds. Additionally, hormone profiling indicated that ABA contributed to the improved seed storability observed in these lines. In summary, these findings provide valuable insights into the regulatory mechanisms of OsCSDs in rice storability, with potential applications for mitigating grain loss and enhancing global food security.

Identification of transcriptional modules linked to the drought response of Brassica napus during seed development and their mitigation by early biotic stress

In order to capture the drought impacts on seed quality acquisition in Brassica napus and its potential interaction with early biotic stress, seeds of the 'Express' genotype of oilseed rape were characterized from late embryogenesis to full maturity from plants submitted to reduced watering (WS) with or without pre-occurring inoculation by the telluric pathogen Plasmodiophora brassicae (Pb + WS or Pb, respectively), and compared to control conditions (C). Drought as a single constraint led to significantly lower accumulation of lipids, higher protein content and reduced longevity of the WS-treated seeds. In contrast, when water shortage was preceded by clubroot infection, these phenotypic differences were completely abolished despite the upregulation of the drought sensor RD20. A weighted gene co-expression network of seed development in oilseed rape was generated using 72 transcriptomes from developing seeds from the four treatments and identified 33 modules. Module 29 was highly enriched in heat shock proteins and chaperones that showed a stronger upregulation in Pb + WS compared to the WS condition, pointing to a possible priming effect by the early P. brassicae infection on seed quality acquisition. Module 13 was enriched with genes encoding 12S and 2S seed storage proteins, with the latter being strongly upregulated under WS conditions. Cis-element promotor enrichment identified PEI1/TZF6, FUS3 and bZIP68 as putative regulators significantly upregulated upon WS compared to Pb + WS. Our results provide a temporal co-expression atlas of seed development in oilseed rape and will serve as a resource to characterize the plant response towards combinations of biotic and abiotic stresses.

Seed Longevity and Ageing: A Review on Physiological and Genetic Factors with an Emphasis on Hormonal Regulation

Upon storage, seeds inevitably age and lose their viability over time, which determines their longevity. Longevity correlates with successful seed germination and enhancing this trait is of fundamental importance for long-term seed storage (germplasm conservation) and crop improvement. Seed longevity is governed by a complex interplay between genetic factors and environmental conditions experienced during seed development and after-ripening that will shape seed physiology. Several factors have been associated with seed ageing such as oxidative stress responses, DNA repair enzymes, and composition of seed layers. Phytohormones, mainly abscisic acid, auxins, and gibberellins, have also emerged as prominent endogenous regulators of seed longevity, and their study has provided new regulators of longevity. Gaining a thorough understanding of how hormonal signalling genes and pathways are integrated with downstream mechanisms related to seed longevity is essential for formulating strategies aimed at preserving seed quality and viability. A relevant aspect related to research in seed longevity is the existence of significant differences between results depending on the seed equilibrium relative humidity conditions used to study seed ageing. Hence, this review delves into the genetic, environmental and experimental factors affecting seed ageing and longevity, with a particular focus on their hormonal regulation. We also provide gene network models underlying hormone signalling aimed to help visualize their integration into seed longevity and ageing. We believe that the format used to present the information bolsters its value as a resource to support seed longevity research for seed conservation and crop improvement.

Genetic Variability in Seed Longevity and Germination Traits in a Tomato MAGIC Population in Contrasting Environments

The stable production of high vigorous seeds is pivotal to crop yield. Also, a high longevity is essential to avoid progressive loss of seed vigour during storage. Both seed traits are strongly influenced by the environment during seed development. Here, we investigated the impact of heat stress (HS) during fruit ripening on tomato seed lifespan during storage at moderate relative humidity, speed (t50) and homogeneity of germination, using a MAGIC population that was produced under optimal and HS conditions. A plasticity index was used to assess the extent of the impact of HS for each trait. HS reduced the average longevity and germination homogeneity by 50% within the parents and MAGIC population. However, there was a high genetic variability in the seed response to heat stress. A total of 39 QTLs were identified, including six longevity QTLs for seeds from control (3) and HS (3) conditions, and six plasticity QTLs for longevity, with only one overlapping with a longevity QTL under HS. Four out of the six longevity QTL co-located with t50 QTL, revealing hotspots for seed quality traits. Twenty-one QTLs with intervals below 3 cM were analyzed using previous transcriptome and gene network data to propose candidate genes for seed vigour and longevity traits.

Tissues and mechanisms associated with Verticillium wilt resistance in tomato using bi-grafted near-isogenic lines

Host resistance is the primary means to control Verticillium dahliae, a soil-borne pathogen causing major losses on a broad range of plants, including tomato. The tissues and mechanisms responsible for resistance remain obscure. In the field, resistant tomato used as rootstocks does not confer resistance. Here, we created bi-grafted plants with near-isogenic lines (NILs) exhibiting (Ve1) or lacking (ve1) resistance to V. dahliae race 1. Ten days after inoculation, scion and rootstock tissues were subjected to differential gene expression and co-expression network analyses. Symptoms only developed in susceptible scions regardless of the rootstock. Infection caused more dramatic alteration of tomato gene expression in susceptible compared with resistant tissues, including pathogen receptor, signaling pathway, pathogenesis-related protein, and cell wall modification genes. Differences were observed between scions and rootstocks, primarily related to physiological processes in these tissues. Gene expression in scions was influenced by the rootstock genotype. A few genes were associated with the Ve1 genotype, which was independent of infection or tissue type. Several were physically clustered, some near the Ve1 locus on chromosome 9. Transcripts mapped to V. dahliae were dominated by secreted candidate effector proteins. These findings advance knowledge of molecular mechanisms underlying the tomato-V. dahliae interaction.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Improvement of Seed Germination under Salt Stress via Overexpressing Caffeic Acid O-methyltransferase 1 (SlCOMT1) in Solanum lycopersicum L

Melatonin (MT) is a phytohormone-like substance and is profoundly involved in modulating nearly all aspects of plant development and acclimation to environmental stressors. However, there remain no studies about the effects of MT on tomato seed germination under salt stress. Here we reported that the overexpression of caffeic acid O-methyltransferase 1 (SlCOMT1) significantly increased both MT content and salt tolerance in the germinated seeds of a transgenic tomato relative to wild type (WT) samples. Physiological investigation showed higher amylase activity in the stressed overexpression seeds than WT, leading to the promoted starch decomposition and enhanced soluble sugar content. The stimulated production of osmolytes and enhanced activities of SOD, POD, and CAT, together with the significant reduction in H₂O₂ and O₂^·- accumulation, were revealed in the stressed overexpression seeds relative to WT, largely accounting for their lower membrane lipid peroxidation. qPCR assay showed that, upon salt stress, the transcript abundance of hub genes related to germination (SlCYP707A1, SlABA1, SlGA3ox2 and SlGA2ox4) and stress tolerance (SlCDPK1, SlWRKY33 and SlMAPK1) were distinctly altered in the overexpression samples when compared to WT, providing a molecular basis for MT-mediated improvement of seed salt tolerance. Altogether, our observations shed new insights into biological functions of SlCOMT1 and could expand its utilization in genetic improvement of tomato salt tolerance in future.

Chromatin dynamics associated with seed desiccation tolerance/sensitivity at early germination in Medicago truncatula

Desiccation tolerance (DT) has contributed greatly to the adaptation of land plants to severe water-deficient conditions. DT is mostly observed in reproductive parts in flowering plants such as seeds. The seed DT is lost at early post germination stage but is temporally re-inducible in 1 mm radicles during the so-called DT window following a PEG treatment before being permanently silenced in 5 mm radicles of germinating seeds. The molecular mechanisms that activate/reactivate/silence DT in developing and germinating seeds have not yet been elucidated. Here, we analyzed chromatin dynamics related to re-inducibility of DT before and after the DT window at early germination in Medicago truncatula radicles to determine if DT-associated genes were transcriptionally regulated at the chromatin levels. Comparative transcriptome analysis of these radicles identified 948 genes as DT re-induction-related genes, positively correlated with DT re-induction. ATAC-Seq analyses revealed that the chromatin state of genomic regions containing these genes was clearly modulated by PEG treatment and affected by growth stages with opened chromatin in 1 mm radicles with PEG (R1P); intermediate openness in 1 mm radicles without PEG (R1); and condensed chromatin in 5 mm radicles without PEG (R5). In contrast, we also showed that the 103 genes negatively correlated with the re-induction of DT did not show any transcriptional regulation at the chromatin level. Additionally, ChIP-Seq analyses for repressive marks H2AK119ub and H3K27me3 detected a prominent signal of H3K27me3 on the DT re-induction-related gene sequences at R5 but not in R1 and R1P. Moreover, no clear H2AK119ub marks was observed on the DT re-induction-related gene sequences at both developmental radicle stages, suggesting that silencing of DT process after germination will be mainly due to H3K27me3 marks by the action of the PRC2 complex, without involvement of PRC1 complex. The dynamic of chromatin changes associated with H3K27me3 were also confirmed on seed-specific genes encoding potential DT-related proteins such as LEAs, oleosins and transcriptional factors. However, several transcriptional factors did not show a clear link between their decrease of chromatin openness and H3K27me3 levels, suggesting that their accessibility may also be regulated by additional factors, such as other histone modifications. Finally, in order to make these comprehensive genome-wide analyses of transcript and chromatin dynamics useful to the scientific community working on early germination and DT, we generated a dedicated genome browser containing all these data and publicly available at https://iris.angers.inrae.fr/mtseedepiatlas/jbrowse/?data=Mtruncatula.

Seed-to-Seedling Transition in Pisum sativum L.: A Transcriptomic Approach

The seed-to-seedling transition is a crucial step in the plant life cycle. The transition occurs at the end of seed germination and corresponds to the initiation of embryonic root growth. To improve our understanding of how a seed transforms into a seedling, we germinated the Pisum sativum L. seeds for 72 h and divided them into samples before and after radicle protrusion. Before radicle protrusion, seeds survived after drying and formed normally developed seedlings upon rehydration. Radicle protrusion increased the moisture content level in seed axes, and the accumulation of ROS first generated in the embryonic root and plumule. The water and oxidative status shift correlated with the desiccation tolerance loss. Then, we compared RNA sequencing-based transcriptomics in the embryonic axes isolated from pea seeds before and after radicle protrusion. We identified 24,184 differentially expressed genes during the transition to the post-germination stage. Among them, 2101 genes showed more prominent expression. They were related to primary and secondary metabolism, photosynthesis, biosynthesis of cell wall components, redox status, and responses to biotic stress. On the other hand, 415 genes showed significantly decreased expression, including the groups related to water deprivation (eight genes) and response to the ABA stimulus (fifteen genes). We assume that the water deprivation group, especially three genes also belonging to ABA stimulus (LTI65, LTP4, and HVA22E), may be crucial for the desiccation tolerance loss during a metabolic switch from seed to seedling. The latter is also accompanied by the suppression of ABA-related transcription factors ABI3, ABI4, and ABI5. Among them, HVA22E, ABI4, and ABI5 were highly conservative in functional domains and showed homologous sequences in different drought-tolerant species. These findings elaborate on the critical biochemical pathways and genes regulating seed-to-seedling transition.

Genetic Aspects and Molecular Causes of Seed Longevity in Plants-A Review

Seed longevity is the most important trait related to the management of gene banks because it governs the regeneration cycle of seeds. Thus, seed longevity is a quantitative trait. Prior to the discovery of molecular markers, classical genetic studies have been performed to identify the genetic determinants of this trait. Post-2000 saw the use of DNA-based molecular markers and modern biotechnological tools, including RNA sequence (RNA-seq) analysis, to understand the genetic factors determining seed longevity. This review summarizes the most important and relevant genetic studies performed in Arabidopsis (24 reports), rice (25 reports), barley (4 reports), wheat (9 reports), maize (8 reports), soybean (10 reports), tobacco (2 reports), lettuce (1 report) and tomato (3 reports), in chronological order, after discussing some classical studies. The major genes identified and their probable roles, where available, are debated in each case. We conclude by providing information about many different collections of various crops available worldwide for advanced research on seed longevity. Finally, the use of new emerging technologies, including RNA-seq, in seed longevity research is emphasized by providing relevant examples.

The Fruit Proteome Response to the Ripening Stages in Three Tomato Genotypes

The tomato is a horticultural crop that appears in various colors as it ripens. Differences in the proteome expression abundance of a tomato depend on its genotype and ripening stage. Thus, this study aimed to confirm the differences in changes in the proteome according to four ripening stages (green, breaker, turning, and mature) of three tomato genotypes, i.e., yellow, black, and red tomatoes, using a gel-based proteomic technique. The number of protein spots shown as two-dimensional electrophoresis (2-DE) gels differed according to tomato genotype and ripening stage. A total of 286 variant proteins were determined using matrix-assisted laser desorption-time of flight (MALDI-TOF) mass spectrometry (MS) analysis, confirming 233 identified protein functions. In three tomato genotypes in each ripening stage, grouping according to the Munich Information Center for Protein Sequences (MIPS) functional categories confirmed the variant proteins involved in the following: energy processes (21%); metabolism (20%); protein fate (15%); protein synthesis (10%); a protein with a binding function or cofactor requirement (8%); cell rescue, defense, and virulence (8%); cellular transport, transport facilitation, and transport routes (6%); the biogenesis of cellular components (5%); cell cycle and DNA processing (2%); others (5%). Among the identified protein spots in the function category, two proteins related to metabolism, four related to energy, four related to protein synthesis, and two related to interaction with the cellular environment showed significantly different changes according to the fruit color by the ripening stage. This study reveals the physiological changes in different types of tomatoes according to their ripening stage and provides information on the proteome for further improvement.

Seed germination and vigor: ensuring crop sustainability in a changing climate

In the coming decades, maintaining a steady food supply for the increasing world population will require high-yielding crop plants which can be productive under increasingly variable conditions. Maintaining high yields will require the successful and uniform establishment of plants in the field under altered environmental conditions. Seed vigor, a complex agronomic trait that includes seed longevity, germination speed, seedling growth, and early stress tolerance, determines the duration and success of this establishment period. Elevated temperature during early seed development can decrease seed size, number, and fertility, delay germination and reduce seed vigor in crops such as cereals, legumes, and vegetable crops. Heat stress in mature seeds can reduce seed vigor in crops such as lettuce, oat, and chickpea. Warming trends and increasing temperature variability can increase seed dormancy and reduce germination rates, especially in crops that require lower temperatures for germination and seedling establishment. To improve seed germination speed and success, much research has focused on selecting quality seeds for replanting, priming seeds before sowing, and breeding varieties with improved seed performance. Recent strides in understanding the genetic basis of variation in seed vigor have used genomics and transcriptomics to identify candidate genes for improving germination, and several studies have explored the potential impact of climate change on the percentage and timing of germination. In this review, we discuss these recent advances in the genetic underpinnings of seed performance as well as how climate change is expected to affect vigor in current varieties of staple, vegetable, and other crops.

Transcriptome profiles revealed molecular mechanisms of alternating temperatures in breaking the epicotyl morphophysiological dormancy of Polygonatum sibiricum seeds

BACKGROUND

To adapt seasonal climate changes under natural environments, Polygonatum sibiricum seeds have a long period of epicotyl morphophysiological dormancy, which limits their wide-utilization in the large-scale plant progeny propagation. It has been proven that the controlled consecutive warm and cold temperature treatments can effectively break and shorten this seed dormancy status to promote its successful underdeveloped embryo growth, radicle emergence and shoot emergence. To uncover the molecular basis of seed dormancy release and seedling establishment, a SMRT full-length sequencing analysis and an Illumina sequencing-based comparison of P. sibiricum seed transcriptomes were combined to investigate transcriptional changes during warm and cold stratifications.

RESULTS

A total of 87,251 unigenes, including 46,255 complete sequences, were obtained and 77,148 unigenes (88.42%) were annotated. Gene expression analyses at four stratification stages identified a total of 27,059 DEGs in six pairwise comparisons and revealed that more differentially expressed genes were altered at the Corm stage than at the other stages, especially Str_S and Eme. The expression of 475 hormone metabolism genes and 510 hormone signaling genes was modulated during P. sibiricum seed dormancy release and seedling emergence. One thousand eighteen transcription factors and five hundred nineteen transcription regulators were detected differentially expressed during stratification and germination especially at Corm and Str_S stages. Of 1246 seed dormancy/germination known DEGs, 378, 790, and 199 DEGs were associated with P. sibiricum MD release (Corm vs Seed), epicotyl dormancy release (Str_S vs Corm), and the seedling establishment after the MPD release (Eme vs Str_S).

CONCLUSIONS

A comparison with dormancy- and germination-related genes in Arabidopsis thaliana seeds revealed that genes related to multiple plant hormones, chromatin modifiers and remodelers, DNA methylation, mRNA degradation, endosperm weakening, and cell wall structures coordinately mediate P. sibiricum seed germination, epicotyl dormancy release, and seedling establishment. These results provided the first insights into molecular regulation of P. sibiricum seed epicotyl morphophysiological dormancy release and seedling emergence. They may form the foundation of future studies regarding gene interaction and the specific roles of individual tissues (endosperm, newly-formed corm) in P. sibiricum bulk seed dormancy.