De novo Taproot Transcriptome Sequencing and Analysis of Major Genes Involved in Sucrose Metabolism in Radish (Raphanus sativus L.)

Integrative omics studies revealed synergistic link between sucrose metabolic isogenes and carbohydrates in poplar roots infected by Fusarium wilt

Advances in carbohydrate metabolism prompted its essential role in defense priming and sweet immunity during plant-pathogen interactions. Nevertheless, upstream responding enzymes in the sucrose metabolic pathway and associated carbohydrate derivatives underlying fungal pathogen challenges remain to be deciphered in Populus, a model tree species. In silico deduction of genomic features, including phylogenies, exon/intron distributions, cis-regulatory elements, and chromosomal localization, identified 59 enzyme genes (11 families) in the Populus genome. Spatiotemporal expression of the transcriptome and the quantitative real-time PCR revealed a minuscule number of isogenes that were predominantly expressed in roots. Upon the pathogenic Fusarium solani (Fs) exposure, dynamic changes in the transcriptomics atlas and experimental evaluation verified Susy (PtSusy2 and 3), CWI (PtCWI3), VI (PtVI2), HK (PtHK6), FK (PtFK6), and UGPase (PtUGP2) families, displaying promotions in their expressions at 48 and 72 h of post-inoculation (hpi). Using the gas chromatography-mass spectrometry (GC-MS)-based non-targeted metabolomics combined with a high-performance ion chromatography system (HPICS), approximately 307 metabolites (13 categories) were annotated that led to the quantification of 46 carbohydrates, showing marked changes between three compared groups. By contrast, some sugars (e.g., sorbitol, L-arabitol, trehalose, and galacturonic acid) exhibited a higher accumulation at 72 hpi than 0 hpi, while levels of α-lactose and glucose decreased, facilitating them as potential signaling molecules. The systematic overview of multi-omics approaches to dissect the effects of Fs infection provides theoretical cues for understanding defense immunity depending on fine-tuned Suc metabolic gene clusters and synergistically linked carbohydrate pools in trees.

Comparative Root Transcriptome Profiling and Gene Regulatory Network Analysis between Eastern and Western Carrot (Daucus carota L.) Cultivars Reveals Candidate Genes for Vascular Tissue Patterning

Carrot (Daucus carota L.) is a highly consumed vegetable rich in carotenoids, known for their potent antioxidant, anti-inflammatory, and immune-protecting properties. While genetic and molecular studies have largely focused on wild and Western carrot cultivars (cvs), little is known about the evolutionary interactions between closely related Eastern and Western cvs. In this study, we conducted comparative transcriptome profiling of root tissues from Eastern (UHSBC-23-1) and Western (UHSBC-100) carrot cv. to better understand differentially expressed genes (DEGs) associated with storage root development and vascular cambium (VC) tissue patterning. Through reference-guided TopHat mapping, we achieved an average mapping rate of 73.87% and identified a total of 3544 DEGs (p < 0.05). Functional annotation and gene ontology classification revealed 97 functional categories, including 33 biological processes, 19 cellular components, 45 metabolic processes, and 26 KEGG pathways. Notably, Eastern cv. exhibited enrichment in cell wall, plant-pathogen interaction, and signal transduction terms, while Western cv. showed dominance in photosynthesis, metabolic process, and carbon metabolism terms. Moreover, constructed gene regulatory network (GRN) for both cvs. obtained orthologs with 1222 VC-responsive genes of Arabidopsis thaliana. In Western cv, GRN revealed VC-responsive gene clusters primarily associated with photosynthetic processes and carbon metabolism. In contrast, Eastern cv. exhibited a higher number of stress-responsive genes, and transcription factors (e.g., MYB15, WRKY46, AP2/ERF TF connected via signaling pathways with NAC036) were identified as master regulators of xylem vessel differentiation and secondary cell wall thickening. By elucidating the comparative transcriptome profiles of Eastern and Western cvs. for the first time, our study provides valuable insights into the differentially expressed genes involved in root development and VC tissue patterning. The identification of key regulatory genes and their roles in these processes represents a significant advancement in our understanding of the evolutionary relations and molecular mechanisms underlying secondary growth of carrot and regulation by vascular cambium.

Insights into the regulation of energy metabolism during the seed-to-seedling transition in marine angiosperm Zostera marina L.: Integrated metabolomic and transcriptomic analysis

Seed development is a crucial phase in the life cycle of seed-propagated plants. As the only group of angiosperms that evolved from terrestrial plants to complete their life cycle submerged in marine environments, the mechanisms underlying seed development in seagrasses are still largely unknown. In the present study, we attempted to combine transcriptomic, metabolomic, and physiological data to comprehensively analyze the molecular mechanism that regulates energy metabolism in Zostera marina seeds at the four major developmental stages. Our results demonstrated that seed metabolism was reprogrammed with significant alteration of starch and sucrose metabolism, glycolysis, the tricarboxylic acid cycle (TCA cycle), and the pentose phosphate pathway during the transition from seed formation to seedling establishment. The interconversion of starch and sugar provided energy storage substances in mature seeds and further acted as energy sources to support seed germination and seedling growth. The glycolysis pathway was active during Z. marina germination and seedling establishment, which provided pyruvate for TCA cycle by decomposing soluble sugar. Notably, the biological processes of glycolysis were severely inhibited during Z. marina seed maturation may have a positive effect on seed germination, maintaining a low level of metabolic activity during seed maturation to preserve seed viability. Increased acetyl-CoA and ATP contents were accompanied with the higher TCA cycle activity during seed germination and seedling establishment, indicating that the accumulations of precursor and intermediates metabolite that can strengthen the TCA cycle and facilitate energy supply for Z. marina seed germination and seedling growth. The large amount of oxidatively generated sugar phosphate promotes fructose 1,6-bisphosphate synthesis to feed back to glycolysis during seed germination, indicating that the pentose phosphate pathway not only provides energy for germination, but also complements the glycolytic pathway. Collectively, our findings suggest these energy metabolism pathways cooperate with each other in the process of seed transformation from maturity to seedling establishment, transforming seed from storage tissue to highly active metabolic tissue to meet the energy requirement seed development. These findings provide insights into the roles of the energy metabolism pathway in the complete developmental process of Z. marina seeds from different perspectives, which could facilitate habitat restoration of Z. marina meadows via seeds.

Metabonomic profiling of clubroot-susceptible and clubroot-resistant radish and the assessment of disease-resistant metabolites

Plasmodiophora brassicae causes a serious threat to cruciferous plants including radish (Raphanus sativus L.). Knowledge on the pathogenic regularity and molecular mechanism of P. brassicae and radish is limited, especially on the metabolism level. In the present study, clubroot-susceptible and clubroot-resistant cultivars were inoculated with P. brassicae Race 4, root hairs initial infection of resting spores (10⁷ CFU/mL) at 24 h post-inoculation and root galls symptom arising at cortex splitting stage were identified on both cultivars. Root samples of cortex splitting stage of two cultivars were collected and used for untargeted metabonomic analysis. We demonstrated changes in metabolite regulation and pathways during the cortex splitting stage of diseased roots between clubroot-susceptible and clubroot-resistant cultivars using untargeted metabonomic analysis. We identified a larger number of differentially regulated metabolites and heavier metabolite profile changes in the susceptible cultivar than in the resistant counterpart. The metabolites that were differentially regulated in both cultivars were mostly lipids and lipid-like molecules. Significantly regulated metabolites and pathways according to the P value and variable important in projection score were identified. Moreover, four compounds, including ethyl α-D-thioglucopyranoside, imipenem, ginsenoside Rg1, and 6-gingerol, were selected, and their anti-P. brassicae ability and effects on seedling growth were verified on the susceptible cultivar. Except for ethyl α-D-thioglucopyranoside, the remaining could inhibit clubroot development of varing degree. The use of 5 mg/L ginsenoside Rg1 + 5 mg/L 6-gingerol resulted in the lowest disease incidence and disease index among all treatments and enhanced seedling growth. The regulation of pathways or metabolites of carbapenem and ginsenoside was further explored. The results provide a preliminary understanding of the interaction between radish and P. brassicae at the metabolism level, as well as the development of measures for preventing clubroot.

Transcriptomic analysis of tuberous root in two sweet potato varieties reveals the important genes and regulatory pathways in tuberous root development

Background

Tuberous root formation and development is a complex process in sweet potato, which is regulated by multiple genes and environmental factors. However, the regulatory mechanism of tuberous root development is unclear.

Results

In this study, the transcriptome of fibrous roots (R0) and tuberous roots in three developmental stages (Rl, R2, R3) were analyzed in two sweet potato varieties, GJS-8 and XGH. A total of 22,914 and 24,446 differentially expressed genes (DEGs) were identified in GJS-8 and XGH respectively, 15,920 differential genes were shared by GJS-8 and XGH. KEGG pathway enrichment analysis showed that the DEGs shared by GJS-8 and XGH were mainly involved in “plant hormone signal transduction” “starch and sucrose metabolism” and “MAPK signal transduction”. Trihelix transcription factor (Tai6.25300) was found to be closely related to tuberous root enlargement by the comprehensive analysis of these DEGs and weighted gene co-expression network analysis (WGCNA).

Conclusion

A hypothetical model of genetic regulatory network for tuberous root development of sweet potato is proposed, which emphasizes that some specific signal transduction pathways like “plant hormone signal transduction” “Ca²⁺signal” “MAPK signal transduction” and metabolic processes including “starch and sucrose metabolism” and “cell cycle and cell wall metabolism” are related to tuberous root development in sweet potato. These results provide new insights into the molecular mechanism of tuberous root development in sweet potato.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08670-x.

Transcriptomics Integrated with Metabolomics Unveil Carotenoids Accumulation and Correlated Gene Regulation in White and Yellow-Fleshed Turnip (Brassica rapa ssp. rapa)

Turnip (Brassica rapa ssp. rapa) is considered to be a highly nutritious and health-promoting vegetable crop, whose flesh color can be divided into yellow and white. It is widely accepted that yellow-fleshed turnips have higher nutritional value. However, reports about flesh color formation is lacking. Here, the white-fleshed inbred line, W21, and yellow-fleshed inbred line, W25, were profiled from the swollen root of the turnip at three developmental periods to elucidate the yellow color formation. Transcriptomics integrated with metabolomics analysis showed that the PSY gene was the key gene affecting the carotenoids formation in W25. The coding sequence of BrrPSY-W25 was 1278 bp and that of BrrPSY-W21 was 1275 bp, and BrrPSY was more highly expressed in swollen roots in W25 than in W21. Transient transgenic tobacco leaf over-expressing BrrPSY-W and BrrPSY-Y showed higher transcript levels and carotenoids contents. Results revealed that yellow turnip formation is due to high expression of the PSY gene rather than mutations in the PSY gene, indicating that a post-transcriptional regulatory mechanism may affect carotenoids formation. Results obtained in this study will be helpful for explaining the carotenoids accumulation of turnips.

Analysis of Phenotypic Characteristics and Sucrose Metabolism in the Roots of Raphanus sativus L

The taproot of radish (Raphanus sativus L.) is an important sink organ; it is morphologically diverse and contains large amounts of secondary metabolites. Sucrose metabolism is believed to be important in the development of sink organs. We measured the amounts of glucose, fructose, and sucrose in the roots of sixty three radish accessions and analyzed the association between the sugar content and the root phenotype. Fructose content correlated with the root color and length characteristics, glucose was the most abundant sugar in the roots, and the sucrose content was very low, compared to that of the hexoses in most of the accessions. Expression analysis of the genes involved in sucrose metabolism, transportation, starch synthesis, and cell wall synthesis was performed through RNA sequencing. The genes encoding sucrose synthases (SUSY) and the enzymes involved in the synthesis of cellulose were highly expressed, indicating that SUSY is involved in cell wall synthesis in radish roots. The positive correlation coefficient (R) between the sucrose content and the expression of cell wall invertase and sugar transporter proteins suggest that hexose accumulation could occur through the apoplastic pathway in radish roots. A positive R score was also obtained when comparing the expression of genes encoding SUSY and fructokinase (FK), suggesting that the fructose produced by SUSY is mostly phosphorylated by FK. In addition, we concluded that sucrose was the most metabolized sugar in radish roots.

De novo Sequencing and Comparative Transcriptome Analyses Provide First Insights Into Polysaccharide Biosynthesis During Fruiting Body Development of Lentinula edodes

Polysaccharides separated from Lentinula edodes are well known for their medicinal properties. However, the precise molecular mechanisms of polysaccharide biosynthesis in L. edodes remain unclear. In this study, the fruiting bodies of L. edodes in four developmental stages with significant differences in polysaccharide yield were collected, and the characteristics of polysaccharides were studied. De novo sequencing and comparative transcriptomic analysis were performed by using high-throughput Illumina RNA-sequencing. KS1P30, KS2P30, KS3P30, and KS4P30 were obtained from the four developmental stages, respectively, by hot water extraction and 30% ethanol precipitation. These four polysaccharides had good immune activity in vitro; all of them were β-glucopyranose with a high molecular weight. Glucose was the main monosaccharide component of these polysaccharides. High-quality clean reads (57.88, 53.17, 53.28, and 47.56 million for different growth stages) and mapping ratios ranging from 84.75 to 90.11% were obtained. In total, 11,493 (96.56%) unigenes and 18,924 (97.46%) transcripts were successfully annotated in five public databases. The biosynthetic pathway and related genes of LEFP30 were mined. The molecular mechanism of LEFP30 yield change in the different developmental stages was predicted. The results provide some insights into the possible mechanisms involved in the biosynthetic pathway of this kind of polysaccharide in L. edodes fruiting bodies. They also indicate that candidate genes can be used as important resources for biotechnology and molecular breeding to regulate L. edodes fruiting body polysaccharide biosynthesis.

Downregulation of the photosynthetic machinery and carbon storage signaling pathways mediate La2O3 nanoparticle toxicity on radish taproot formation

The molecular and physiological mechanisms of how rare earth oxide nanoparticles (NPs) alter radish (Raphanus sativus L.) taproot formation and cracking were investigated in the present study. We compared plants that received suspensions of 10, 50, 100, 300 mg L^-1 of La₂O₃ NPs, 300 m L^-1 La₂O₃ bulk-particles (BPs), 0.8 m L^-1 La³⁺, or only water for six days during their tuber formation period. 100 and 300 mg L^-1 La₂O₃ NPs exposure decreased storage root biomass by 38% and 60%, respectively, and they both induced visible root cracking. Physiological analyses showed that La₂O₃ NPs exposure (>100 mg L^-1) significantly inhibited leaf net photosynthetic rate, cell wall pectin synthesis of both storage root epidermis and xylem parenchyma tissues, but increased the contents of cellulose and hemicellulose 1 in root epidermis cell walls. Moreover, transcriptome analysis further found that La₂O₃ NPs changed root cell wall structure by down-regulating core genes involved in cell wall pectin and IAA biosynthesis, which coincided with the observed La₂O₃ NPs-induced root cracking. Our results revealed the molecular mechanisms related to cell wall carbohydrate metabolism in response to NPs stress, providing a step forward for understanding the causes of NPs phytotoxicity on edible plant taproot formation and cracking.

Change in Sucrose Cleavage Pattern and Rapid Starch Accumulation Govern Lily Shoot-to-Bulblet Transition in vitro

In bulb crops, bulbing is a key progress in micropropagation and is the feature that most distinguishes bulbous crops from other plants. Generally, bulbing involves a shoot-to-bulblet transition; however, the underlying mechanism remains elusive. We explored this process by tracking the shoot-to-bulblet transition under different culture conditions. Rapid starch accumulation occurred at 15 days after transplanting (DAT) in the bulblet-inducing treatments as confirmed via histological observations and the significant elevation of starch synthesis related-gene transcription, including LohAGPS, LohAGPL, LohGBSS, LohSS, and LohSBE. However, for shoots that did not transition to bulblets and maintained the shoot status, much higher soluble sugars were detected. Interestingly, we observed a clear shift from invertase-catalyzed to sucrose synthase-catalyzed sucrose cleavage pattern based on the differential expression of LohCWIN and LohSuSy during the key transition stage (prior to and after bulbing at 0-15 DAT). Shoots that transitioned into bulblets showed significantly higher LohSuSy expression, especially LohSuSy4 expression, than shoots that did not transition. A symplastic phloem unloading pathway at the bulblet emergence stage (15 DAT) was verified via the 6(5)-carboxyfluorescein diacetate fluorescent tracer. We propose that starch is the fundamental compound in the shoot-to-bulblet transition and that starch synthesis is likely triggered by the switch from apoplastic to symplastic sucrose unloading, which may be related to sucrose depletion. Furthermore, this study is the first to provide a complete inventory of the genes involved in starch metabolism based on our transcriptome data. Two of these genes, LohAGPS1.2b and LohSSIIId, were verified by rapid amplification of cDNA ends cloning, and these data will provide additional support for Lilium research since whole genome is currently lacking.

Effects of heat stress on the biological Maillard reaction, oxidative stress, and occurrence of internal browning in Japanese radish (Raphanus sativus L.)

Internal browning (IB) of Japanese radish, in which the entire central part of the root changes from light tan to dark brown, is a physiological disorder caused by heat stress at the root maturation stage and represents significant economic losses to farmers. In this study, we cultivated three cultivars with different tolerance to IB disorder under high-temperature conditions, and examined the physiological factors involved in IB disorder differences among cultivars. There was a close correlation between the severity of IB disorder among cultivars and the size of root parenchymal cells. Significantly smaller cells were observed in the IB-resistant cultivars compared to the IB-susceptible one. Glucose concentration and invertase gene expression level tended to be higher in the normal roots of the IB-susceptible cultivar than in those of the resistant cultivars. Results of immunohistochemical analysis revealed that hexanoyl-lysine and acrolein, which are indicators of oxidative stress, in parenchymal cells were significantly elevated in the normal roots of the IB-susceptible cultivar than in those of the resistant ones. Similarly, the normal roots of the IB-susceptible cultivar showed elevated levels of Nε-(carboxymethyl)lysine, which is an indicator of the Maillard reaction, compared with those of the resistant cultivars. In addition, the immunohistochemistry of these indicators of oxidative stress and the Maillard reaction were more strongly observed in the damaged roots of the IB-susceptible cultivar. These results indicate that the occurrence of IB in Japanese radish is closely related to the increase in cell size, enhancement of the Maillard reaction by elevation in reducing sugar concentration, and increase in intracellular oxidative stress.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size.

RESULTS

A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis.

CONCLUSION

This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

Early-stage sugar beet taproot development is characterized by three distinct physiological phases

Despite the agronomic importance of sugar beet (Beta vulgaris L.), the early-stage development of its taproot has only been poorly investigated. Thus, the mechanisms that determine growth and sugar accumulation in sugar beet are largely unknown. In the presented study, a physiological characterization of early-stage sugar beet taproot development was conducted. Activities were analyzed for fourteen key enzymes of carbohydrate metabolism in developing taproots over the first 80 days after sowing. In addition, we performed in situ localizations of selected carbohydrate-metabolic enzyme activities, anatomical investigations, and quantifications of soluble carbohydrates, hexose phosphates, and phytohormones. Based on the accumulation dynamics of biomass and sucrose, as well as on anatomical parameters, the early phase of taproot development could be subdivided into three stages-prestorage, transition, secondary growth and sucrose accumulation stage-each of which was characterized by distinct metabolic and phytohormonal signatures. The enzyme activity signatures corresponding to these stages were also shown to be robustly reproducible in experiments conducted in two additional locations. The results from this physiological phenotyping approach contribute to the identification of the key regulators of sugar beet taproot development and open up new perspectives for sugar beet crop improvement concerning both physiological marker-based breeding and biotechnological approaches.

Comparative transcriptome and metabolome analyses provide new insights into the molecular mechanisms underlying taproot thickening in Panax notoginseng

BACKGROUND

Taproot thickening is a complex biological process that is dependent on the coordinated expression of genes controlled by both environmental and developmental factors. Panax notoginseng is an important Chinese medicinal herb that is characterized by an enlarged taproot as the main organ of saponin accumulation. However, the molecular mechanisms of taproot enlargement are poorly understood.

RESULTS

A total of 29,957 differentially expressed genes (DEGs) were identified during the thickening process in the taproots of P. notoginseng. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment revealed that DEGs associated with "plant hormone signal transduction," "starch and sucrose metabolism," and "phenylpropanoid biosynthesis" were predominantly enriched. Further analysis identified some critical genes (e.g., RNase-like major storage protein, DA1-related protein, and Starch branching enzyme I) and metabolites (e.g., sucrose, glucose, fructose, malate, and arginine) that potentially control taproot thickening. Several aspects including hormone crosstalk, transcriptional regulation, homeostatic regulation between sugar and starch, and cell wall metabolism, were identified as important for the thickening process in the taproot of P. notoginseng.

CONCLUSION

The results provide a molecular regulatory network of taproot thickening in P. notoginseng and facilitate the further characterization of the genes responsible for taproot formation in root medicinal plants or crops.

Differential Regulation of Anthocyanins in Green and Purple Turnips Revealed by Combined De Novo Transcriptome and Metabolome Analysis

Purple turnip Brassica rapa ssp. rapa is highly appreciated by consumers but the metabolites and molecular mechanisms underlying the root skin pigmentation remain open to study. Herein, we analyzed the anthocyanin composition in purple turnip (PT) and green turnip (GT) at five developmental stages. A total of 21 anthocyanins were detected and classified into the six major anthocynanin aglycones. Distinctly, PT contains 20 times higher levels of anthocyanins than GT, which explain the difference in the root skin pigmentation. We further sequenced the transcriptomes and analyzed the differentially expressed genes between the two turnips. We found that PT essentially diverts dihydroflavonols to the biosynthesis of anthocyanins over flavonols biosynthesis by strongly down-regulating one flavonol synthase gene, while strikingly up-regulating dihydroflavonol 4-reductase (DFR), anthocyanidin synthase and UDP-glucose: flavonoid-3-O-glucosyltransferase genes as compared to GT. Moreover, a nonsense mutation identified in the coding sequence of the DFR gene may lead to a nonfunctional protein, adding another hurdle to the accumulation of anthocyanin in GT. We also uncovered several key members of MYB, bHLH and WRKY families as the putative main drivers of transcriptional changes between the two turnips. Overall, this study provides new tools for modifying anthocyanin content and improving turnip nutritional quality.

The radish genome database (RadishGD): an integrated information resource for radish genomics

Radish (Raphanus sativus L.) is an important root vegetable crop in the family Brassicaceae, which provides diverse nutrients for human health and is closely related to the Brassica crop species. Recently, we sequenced and assembled the radish genome into nine chromosome pseudomolecules. In addition, we developed diverse genomic resources, including genetic maps, molecular markers, transcriptome, genome-wide methylation and variome data. In this study, we describe the radish genome database (RadishGD), including details of data sets that we generated and the web interface that allows access to these data. RadishGD comprises six major units that enable researchers and general users to search, browse and analyze the radish genomic data in an integrated manner. The Search unit provides gene structures and sequences for gene models through keyword or BLAST searches. The Genome browser displays graphic representations of gene models, mRNAs, repetitive sequences, genome-wide methylation and variomes among various genotypes. The Functional annotation unit offers gene ontology, plant ontology, pathway and gene family information for gene models. The Genetic map unit provides information about markers and their genetic locations using two types of genetic maps. The Expression unit presents transcriptional characteristics and methylation levels for each gene in 18 tissues. All sequence data incorporated into RadishGD can be downloaded from the Data resources unit. RadishGD will be continually updated to serve as a community resource for radish genomics and breeding research.

Comparative transcriptome profiling of genes and pathways involved in leaf-patterning of Clivia miniata var. variegata

Clivia miniata var. variegata (Cmvv) typically possesses yellow- and green-striped leaves. The striped plant not only has a high ornamental value but also be suitable for photosynthesis and chloroplast development research. Our previous study had revealed that yellow stripes (YSs) of Cmvv leaves contain chlorophyll-less ineffective chloroplasts. However, mechanism of Cmvv variegation is yet to be investigated. In the study, transcriptomes of both the YSs and green stripes (GSs) from single Cmvv leaves were compared using high-throughput sequencing. A total of 688 differential expression genes (DEGs) were identified based on biological replications. The qRT-PCR results indicated that transcriptome profiles accurately reflected global transcriptome differences between YSs and GSs. Subcellular localization analysis suggested that 56 DEG proteins were targeted to chloroplasts, and might be involved in anterograde signaling and leaf patterning. Moreover, the DEGs were mostly enriched in photosynthesis and plant-pathogen interaction KEGG pathways. Meanwhile, there should be coordination interaction between the two pathways. Seven of the eight DEGs involved in photosynthesis KEGG pathway were chloroplast-encoded genes and distributed among different cistrons of chloroplast DNA (cpDNA) large single copy regions (LSC) which are more prone to mutation. It was proposed that the YSs were caused by mutation(s) in cpDNA LSC. Thus, when the primary zygote of Cmvv was chimeric in LSC, leaf might be yellow- and green-striped. The study would give new insights into plant variegation and offers candidate genes to guide future research attempting to breed variegated plants.

Transcriptome sequencing and analysis during seed growth and development in Euryale ferox Salisb

BACKGROUND

Euryale ferox Salisb., an annual aquatic plant, is the only species in the genus Euryale in the Nymphaeaceae. Seeds of E. ferox are a nutritious food and also used in traditional Chinese medicine (Qian Shi in Mandarin). The molecular events that occurred during seed development in E. ferox have not yet been characterized. In this study, we performed transcriptomic analysis of four developmental stages (T1, T2, T3, and T4) in E. ferox seeds with three biological replicates per developmental stage to understand the physiological and biochemical processes during E. ferox seeds development.

RESULTS

313,844,425 clean reads were assembled into 160,107 transcripts and 85,006 unigenes with N50 lengths of 2052 bp and 1399 bp, respectively. The unigenes were annotated using five public databases (NR, COG, Swiss-Prot, KEGG, and GO). In the KEGG database, all of the unigenes were assigned to 127 pathways, of which phenylpropanoid biosynthesis was associated with the synthesis of secondary metabolites during E. ferox seed growth and development. Phenylalanine ammonia-lyase (PAL) as the first key enzyme catalyzed the conversion of phenylalanine to trans-cinnamic acid, then was related to the synthesis of flavonoids, lignins and alkaloid. The expression of PAL1 reached its peak at T3 stage, followed by a slight decrease at T4 stage. Cytochrome P450 (P450), encoded by CYP84A1 (which also called ferulate-5-hydroxylase (F5H) in Arabidopsis), was mainly involved in the biosynthesis of lignins.

CONCLUSIONS

Our study provides a transcriptomic analysis to better understand the morphological changes and the accumulation of medicinal components during E. ferox seed development. The increasing expression of PAL and P450 encoded genes in phenylpropanoid biosynthesis may promote the maturation of E. ferox seed including size, color, hardness and accumulation of medicinal components.

Transcript profiling of sucrose synthase genes involved in sucrose metabolism among four carrot (Daucus carota L.) cultivars reveals distinct patterns

BACKGROUND

Carrot which contains lots of nutrients has a large demand around the world. The soluble sugar content in fleshy root of carrot directly influences its taste and quality. Sucrose, as an important member of soluble sugar, is the main product of photosynthesis in higher plants and it plays pivotal roles in physiological processes including energy supply, signal transduction, transcriptional regulation, starch and cellulose synthesis, and stress tolerance. Sucrose synthase is a key enzyme involved in sucrose metabolism and is closely related to sucrose content. However, the molecular mechanism involved in sucrose metabolism in carrot has lagged behind.

RESULTS

Here, carrot roots of five developmental stages from four carrot cultivars were collected, and the contents of soluble sugar and sucrose in different stages and cultivars were surveyed. Three DcSus genes (DcSus1, DcSus2, and DcSus3), with lengths of 2427 bp, 2454 bp and 2628 bp, respectively, were identified and cloned in carrot. Phylogenetic analysis from the deduced amino acid sequences suggested that three DcSus were clustered into three distinct groups (SUSI, II and III). Results of enzymatic profiles demonstrated that the DcSus activities showed decrease trends during taproot development. Correlation analysis indicated that the DcSus activity showed negative correlation with soluble sugar content and strong negative correlation with sucrose concentration. Quantitative real-time PCR analysis showed that the expression profiles of the DcSus genes are significantly different in carrot tissues (root, leaf blade, and petiole), and the expression levels of the DcSus genes in the leaf blade were much higher than that in the root and petiole. The expression profiles of DcSus genes showed strong negative correlation with both sucrose content and soluble sugar content.

CONCLUSIONS

During carrot root development, the soluble sugar content and sucrose content showed increasing trends, while DcSus activities had persisting declinations, which may be due to the decreasing expression levels of genes encoding sucrose synthase. Our data demonstrate that synthesis of sucrose in carrot tissue is closely related with DcSus genes. The results from our study would not only provide effective insights of sucrose metabolism in carrot, but also are beneficial for biologists to improve carrot quality.

Comparative proteomic analysis provides insight into a complex regulatory network of taproot formation in radish (Raphanus sativus L.)

The fleshy taproot of radish is an important storage organ determining its yield and quality. Taproot thickening is a complex developmental process in radish. However, the molecular mechanisms governing this process remain unclear at the proteome level. In this study, a comparative proteomic analysis was performed to analyze the proteome changes at three developmental stages of taproot thickening using iTRAQ approach. In total, 1862 differentially expressed proteins (DEPs) were identified from 6342 high-confidence proteins, among which 256 up-regulated proteins displayed overlapped accumulation in S1 (pre-cortex splitting stage) vs. S2 (cortex splitting stage) and S1 vs. S3 (expanding stage) pairs, whereas 122 up-regulated proteins displayed overlapped accumulation in S1 vs. S3 and S2 vs. S3 pairs. Gene Ontology (GO) and pathway enrichment analysis showed that these DEPs were mainly involved in several processes such as "starch and sucrose metabolism", "plant hormone signal transduction", and "biosynthesis of secondary metabolites". A high concordance existed between iTRAQ and RT-qPCR at the mRNA expression levels. Furthermore, association analysis showed that 187, 181, and 96 DEPs were matched with their corresponding differentially expressed genes (DEGs) in S1 vs. S2, S1 vs. S3, and S2 vs. S3 comparison, respectively. Notably, several functional proteins including cell division cycle 5-like protein (CDC5), expansin B1 (EXPB1), and xyloglucan endotransglucosylase/hydrolase protein 24 (XTH24) were responsible for cell division and expansion during radish taproot thickening process. These results could facilitate a better understanding of the molecular mechanism underlying taproot thickening, and provide valuable information for the identification of critical genes/proteins responsible for taproot thickening in root vegetable crops.

Comparative transcriptomics uncovers alternative splicing and molecular marker development in radish (Raphanus sativus L.)

Background

Alternative splicing (AS) plays important roles in gene expression and proteome diversity. Single nucleotide polymorphism (SNP) and insertion/deletion (InDel) are abundant polymorphisms and co-dominant inheritance markers, which have been widely used in germplasm identification, genetic mapping and marker-assisted selection in plants. So far, however, little information is available on utilization of AS events and development of SNP and InDel markers from transcriptome in radish.

Results

In this study, three radish transcriptome datasets were collected and aligned to the reference radish genome. A total of 56,530 AS events were identified from three radish genotypes with intron retention (IR) being the most frequent AS type, which accounted for 59.4% of the total expressed genes in radish. In all, 22,412 SNPs and 9436 InDels were identified with an average frequency of 1 SNP/17.9 kb and 1 InDel/42.5 kb, respectively. A total of 43,680 potential SSRs were identified in 31,604 assembled unigenes with a density of 1 SSR/2.5 kb. The ratio of SNPs with nonsynonymous/synonymous mutations was 1.05:1. Moreover, 35 SNPs and 200 InDels were randomly selected and validated by Sanger sequencing, 83.9% of the SNPs and 70% of the InDels exhibited polymorphism among these three genotypes. In addition, the 15 SNPs and 125 InDels were found to be unevenly distributed on 9 linkage groups. Furthermore, 40 informative InDel markers were successfully used for the genetic diversity analysis on 32 radish accessions.

Conclusions

These results would not only provide new insights into transcriptome complexity and AS regulation, but also furnish large amount of molecular marker resources for germplasm identification, genetic mapping and further genetic improvement of radish in breeding programs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3874-4) contains supplementary material, which is available to authorized users.

Differential Gene Expression Analysis in Polygonum minus Leaf upon 24 h of Methyl Jasmonate Elicitation

Polygonum minus is an herbal plant that grows in Southeast Asian countries and traditionally used as medicine. This plant produces diverse secondary metabolites such as phenolic compounds and their derivatives, which are known to have roles in plant abiotic and biotic stress responses. Methyl jasmonate (MeJA) is a plant signaling molecule that triggers transcriptional reprogramming in secondary metabolism and activation of defense responses against many biotic and abiotic stresses. However, the effect of MeJA elicitation on the genome-wide expression profile in the leaf tissue of P. minus has not been well-studied due to the limited genetic information. Hence, we performed Illumina paired-end RNA-seq for de novo reconstruction of P. minus leaf transcriptome to identify differentially expressed genes (DEGs) in response to MeJA elicitation. A total of 182,111 unique transcripts (UTs) were obtained by de novo assembly of 191.57 million paired-end clean reads using Trinity analysis pipeline. A total of 2374 UTs were identified to be significantly up-/down-regulated 24 h after MeJA treatment. These UTs comprising many genes related to plant secondary metabolite biosynthesis, defense and stress responses. To validate our sequencing results, we analyzed the expression of 21 selected DEGs by quantitative real-time PCR and found a good correlation between the two analyses. The single time-point analysis in this work not only provides a useful genomic resource for P. minus but also gives insights on molecular mechanisms of stress responses in P. minus.

Comparative Transcriptome Analysis Reveal Candidate Genes Potentially Involved in Regulation of Primocane Apex Rooting in Raspberry (Rubus spp.)

Raspberries (Rubus spp.) exhibit a unique rooting process that is initiated from the stem apex of primocane, conferring an unusual asexual mode of reproduction to this plant. However, the full complement of genes involved in this process has not been identified. To this end, the present study analyzed the transcriptomes of the Rubus primocane and floricane stem apex at three developmental stages by Digital Gene Expression profiling to identify genes that regulate rooting. Sequencing and de novo assembly yielded 26.82 Gb of nucleotides and 59,173 unigenes; 498, 7,346, 4,110, 7,900, 9,397, and 4,776 differently expressed genes were identified in paired comparisons of SAF1 (floricane at developmental stage 1) vs. SAP1 (primocane at developmental stage 1), SAF2 vs. SAP2, SAF3 vs. SAP3, SAP1 vs. SAP2, SAP1 vs. SAP3, and SAP2 vs. SAP3, respectively. SAP1 maintains an extension growth pattern; SAP2 then exhibits growth arrest and vertical (downward) gravitropic deflection; and finally, short roots begin to form on the apex of SAP3. The Kyoto Encyclopedia of Genes and Genomes enrichment analysis of SAP1 vs. SAP2 revealed 12 pathways that were activated in response to shoot growth arrest and root differentiation, including circadian rhythm-plant (ko04712) and plant hormone signal transduction (ko04075). Our results indicate that genes related to circadian rhythm, ethylene and auxin signaling, shoot growth, and root development are potentially involved in the regulation of primocane apex rooting in Rubus. These findings provide a basis for elucidating the molecular mechanisms of primocane apex rooting in this economically valuable crop.

Comparative Transcriptomic Analysis Reveals Novel Insights into the Adaptive Response of Skeletonema costatum to Changing Ambient Phosphorus

Phosphorus (P) is a limiting macronutrient for diatom growth and productivity in the ocean. Much effort has been devoted to the physiological response of marine diatoms to ambient P change, however, the whole-genome molecular mechanisms are poorly understood. Here, we utilized RNA-Seq to compare the global gene expression patterns of a marine diatom Skeletonema costatum grown in inorganic P-replete, P-deficient, and inorganic- and organic-P resupplied conditions. In total 34,942 unique genes were assembled and 20.8% of them altered significantly in abundance under different P conditions. Genes encoding key enzymes/proteins involved in P utilization, nucleotide metabolism, photosynthesis, glycolysis, and cell cycle regulation were significantly up-regulated in P-deficient cells. Genes participating in circadian rhythm regulation, such as circadian clock associated 1, were also up-regulated in P-deficient cells. The response of S. costatum to ambient P deficiency shows several similarities to the well-described responses of other marine diatom species, but also has its unique features. S. costatum has evolved the ability to re-program its circadian clock and intracellular biological processes in response to ambient P deficiency. This study provides new insights into the adaptive mechanisms to ambient P deficiency in marine diatoms.

Transcriptome Profiling of Taproot Reveals Complex Regulatory Networks during Taproot Thickening in Radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is one of the most important vegetable crops worldwide. Taproot thickening represents a critical developmental period that determines yield and quality in radish life cycle. To isolate differentially expressed genes (DGEs) involved in radish taproot thickening process and explore the molecular mechanism underlying taproot development, three cDNA libraries from radish taproot collected at pre-cortex splitting stage (L1), cortex splitting stage (L2), and expanding stage (L3) were constructed and sequenced by RNA-Seq technology. More than seven million clean reads were obtained from the three libraries, from which 4,717,617 (L1, 65.35%), 4,809,588 (L2, 68.24%) and 4,973,745 (L3, 69.45%) reads were matched to the radish reference genes, respectively. A total of 85,939 transcripts were generated from three libraries, from which 10,450, 12,325, and 7392 differentially expressed transcripts (DETs) were detected in L1 vs. L2, L1 vs. L3, and L2 vs. L3 comparisons, respectively. Gene Ontology and pathway analysis showed that many DEGs, including EXPA9, Cyclin, CaM, Syntaxin, MADS-box, SAUR, and CalS were involved in cell events, cell wall modification, regulation of plant hormone levels, signal transduction and metabolisms, which may relate to taproot thickening. Furthermore, the integrated analysis of mRNA-miRNA revealed that 43 miRNAs and 92 genes formed 114 miRNA-target mRNA pairs were co-expressed, and three miRNA-target regulatory networks of taproot were constructed from different libraries. Finally, the expression patterns of 16 selected genes were confirmed using RT-qPCR analysis. A hypothetical model of genetic regulatory network associated with taproot thickening in radish was put forward. The taproot formation of radish is mainly attributed to cell differentiation, division and expansion, which are regulated and promoted by certain specific signal transduction pathways and metabolism processes. These results could provide new insights into the complex molecular mechanism underlying taproot thickening and facilitate genetic improvement of taproot in radish.

Comparative Transcriptome Analyses Reveal a Special Glucosinolate Metabolism Mechanism in Brassica alboglabra Sprouts

Brassica sprouts contain abundant phytochemicals, especially glucosinolates (GSs). Various methods have been used to enhance GS content in sprouts. However, the molecular basis of GS metabolism in sprouts remains an open question. Here we employed RNA-seq analysis to compare the transcriptomes of high-GS (JL-08) and low-GS (JL-09) Brassica alboglabra sprouts. Paired-end Illumina RNA-seq reads were generated and mapped to the Brassica oleracea reference genome. The differentially expressed genes were analyzed between JL-08 and JL-09. Among these, 1477 genes were up-regulated and 1239 down-regulated in JL-09 compared with JL-08. Enrichment analysis of these differentially expressed genes showed that the GS biosynthesis had the smallest enrichment factor and the highest Q-value of all metabolic pathways in Kyoto Encyclopedia of Genes and Genomes database, indicating the main metabolic difference between JL-08 and JL-09 is the GS biosynthetic pathway. Thirty-seven genes of the sequenced data were annotated as putatively involved in GS biosynthesis, degradation, and regulation, of which 11 were differentially expressed in JL-08 and JL-09. The expression level of GS degradation enzyme myrosinase in high-GS JL-08 was lower compared with low-GS JL-09. Surprisingly, in high-GS JL-08, the expression levels of GS biosynthesis genes were also lower than those in low-GS JL-09. As the GS contents in sprouts are determined by dynamic equilibrium of seed stored GS mobilization, de novo synthesis, degradation, and extra transport, the result of this study leads us to suggest that efforts to increase GS content should focus on either raising GS content in seeds or decreasing myrosinase activity, rather than improving the expression level of GS biosynthesis genes in sprouts.