Transcription Profile of Potato (Solanum tuberosum L.) Growing In Vitro

Despite countless papers on plant tissue culture, few have assessed the transcriptional changes that occur in a developing explant from the moment it is cut and plated, and during early stages of development. In this study, the mRNA profile of in vitro potato (Solanum tuberosum L.) plantlets derived from internode stems was assessed via four comparisons during five time intervals (0 h, 24 h, 48 h, 1 week, and 4 weeks) using the Kyoto Encyclopedia of Genes and Genomes (KEGG) classification. After screening KEGG metabolic groups, 35 processes were shown to be either up- or down-regulated. The vast majority (28 processes) were related to growth or development. The expression intensity of 40,430 genes was assessed and 158, 107, 163, and 142 sequences were either significantly up- or down-regulated at 0 h vs 24 h, 24 h vs 48 h, 48 h vs 1 week, and 1 week vs 4 weeks, respectively. A total of 10 DEGs coding for transcription factors were significantly down-regulated, including five ethylene-responsive transcription factors, and two probable WRKY transcription factors. RT-qPCR was used to validate RNA-seq data: Spearman and Pearson correlation coefficients between SeqMonk LFC and RT-qPCR LFC were 0.748 and 0.879, respectively. To the best of our knowledge, this is the first transcriptomic assessment of the changes taking place in the mRNA profile of a developing potato stem explant.


Introduction
A search on any major database for studies related to in vitro plant culture will reveal thousands of papers that describe the response of explants to different growth conditions. What is not so common are studies that have assessed the molecular or transcriptional changes that occur in a developing explant from the moment it is cut and plated, through initial stages of development until the end of subculture, even though considerable evidence exists about the genetic switch of cell fate during in vitro organogenesis (Zhao et al. 2008). Che et al. (2006) built a transcriptomic map of the early developmental events in Arabidopsis thaliana during callus, root, or shoot regeneration after exposure to an auxin-rich callus induction phase followed by a cytokinin-rich shoot induction phase. They found that the biological, cellular, and molecular profile changed depending on the organogenic pathway that had been selected. An earlier study by the same group found that hundreds of genes were expressed during early stages of shoot development (Che et al. 2002). Pischke et al. (2006) noted that histidine kinases were over-regulated in habituated callus, i.e., callus that develops independent of cytokinin during successive subcultures, while several genes involved in cytokinin signaling also showed altered expression. The Arabidopsis Branching Enzyme 1 (BE1) gene, which encodes a glycoside hydrolase that is involved in carbohydrate metabolism, determines the fate of organogenesis and somatic embryogenesis in A. thaliana (Wang et al. 2014). The regenerative ability of cells that underlies totipotency and dedifferentiation lies in the regulation of Electronic supplementary material The online version of this article (https ://doi.org/10.1007/s0034 4-020-10133 -0) contains supplementary material, which is available to authorized users. RNA metabolism (Ohtani 2015). Morphogenesis is under the control of a number of mechanisms that regulate gene activity, including the activity of transcription factors (TFs) and microRNAs (Zhuravlev and Omelko 2008).
The ability to identify processes that are activated or repressed during different stages of growth and development of an in vitro plant is of basic biological interest, in terms of development, but could also provide valuable clues for specific metabolic engineering that would target specific processes for up-or down-regulation, either to manipulate growth or physiology, or to overexpress desired metabolites. In this study, we wished to understand the metabolic changes taking place in in vitro stem explants of potato (Solanum tuberosum L.) during a four-week interval, which is the time it takes for axillary shoots and roots to develop. Recent advances using this model plant revealed that when a stem explant is cut, 189 genes are up-regulated, including the increased production of products related to defense, stress response, and wound healing (Teixeira da . Although several stress-related genes were up-regulated in response to ultrasound (US), even during a 4-week developmental period following explant plating, resulting plantlets survived, both after US transmitted by air ) and liquid-transmitted US (PE-US; Teixeira da Silva et al. 2020). A better understanding of the transcriptomic profile would also allow for the developmental fate to be better understood  if the same explants were to be placed on different media with different plant growth regulators, thereby inducing different organs.

Plant material and Growth Conditions
Plant material as well as culture and growth conditions follow those employed by Dobránszki et al. (2017Dobránszki et al. ( , 2019 and Teixeira da Silva et al. ( , 2020. Briefly, singlenode segments containing a single leaf were excised from 4-week-old in vitro potato plantlets of 'Desirée' and plated for 4 weeks on plant growth regulator-free Murashige and Skoog (1962) medium. Plant material was sampled at 0 h (explants); 24 h, 48 h (explants with a developing axillary shoot and adventitious roots); 1 week and 4 weeks (developing explants and growing plantlets). Sampling times based on the growing phases were as determined in Dobránszki et al. (2017) although an additional sampling time was inserted at 48 h after explants were subcultured. By 48 h after the start of subculture, rooting percentage of explants exceeded 50%. Samples (explants at 0, 24 and 48 h; leaves, stems, and roots (i.e., plantlets) at 1 week and 4 weeks; Fig. 1) were stored immediately at -80 °C until further analysis.

mRNA Isolation, Library Construction, and Sequencing
Total RNA was purified from the samples as three biological replicates. Total RNA isolation and removal of rRNA were made as described earlier Teixeira da Silva et al. 2019, 2020. The 15 mRNA libraries derived from the three biological replicates of five treatments were pooled into five qualified libraries. Deep sequencing (150 bp paired-end reads, expected data volume was 100 M reads/ sample) was conducted as in our earlier studies Teixeira da Silva et al. 2019, 2020.

Validation of DEGs by RT-qPCR
Total RNA was isolated from five samples as three biological replicates as described earlier Teixeira da Silva et al. 2019, 2020. For the RT-qPCR analysis, we selected eight DEGs from the RNA-seq datasets (PGSC0003DMG400003058, PGSC0003DMG40006825, ENSRNA49463513, PGSC0003DMG400010762, PGSC0003DMG400027054, PGSC0003DMG400028022, PGSC0003DMG400026885, PGSC0003DMG400046976) based on the most negative and positive changes in intensity in the SeqMonk logarithmic fold change (LFC) values. Five reference genes (EF1α, elongation factor-1alpha; actin; tubulin; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; sec3, exocyst complex component sec3) were selected for RT-qPCR validation based on a study by   (Teixeira da Silva et al. 2020). RT-qPCR primers (Suppl. Table 4) were developed for the chosen DEGs and the five normalizing (reference) genes (GAPDH, Actin, sec3, EF1α, tubulin) with the CLC Main Workbench 7.9.2 (Qiagen, Hilden, Germany). geNorm (Vandesompele et al. 2002), NormFinder (Andersen et al. 2004), and Best-Keeper (Pfaffl et al. 2004) statistical methods were used to compare the stability of expression intensity among candidate normalizing genes based on the cycle quantification (Cq) value. The results were compared in RefFinder (Pfaffl et al. 2004) based on the geometric mean of the rankings of every single gene calculated by geNorm, NormFinder, and BestKeeper. RT-qPCR was performed and evaluated as in our previous studies Teixeira da Silva et al. 2019, 2020. Correlation coefficients (Spearman and Pearson) were calculated using Excel in Microsoft Office 2018 (Microsoft, Redmond, WA, USA).

Global Changes in RNA Expression Profiles
The expression intensity of 40,430 genes was assessed. A total of 158, 107, 163, and 142 sequences were either significantly up-or down-regulated at 0 h vs 24 h, 24 h vs 48 h, 48 h vs 1 week, and 1 week vs 4 weeks, respectively (Suppl.   up-regulated genes (Table 1). Scatter plots show the correlation between any two sets of treatment pairs (Suppl. Fig. 1).
Heat maps based on a per-probe normalization (Suppl. Fig. 2) were generated.

RNA-seq/DEG Analysis of Biological, Cellular, and Molecular Processes
When comparing the significantly up-and down-regulated biological, cellular, and molecular processes (Suppl. Fig. 3; Suppl. Table 2), a cut-off value of ≥ 50% of DEGs was applied in the comparison across all treatments so as to focus on the most important (weighted) processes. Among the 23 biological processes, the vast majority of significantly upand down-regulated DEGs were limited to six categories (carbohydrate metabolism, catabolism, cellular component organization, cellular protein modification, response to biotic stimulus, and response to stress), the latter process displaying significantly up-and down-regulated DEGs in 100% of the comparisons. By 24 h after subculture, six processes were down-regulated while seven processes were upregulated. In the transition between different growth periods, 12 biological processes were down-regulated and seven processes were up-regulated in the transition between 24 and 48 h, while 10 processes each for the transition between 48 h and 1 week, as well as between 1 and 4 weeks, were up-regulated and 7 processes were down-regulated. Only one significantly up-regulated DEG accounted for a high relative percentage value (≥ 30%), namely the response to stress at the transition between 0 and 24 h (Suppl. Table 2). Among the 13 cellular processes, the vast majority of significantly up-and down-regulated DEGs were limited to three locations (cell wall, extracellular region, and plastids), the plastids displaying significantly up-and down-regulated DEGs in 50% of the comparisons. As plants developed in vitro, the number of cellular processes that were up-or down-regulated accounted for less than 50% of all cellular processes, with the highest number being 5 up-regulated processes between 0 and 24 h and between 24 and 48 h, 6 down-regulated processes between 24 and 48 h, and 5 downregulated processes between 1 and 4 weeks. When considering a cut-off value of ≥ 35%, six developmental transitions in three locations in in vitro development showed a significantly high proportion of down-regulated DEGs, namely in the extracellular region, membrane, and nucleus, as well as a significantly high proportion of up-regulated DEGs, namely in four developmental transitions in two locations (extracellular region and membrane) (Suppl. Table 2).
Among the 15 molecular processes, the vast majority of significantly up-and down-regulated DEGs were limited to four activities (enzyme regulator activity, hydrolase activity, nucleotide binding, and protein binding), nucleotide binding accounting for 100% of significantly up-and down-regulated DEGs. Only in the transition between 0 and 24 h and between 24 and 48 of plant development in vitro did the number of molecular processes that were up-or down-regulated accounted for more than 50% of all molecular processes, specifically 8 up-regulated processes and 9 down-regulated processes, respectively. When considering the same relative level as biological processes (i.e., ≥ 30%), hydrolase activity was significantly up-regulated during in vitro development in the transition between 0 and 24 h and between 24 and 48 h, but significantly down-regulated in the transition between 48 h and 1 week. Protein binding was significantly up-regulated during the transition between 48 h and 1 week, but significantly down-regulated in the transition between 1 and 4 weeks (Suppl. Table 2).

Up-and Down-Regulated DEGs Related to Metabolic Processes and Transcription Factors
When comparing five sampling times, i.e., four comparisons (0 h vs 24 h, 24 h vs 48 h, 48 h vs 1 week, 1 week vs 4 weeks), the focus was only on significantly up-or downregulated DEGs related to growth and development (Table 2; Suppl. Table 3). In amino acid metabolism, DEGs related to nine assumed metabolic processes were either up-or downregulated between 0 h and 4 weeks (two processes for biosynthesis of other secondary metabolites, eight processes for carbohydrate metabolism, three processes for energy metabolism, five processes for lipid metabolism, four processes for metabolism of terpenoids and polyketides, and two processes for vitamins). DEGs related to aminobenzoate degradation were also up-and down-regulated over the same time interval.
A total of 10 DEGs coding for transcription factors (TFs) were significantly down-regulated, including five ethyleneresponsive TFs, and two probable WRKY TFs (Suppl. Table 5). The vast majority (80%) were down-regulated between 0 and 24 h, with the remainder (20%) down-regulated between 48 h and 1 week.

Validity of DEG Analysis by RT-qPCR
RT-qPCR was used to validate RNA-seq data. The Ref-Finder results of the reference genes are shown in Suppl. Table 6. Under in vitro growing conditions, the ranking order (from the most stable to the least stable) was GAPDH > actin > tubulin > sec3 > EF1α. We selected the GAPDH gene as the reference gene for DEG validation, as described previously (Teixeira da Silva et al. 2020). According to the RT-qPCR results, all of the chosen DEGs were detected as true-positive up-or down-regulated DEGs (Suppl. Fig. 4; Suppl. Table 6). The correlation coefficients between SeqMonk LFC and RT-qPCR LFC were 0.748 and 0.879 according to Spearman and Pearson correlations,  Fig. 4). These high correlation coefficients indicate a high positive correlation between SeqMonk LFC and RT-qPCR LFC.

Discussion
This study aimed to assess the changes taking place in freshly prepared stem internode explants that are then grown in medium free of plant growth regulators over 4 weeks, when roots and axillary shoots form. To achieve this, we used KEGG-based classification of metabolic processes based on RNA-seq analysis and DEGs, i.e., a transcriptomic analysis. Transcriptomic profiling of a cut stem explant (Teixeira da Silva et al. 2019) and of an explant that had been exposed to air-based US  or liquid-based PE-US (Teixeira da Silva et al. 2020) and then left to grow for 4 weeks revealed the up-regulation of several stress-related genes. Plantlets exposed to PE-US showed significantly longer roots, shoot and root fresh weight, and chlorophyll a/b ratio than control plantlets (Teixeira da Silva et al. 2020).

Possible Metabolic Changes Related to Upand Down-Regulated DEGs Between 0 and 24 h of the Subculture
The transcriptomic changes during the first 24 h of subculture, after placing the explant onto medium, were connected to amino acid, carbohydrate, energy, lipid, terpenoid, and polyketide metabolism and biosynthesis of flavonoids and phenylpropanoids.
Polyphenol oxidases (PPOs) can assist in the oxidation of l-3,4-dihydroxyphenylalanine (L-DOPA)-a non-protein amino acid and inhibitor of electron transport in plants (Mustaq et al. 2013)-to form melanin, producing reactive oxygen species in the process (Soares et al. 2014). While most PPOs are associated with the wounding response, usually in response to abiotic stress, some PPOs are linked with the biosynthesis of metabolites, such as L-DOPA in betalain biosynthesis, tyrosine metabolism in walnut, and formation of 8-8ʹ-linked lignans in creosote bush (Sullivan 2015). The phenylpropanoid pathway has metabolic cross-talk with the biosynthesis of lignins and flavonoids in plants (Li et al. 2010;Mouradov and Spangenberg 2014). The shikimate O-hydroxycinnamoyltransferase enzyme, p-coumaroyl-CoA, was down-regulated between 0 and 24 h. In the same period, lactoperoxidase (EC. 1.11.1.7) was similarly up-regulated as was observed immediately after cutting plantlets into explants (Teixeira da Silva et al. 2019). This might increase the production of lignins provided that the precursors are available. The substrate of p-Coumaroyl-CoA has a role in lignin and flavonoid biosynthesis (Vogt, 2010).
When p-coumaroyl-CoA is down-regulated, this results in a decrease in caffeoyl CoA, which is a precursor of S (syringyl) and G (guaiacyl) units of lignin (Hoffmann et al. 2004). In Arabidopsis, reduction in the activity of this enzyme results in reduced lignin production, but the accumulation of flavonoids, and therefore affects the transport of auxin (Besseau et al. 2007;Li et al. 2010). Up-regulated β-glucosidase at the beginning of subculture (0 to 24 h) might lead to an increase in coumarin, which is an anti-auxin component in plants (Chattha et al. 2016). In later phases of subculture, β-glucosidase was down-regulated, which might result in the decreased production of coumarin.
Chitinase was the only enzyme that was down-regulated immediately after explant preparation (i.e., when whole plantlets were cut into explants) causing a decrease in NAG production immediately after cutting (Teixeira da Silva et al. 2019). Thereafter, the up-regulated chitinase, not only between 0 and 24 h, but in all comparisons, i.e., throughout the whole 4-week subculture period, might result in the increased production of N-acetylglucosamine (NAG), which plays a signaling role during plant development (Konopka 2012).
When explants were prepared from whole plantlets, the level of galactinol increased after the up-regulation of DEG for inositol 3-α-galactosyltransferase, as described earlier (Teixeira da Silva et al. 2019). However, it was down-regulated between 0 and 24 h, causing decreased galactinol production by 24 h.
Down-regulation of Rubisco might decrease the amount of glycolate, an intermediate of the glycolate pathway (photorespiration) (Allan et al. 2009). Down-regulation of the Rubisco sequence implies not only the modification of glyoxalate and dicarboxylate metabolism, but also the downregulation of carbon fixation (Allan et al. 2009). Catalase, which catalyzes the glycolate glyoxylate conversion, was upregulated at the same time. The β-glucosidase in cyanoamino acid metabolism is involved in the conversion of β-d-glucose residues to β-d-glucose. A closely related homolog sequence for the potato β-glucosidase and involved in cyanoamino acid metabolism was down-regulated in tomato in response to salt stress ). The up-regulation of DEGs for β-glucosidase, glucan endo-1-3-β-d-glucosidase causes an increase in sugar (glucose and sucrose) production.
Plant oxylipins are oxygenated fatty acids that have a wide range of biological functions in plant growth, development, and protection against stresses (Griffiths 2015). Cytochrome P450, a monooxygenase, catalyzes the metabolism of oxylipins (Griffiths 2015). Lipoxygenases catalyze the oxidation of lipids, which serve as a storage and energy source (Andreou and Feussner 2009). Phospholipase A1 catalyzes the first reaction of the hydrolysis of arachidonic acid (Hanna and Hafez 2018). Increased production of 2-oxobutanoate would result in increased levels of isoleucine and jasmonoyl-L-isoleucine (JA-Ile), which is a key hormonal signal (Widemann et al. 2013). Down-regulation of linoleate 13-S-lipoxygenase from 0 to 24 h caused the decreased production of death acid (10-oxo-11,15-phytodienoic acid, 10-OPDA), which increased in response to the stress caused by explant cutting and preparation at 0 h (Teixeira da . Similarly, as their precursor, its down-regulation decreased the production of methyl jasmonate and traumatic acid, both of which participate in wound healing. Threonine ammonia lyase, or threonine deaminase, which is responsible for maintaining the equilibrium of isoleucine in plants together with methionine γ-lyase (Joshi et al. 2010), was up-regulated in glycine, serine and threonine metabolism, as well as in valine, leucine, and isoleucine biosynthesis, between 0 and 24 h. Thioredoxin reductase regulates photosynthetic enzymes in chloroplasts (Arnér and Holmgren 2000). The reduction of thioredoxin by NADPHdependent thioredoxin reductases may play a part in embryo, root, and shoot apical meristem development (Reichheld et al. 2007;Yano 2014). Catalase, involved in tryptophan metabolism, and which breaks down hydrogen peroxide, for example, in response to abiotic stress (Caverzan et al. 2012), was up-regulated between 0 and 24 h, a similar response that was observed for glyoxylate and dicarboxylate metabolism (carbohydrate metabolism).
The cleavage of 9-cis-epoxycarotenoids to apocarotenoid and xanthoxin, catalyzed by 9-cis-epoxycarotenoid dioxygenase (NCED), is the precursor to the abscisic acid (ABA) biosynthetic pathway (Han et al. 2004). NCED was up-regulated in response to stress caused by explant cutting at 0 h (Teixeira da Silva et al. 2019). However, 24 h after placing explants onto regeneration medium, down-regulation of the DEG responsible for NCED caused decreased production of xanthoxin and thus ABA. The A. thaliana NCED3 mutant is unable to accumulate ABA in water-stressed plants, possibly as a result of increased water loss due to impaired stomatal function (Ruggiero et al. 2004), and enhances resistance to multiple abiotic stresses in rice (Huang et al. 2018). Downregulation of NCED results in decreased ABA production. 2-Oxoglutarate-dependent dioxygenase, such as gibberellin 2β-dioxygenase, is involved in the catabolism and inactivation of gibberellin (Rieu et al. 2008). When working in conjunction with a transcription factor, BRANCHED1, NCED3 suppresses bud development in A. thaliana as a result of the accumulation of ABA (González-Grandío et al. 2016). Down-regulation of gibberellin 2β-dioxygenase results in increased level of gibberellins.
The above-mentioned changes, which are summarized in Fig. 1, including the accumulation of sugar, modified auxin transport, and an increase in the level of coumarin (an antiauxin agent), an increase in the level of 2-oxobutanoate leading to the increased production of isoleucine and jasmonoyll-isoleucine, an increase in the level of gibberellin, increased production of NAG, fatty acid, L-DOPA but a decrease of thioredoxin and a decrease in the production of methyl jasmonate, traumatic acid, and xanthoxin (and thereby ABA) collectively indicate that explants mitigated the stress caused by explant cutting by 24 h, and initiated the development and growth of new axillary shoots and roots or root meristems from single-node explants. Shoot development was observed in 15% of explants (vs. root development in 6% of explants).

Possible Metabolic Changes Related to Upand Down-Regulated DEGs Between 24 and 48 h of the Subculture
Compared to the transcriptomic changes of the first 24 h, between 24 and 48 h, also the metabolisms of vitamins, nucleotide, and xenobiotics were affected but the metabolism of terpenoids and polyketides was not affected.
Asparagine synthetase (AS) catalyzes the formation of asparagine and glutamate (Gaufichon et al. 2016), the former being important for the release of nitrogen for amino acid and protein synthesis (Sieciechowicz et al. 1998). Its downregulation decreases the levels of asparagine and glutamate. Acetate-CoA ligase or acetate-CoA synthetase (ACS) sets the genetic code in the first step of protein synthesis by linking amino acids and tRNA (Pang et al. 2014).
The level of galactinol may have increased again by 48 h if one considers the up-regulation of the DEG responsible for inositol 3-α-galactosyltransferase in the 24 h vs. 48 h comparison. The raffinose family group of oligosaccharides in plants participates in sugar transport and sugar storage, or they act as signaling molecules, its synthesis starting with the formation of galactinol catalyzed by inositol 3-α-galactosyltransferase (Sengupta et al. 2015).
Up-regulation of PI3P synthase and 1-phosphatidylinositol 4-kinase might lead to increased production of myoinositol and different forms of phosphatidyl-inositols. Inositol phospholipids are important regulators of plant growth (Stevenson et al. 2000), and play a role as signals in intracellular communication and in response to environmental changes (Adepoju et al. 2017). Up-regulation of the 1-phosphatidylinositol 4-kinase sequence might induce changes in energy metabolism and signaling in the plants by overproduction of phosphatidyl-inositol-4-phosphate (PI4P). PI4P has a regulatory role as part of the phosphatidyl-inositol signaling system, otherwise it is the precursor of different phosphatidyl-inositol-biphosphates (PI(4,5)P 2 and PI(3,4)P 2 ) (Stevenson et al. 2000;Williams et al. 2015;Heilmann 2016). Different forms of phosphatidyl-inositol phosphates (mainly 3-and 4-phosphates (PI3P, PI4P)) play specific roles in membrane biogenesis, as well as vesicle and membrane trafficking from the endoplasmic reticulum to the Golgi-apparatus and plasmalemma (Stevenson et al. 2000). PI(4,5)P 2 is a precursor of second messengers (Stevenson et al. 2000). Myo-inositol plays an important role in the biosynthesis of the cell wall, storage of phosphate, inositol and phytic acid, intercellular communication, as part of the phosphatidyl-inositol signaling pathway, and in the storage and transport of auxins (Loewus and Murthy 2000;Williams et al. 2015). Mannanase is linked to the onset of intensive growth from explants by cutting structural mannans in the cell wall, thereby loosening the cell wall during growth (Schröder et al. 2009). Up-regulated pectin demethoxylase and pectin depolymerase might result in the overproduction of d-galacturonate whose pathway can serve as an alternative pathway for the production of ascorbic acid, which is necessary for cell expansion during plant growth, in the Solanaceae (Rigano et al. 2018). Pectate lyase-like genes are necessary for pectin degradation during growth in plants (Leng et al. 2017).
Nucleoside-triphosphate phosphatase was up-regulated at 0 h in response to explant preparation (cutting) (Teixeira da Silva et al. 2019) but its down-regulation started at 48 h of tissue culture and increased during subculture causing a change in thiamine production. Down-regulated nucleosidetriphosphate phosphatase in vitamin B1 metabolism (thiamine metabolism) might cause the increased conversion of thiamine diphosphate to thiamine phosphate (Makarchikov 2009). However, sequences coding for adenosine diphosphatase were not affected. Down-regulated acid phosphatase might inhibit the conversion of flavin mononucleotide (FMN) to riboflavin and promote the highest level of FMN and FAD as both are co-factors of several oxidoreductases (Gerdes et al. 2012). However, sequences responsible for FMN hydrolase, which catalyzes the same conversion, were not affected.
Transcriptomic changes detected between 24 and 48 h ( Fig. 1) are mainly related to increased signaling and regulatory activities (increased production of galactinol, myoinositol, phosphatidyl-inositols) pertaining to plant growth and development, increased inter-and intracellular communication (increased production of phosphatidyl-inositols, myo-inositol), and cell wall biosynthesis and secondary cell wall formation (increased activity of mannases, increased production of ascorbic acid and pectin). Increased glucose and sucrose production and starch degradation serve as sources of energy for developmental and growth processes, as was demonstrated for the previous time period (between 0 and 24 h). Down-regulation of asparagine-synthase affects both auxin production and asparagine and glutamate levels. Aminoacyl-tRNA biosynthesis and thus protein synthesis increased, partially expressed as increased growth and development. All explants formed axillary shoots, and almost 70% of these also developed adventitious roots (Fig. 1).

Possible Metabolic Changes Related to Upand Downregulated DEGs Between 48 h and 1st
Week of the Subculture Between 48 h and 1 week of subculture, during the intensive growth of axillary shoots and adventitious roots when all explants had already formed shoots that had rooted (Fig. 1), only the metabolism of terpenoids and polyketides was not affected out of processes studied and listed in the Table 2. The metabolism of cyanoamino acid and tryptophan and the biosynthesis of aminoacyl tRNA, which were up-regulated earlier, were down-regulated during this time interval, similarly to the phenylpropanoid biosynthesis, arachidonic acid and linoleic acid metabolisms, fatty acid degradation.
Down-regulated NADH dehydrogenase and NADH:ubiquinone reductase (H + -translocating) might result in a decrease or inhibition of H + translocation and thereby the electrochemical proton gradient through complex I of the electron transport chain in the inner membrane of mitochondria (Hüttemann et al. 2007;Subrahmanian et al. 2016). However, alternative pathways in plants, such as alternative oxidases, or alternative NAD(P)H dehydrogenases, can transport electrons to the electron transport chain by bypassing the route through complex I. The by-passing pathways for electron transfer are dynamically regulated in plants and are based on the plant's metabolic state (Schertl and Braun 2014;Toro and Pinto 2015;Subrahmanian et al. 2016).
Down-regulation of pectin demethoxylase and pectate lyase decreased the degradation of pectin. DEGs responsible for β-glucosidase, glucan endo-1-3-β-d-glucosidase, sucrose synthase, and glycogen phosphorylase, which all were upregulated during the first 48 h of the subculture, were downregulated between 48 h and 1 week. Down-regulation of glucose-1-phosphate adenylyltransferase results in decrease in the production of ADP-glucose and cross-talks to starch and sugar metabolism. Sugar can modify the growth and development of plants including the expansion and division of cells at many levels. For example, it affects the transcription of many genes, maintenance of meristem function, participates in signaling, provides energy, affects biomass, and modulates the level, transport, and/or signal transduction of auxins (Wang and Ruan 2013;Lastdrager et al. 2014;Ivakov et al. 2017), and participates and cooperates in signaling of other hormones like abscisic acid and ethylene (Eveland and Jackson 2012).

Possible Metabolic Changes Related to Upand Downregulated DEGs Between 1st and 4th
Week of the Subculture Cysteine and cysteine-related molecules are involved in photosystem protection, as well as in root hair development (Romero et al. 2014). DNA (cytosine-5-)-methyltransferase, involved in cysteine and methionine metabolism mediating the [S-adenosyl-l-methionine + DNA containing cytosine = S-adenosyl-l-homocysteine + DNA containing 5-methylcytosine] reaction, was down-regulated between 1 and 4 weeks. Up-regulation of linoleate 13-S-lipoxygenase might cause change in the lipid oxidation by the end of the subculture. Vetispiradiene synthase is involved in cyclization during phytoalexin biosynthesis (Bohlmann et al. 1998;Ginzberg et al. 2009) was up-regulated. Cytokinin dehydrogenase, which was up-regulated, catalyzes the irreversible degradation, via oxidation, and thus inactivation of cytokinins (Frébort et al. 2011;Schäfer et al. 2015). As a result of the down-regulation of β-glucosidase, the production of sugars (glucose, sucrose) decreased, storing energy in macromolecules.

Transcription Factors
In our study, DEGs for 10 TFs were significantly and mostly down-regulated. Eight of them (ethylene-responsive TF ABR1-like isoform X1, ERF020-like, ERF027-like, ERF061, ERF109-like, probable WRKY TF 40 and 53, and TF bHLH35) were significantly up-regulated immediately after explant preparation (0 h) in response to cutting (Teixeira da Silva et al. 2019) but down-regulated after explants were cultured on medium for 24 h. Potato explants showed the expression of 21 TF-related DEGs in response to PE-US (Teixeira da Silva et al. 2020). Che et al. (2006) found that an ethylene-responsive factor gene was up-regulated during early adventitious shoot development in in vitro Arabidopsis thaliana.

Conclusion
A transcriptomic profile now exists for a cut potato stem explant (Teixeira da Silva et al. 2019) and for an explant that had been exposed to US (ultrasound transmitted by air)  or PE-US (ultrasound transmitted by liquid) (Teixeira da Silva et al. 2020) and then developed for 4 weeks. This study aimed to better understand the transcriptomic changes in a developing single-node stem potato explant during in vitro growth until 4 weeks of age. To achieve this, DEGs pertaining to metabolic processes and their respective theoretical enzymes and/or products were assessed (Fig. 1). We found that 80% of these processes were clearly linked to plant growth and development as described in detail in the discussion section. Compared to our earlier finding (Teixeira da Silva et al. 2019), we observed changes in DEGs of stress-related enzymes, since the stress caused by explant preparation ceased after explants were cultured for 24 h on medium. Between 0 and 24 h, when 15% and 6% of explants had formed shoots and roots, respectively, production of 2-oxobutanoate, L-DOPA, β-glucosidase, coumarin, NAG, fatty acid from phospholipids, gibberellin (GA), and the activity of PPO increased, whereas production of lignin, galactinol, Rubisco, lipid oxidation, death acid, and ABA decreased. Between 24 and 48 h, when 100% and 69% of explants had formed shoots and roots, respectively, production of NAG, D-galacturonate (ascorbic acid), galactinol, myo-inositol, and different forms of phosphatidyl inositols (PIs), BG, and oxylipins increased, as did growth, and formation of secondary cell walls, whereas production of asparagine, glutamate, IAA, coumarin, and lipid oxidation decreased. Between 48 h and 1 week (at 1 week, 100% of explants had also already formed roots), production of NAG, FMN, and FAD increased, whereas production of ACS, coumarin and BG decreased as did the degradation of pectin, auxin biosynthesis, H + translocation, and lipid oxidation. Between 1 and 4 weeks, when several nodes developed with fully developed leaves and an adventitious root system, production of NAG, ACS, JA-precursors, phytoalexin (vetispiradiene) increased, as did aminoacyl (tryptophanyl)-tRNA synthesis, degradation of cytokinins, and lipid oxidation, whereas production of coumarin and BG decreased, as did PPO activity. GSE123176, BioProject ID PRJNA507769, and SRA ID SRP171630 for the five samples: GSM3498049, GSM3498050, GSM3498051, GSM3498052, GSM3498053.

Compliance with Ethical Standards
Conflicts of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.