Histological and Transcriptomic Analysis during Bulbil Formation in Lilium lancifolium

Endogenous Hormone Levels and Transcriptomic Analysis Reveal the Mechanisms of Bulbil Initiation in Pinellia ternata

Pinellia ternata is a medicinal plant that has important pharmacological value, and the bulbils serve as the primary reproductive organ; however, the mechanisms underlying bulbil initiation remain unclear. Here, we characterized bulbil development via histological, transcriptomic, and targeted metabolomic analyses to unearth the intricate relationship between hormones, genes, and bulbil development. The results show that the bulbils initiate growth from the leaf axillary meristem (AM). In this stage, jasmonic acid (JA), abscisic acid (ABA), isopentenyl adenosine (IPA), and salicylic acid (SA) were highly enriched, while indole-3-acetic acid (IAA), zeatin, methyl jasmonate (MeJA), and 5-dexoxystrigol (5-DS) were notably decreased. Through OPLS-DA analysis, SA has emerged as the most crucial factor in initiating and positively regulating bulbil formation. Furthermore, a strong association between IPA and SA was observed during bulbil initiation. The transcriptional changes in IPT (Isopentenyltransferase), CRE1 (Cytokinin Response 1), A-ARR (Type-A Arabidopsis Response Regulator), B-ARR (Type-B Arabidopsis Response Regulator), AUX1 (Auxin Resistant 1), ARF (Auxin Response Factor), AUX/IAA (Auxin/Indole-3-acetic acid), GH3 (Gretchen Hagen 3), SAUR (Small Auxin Up RNA), GA2ox (Gibberellin 2-oxidase), GA20ox (Gibberellin 20-oxidase), AOS (Allene oxide synthase), AOC (Allene oxide cyclase), OPR (Oxophytodienoate Reductase), JMT (JA carboxy l Methyltransferase), COI1 (Coronatine Insensitive 1), JAZ (Jasmonate ZIM-domain), MYC2 (Myelocytomatosis 2), D27 (DWARF27), SMAX (Suppressor of MAX2), PAL (Phenylalanine Ammonia-Lyase), ICS (Isochorismate Synthase), NPR1 (Non-expressor of Pathogenesis-related Genes1), TGA (TGACG Sequence-specific Binding), PR-1 (Pathogenesis-related), MCSU (Molybdenium Cofactor Sulfurase), PP2C (Protein Phosphatase 2C), and SnRK (Sucrose Non-fermenting-related Protein Kinase 2) were highly correlated with hormone concentrations, indicating that bulbil initiation is coordinately controlled by multiple phytohormones. Notably, eight TFs (transcription factors) that regulate AM initiation have been identified as pivotal regulators of bulbil formation. Among these, WUS (WUSCHEL), CLV (CLAVATA), ATH1 (Arabidopsis Thaliana Homeobox Gene 1), and RAX (Regulator of Axillary meristems) have been observed to exhibit elevated expression levels. Conversely, LEAFY demonstrated contrasting expression patterns. The intricate expression profiles of these TFs are closely associated with the upregulated expression of KNOX(KNOTTED-like homeobox), suggesting a intricate regulatory network underlying the complex process of bulbil initiation. This study offers a profound understanding of the bulbil initiation process and could potentially aid in refining molecular breeding techniques specific to P. ternata.

Starch and sucrose metabolism plays an important role in the stem development in Medicago sativa

The forage quality of alfalfa (Medicago sativa ) stems is greater than the leaves. Sucrose hydrolysis provides energy for stem development, with starch being enzymatically converted into sucrose to maintain energy homeostasis. To understand the physiological and molecular networks controlling stem development, morphological characteristics and transcriptome profiles in the stems of two alfalfa cultivars (Zhungeer and WL168) were investigated. Based on transcriptome data, we analysed starch and sugar contents, and enzyme activity related to starch-sugar interconversion. Zhungeer stems were shorter and sturdier than WL168, resulting in significantly higher mechanical strength. Transcriptome analysis showed that starch and sucrose metabolism were significant enriched in the differentially expressed genes of stems development in both cultivars. Genes encoding INV , bglX , HK , TPS and glgC downregulated with the development of stems, while the gene encoding was AMY upregulated. Weighted gene co-expression network analysis revealed that the gene encoding glgC was pivotal in determining the variations in starch and sucrose contents between the two cultivars. Soluble carbohydrate, sucrose, and starch content of WL168 were higher than Zhungeer. Enzyme activities related to sucrose synthesis and hydrolysis (INV, bglX, HK, TPS) showed a downward trend. The change trend of enzyme activity was consistent with gene expression. WL168 stems had higher carbohydrate content than Zhungeer, which accounted for more rapid growth and taller plants. WL168 formed hollow stems were formed during rapid growth, which may be related to the redistribution of carbohydrates in the pith tissue. These results indicated that starch and sucrose metabolism play important roles in the stem development in alfalfa.

Bulbil initiation: a comprehensive review on resources, development, and utilisation, with emphasis on molecular mechanisms, advanced technologies, and future prospects

Bulbil is an important asexual reproductive structure of bulbil plants. It mainly grows in leaf axils, leaf forks, tubers and the upper and near ground ends of flower stems of plants. They play a significant role in the reproduction of numerous herbaceous plant species by serving as agents of plant propagation, energy reserves, and survival mechanisms in adverse environmental conditions. Despite extensive research on bulbil-plants regarding their resources, development mechanisms, and utilisation, a comprehensive review of bulbil is lacking, hindering progress in exploiting bulbil resources. This paper provides a systematic overview of bulbil research, including bulbil-plant resources, identification of development stages and maturity of bulbils, cellular and molecular mechanisms of bulbil development, factors influencing bulbil development, gene research related to bulbil development, multi-bulbil phenomenon and its significance, medicinal value of bulbils, breeding value of bulbils, and the application of plant tissue culture technology in bulbil production. The application value of the Temporary Immersion Bioreactor System (TIBS) and Terahertz (THz) in bulbil breeding is also discussed, offering a comprehensive blueprint for further bulbil resource development. Additionally, additive, seven areas that require attention are proposed: (1) Utilization of modern network technologies, such as plant recognition apps or websites, to collect and identify bulbous plant resources efficiently and extensively; (2) Further research on cell and tissue structures that influence bulb cell development; (3) Investigation of the network regulatory relationship between genes, proteins, metabolites, and epigenetics in bulbil development; (4) Exploration of the potential utilization value of multiple sprouts, including medicinal, ecological, and horticultural applications; (5) Innovation and optimization of the plant tissue culture system for bulbils; (6) Comprehensive application research of TIBS for large-scale expansion of bulbil production; (7) To find out the common share genetics between bulbils and flowers.

Auxin regulates bulbil initiation by mediating sucrose metabolism in Lilium lancifolium

Lily bulbils, which serve as advantageous axillary organs for vegetative propagation, have not been extensively studied in terms of the mechanism of bulbil initiation. The functions of auxin and sucrose metabolism have been implicated in axillary organ development, but their relationship in regulating bulbil initiation remains unclear. In this study, exogenous indole-3-acetic acid (IAA) treatment increased the endogenous auxin levels at leaf axils and significantly decreased bulbil number, whereas treatment with the auxin polar transport inhibitor N-1-naphthylphthalamic acid (NPA), which resulted in a low auxin concentration at leaf axils, stimulated bulbil initiation and increased bulbil number. A low level of auxin caused by NPA spraying or silencing of auxin biosynthesis genes YUCCA FLAVIN MONOOXYGENASE-LIKE 6 (LlYUC6) and TRYPTOPHAN AMINOTRANSFERASERELATED 1 (LlTAR1) facilitated sucrose metabolism by activating the expression of SUCROSE SYNTHASES 1 (LlSusy1) and CELL WALL INVERTASE 2 (LlCWIN2), resulting in enhanced bulbil initiation. Silencing LlSusy1 or LlCWIN2 hindered bulbil initiation. Moreover, the transcription factor BASIC HELIX-LOOP-HELIX 35 (LlbHLH35) directly bound the promoter of LlSusy1, but not the promoter of LlCWIN2, and activated its transcription in response to the auxin content, bridging the gap between auxin and sucrose metabolism. In conclusion, our results reveal that an LlbHLH35-LlSusy1 module mediates auxin-regulated sucrose metabolism during bulbil initiation.

Comprehensive metabolite and biological profile of "Sulmona Red Garlic" ecotype's aerial bulbils

"Sulmona Red Garlic" is a well-known Italian traditional product. Bulbs, used for culinary purposes, have been largely investigated for their medicinal properties whereas aerial bulbils are usually removed as waste material. Here, for the first time, chemical composition and biological properties of the hydroalcoholic extract from aerial bulbils were investigated. Complementary information on metabolite composition were obtained using both NMR based untargeted and HPLC-DAD targeted methodologies. The NMR analysis revealed the presence of sugars, organic acids, amino acids, organosulphur compounds (methiin, alliin, allicin and cycloalliin), and other secondary metabolites. In particular, methiin and alliin were identified for the first time in the NMR spectra of aerial bulbil garlic extracts. Polyphenol content was determined by HPLC-DAD analysis: catechin, chlorogenic acid, and gallic acid turned out to be the most abundant phenolics. Hydroalcoholic extract blocked cell proliferation of colon cancer cell line HCT116 with an IC50 of 352.07 µg/mL, while it was non-toxic to myoblast cell line C2C12. In addition, it caused seedling germination reduction of two edible and herbaceous dicotyledon species, namely Cichorium intybus and C. endivia. Moreover, the same extract reduced the gene expression of TNF-α (tumor necrosis factor), HIF1-α (hypoxia-inducible factor), VEGFA (vascular endothelial growth factor), and transient receptor potential (TRP) M8 (TRPM8) indicating the ability to contrast cancer development through the angiogenic pathway. Final, in silico experiments were also carried out supporting the biological effects of organosulphur compounds, particularly alliin, which may directly interact with TRPM8. The results here reported suggest the potential use of garlic aerial bulbils often considered a waste product as a source in phytotherapeutic remedies.

Full-length transcriptome analysis revealed that 2,4-dichlorophenoxyacetic acid promoted in vitro bulblet initiation in lily by affecting carbohydrate metabolism and auxin signaling

Bulblet initiation, including adventitious bud initiation and bulblet formation, is a crucial process for lily and other bulbous flowers that are commercially propagated by vegetative means. Here, by a hybrid strategy combining Pacific Biosciences (PacBio) full-length sequencing and Illumina RNA sequencing (RNA-seq), high-quality transcripts of L. brownii (Lb) and its variety, L. brownii var. giganteum (Lbg), during in vitro bulblet initiation were obtained. A total of 53,576 and 65,050 high-quality non-redundant full-length transcripts of Lbg and Lb were generated, respectively. Morphological observation showed that Lbg possessed a stronger capacity to generate bulblets in vitro than Lb, and 1 mg L-1 2,4-dichlorophenoxyacetic acid (2,4-D) significantly increased bulblet regeneration rate in two lilies. Screening of differentially expressed transcripts (DETs) between different stages and Mfuzz analysis showed 0 DAT to 1 DAT was the crucial stage with the most complex transcriptional change, with carbohydrate metabolism pathway was significantly enriched. In addition, 6,218 and 8,965 DETs were screened between the 2,4-D-treated group and the control group in Lbg and Lb, respectively. 2,4-D application had evident effects on the expression of genes involved in auxin signaling pathway, such as TIRs, ARFs, Aux/IAAs, GH3s and SAURs. Then, we compared the expression profiles of crucial genes of carbohydrate metabolism between different stages and different treatments. SUSs, SUTs, TPSs, AGPLs, GBSSs and SSs showed significant responses during bulblet initiation. The expression of CWINs, SUTs and SWEETs were significantly upregulated by 2,4-D in two lilies. In addition, 2,4-D increased the expression of starch degradation genes (AMYs and BAMs) and inhibited starch synthesis genes (AGPLs, GBSSs and SSs). SBEs were significantly upregulated in Lbg but not in Lb. Significant co-expression was showed between genes involved in carbohydrate metabolism and auxin signaling, together with transcription factors such as bHLHs, MYBs, ERFs and C3Hs. This study indicates the coordinate regulation of bulblet initiation by carbohydrate metabolism and auxin signaling, serving as a basis for further studies on the molecular mechanism of bulblet initiation in lily and other bulbous flowers.

mRNA-Seq and miRNA-Seq Analyses Provide Insights into the Mechanism of Pinellia ternata Bulbil Initiation Induced by Phytohormones

Pinellia ternata (Thunb.) Breit (abbreviated as P. ternata) is a plant with an important medicinal value whose yield is restricted by many factors, such as low reproductive efficiency and continuous cropping obstacles. As an essential breeding material for P. ternata growth and production, the bulbils have significant advantages such as a high survival rate and short breeding cycles. However, the location effect, influencing factors, and molecular mechanism of bulbil occurrence and formation have not been fully explored. In this study, exogenously applied phytohormones were used to induce in vitro petiole of P. ternata to produce bulbil structure. Transcriptome sequencing of mRNA and miRNA were performed in the induced petiole (TCp) and the induced bulbil (TCb). Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed for the identification of key genes and pathways involved in bulbil development. A total of 58,019 differentially expressed genes (DEGs) were identified. The GO and KEGG analysis indicated that DEGs were mainly enriched in plant hormone signal transduction and the starch and sucrose metabolism pathway. The expression profiles of miR167a, miR171a, and miR156a during bulbil induction were verified by qRT-PCR, indicating that these three miRNAs and their target genes may be involved in the process of bulbil induction and play an important role. However, further molecular biological experiments are required to confirm the functions of the identified bulbil development-related miRNAs and targets.

Bottom-Up Effects of Drought-Stressed Cotton Plants on Performance and Feeding Behavior of Aphis gossypii

Drought, a major stress for crop plants, is expected to increase in frequency due to climate change. Drought can alter crop growth and levels of secondary plant metabolites, which in turn can affect herbivores, but this latter point is still controversial. This study used three different polyethylene glycol (PEG-6000) levels (0%, 1%, and 3%) to simulate drought stress and evaluated their effects on cotton plants and the impacts on the performance of the cotton aphid Aphis gossypii. Cotton plants under drought stress showed decreased water content, above-ground biomass, and nitrogen content and increased soluble protein, soluble sugar, and tannin contents. Based on analysis of the developmental time and fecundity data from individuals and at the population level, a significantly lower fecundity and population abundance of A. gossypii were detected on cotton plants with drought stress, which supports the "plant vigor hypothesis". The poor development of A. gossypii is possibly related to lower xylem sap and phloem ingestion under drought stress. In addition, the increased tannin content of cotton plants induced by drought and lower detoxification enzyme activities of A. gossypii may have affected the responses of aphids to drought-stressed plants. Overall, the results showed that drought stress altered the physiological characteristics of the cotton plants, resulting in adverse bottom-up effects on cotton aphid performances. This implies that the adoption of drip irrigation under plastic film that can help alleviate drought stress may favor the population growth of cotton aphids.

Transcriptome Analysis Reveals the Molecular Regularity Mechanism Underlying Stem Bulblet Formation in Oriental Lily 'Siberia'; Functional Characterization of the LoLOB18 Gene

The formation of underground stem bulblets in lilies is a complex biological process which is key in their micropropagation. Generally, it involves a stem-to-bulblet transition; however, the underlying mechanism remains elusive. It is important to understand the regulatory mechanism of bulblet formation for the reproductive efficiency of Lilium. In this study, we investigated the regulatory mechanism of underground stem bulblet formation under different conditions regarding the gravity point angle of the stem, i.e., vertical (control), horizontal, and slanting. The horizontal and slanting group displayed better formation of bulblets in terms of quality and quantity compared with the control group. A transcriptome analysis revealed that sucrose and starch were key energy sources for bulblet formation, auxin and cytokinin likely promoted bulblet formation, and gibberellin inhibited bulblet formation. Based on transcriptome analysis, we identified the LoLOB18 gene, a homolog to AtLOB18, which has been proven to be related to embryogenic development. We established the stem bud growth tissue culture system of Lilium and silenced the LoLOb18 gene using the VIGS system. The results showed that the bulblet induction was reduced with down-regulation of LoLOb18, indicating the involvement of LoLOb18 in stem bulblet formation in lilies. Our research lays a solid foundation for further molecular studies on stem bulblet formation of lilies.

Comparative transcriptome and metabolome analyses identified the mode of sucrose degradation as a metabolic marker for early vegetative propagation in bulbs of Lycoris

Vegetative propagation (VP) is an important practice for production in many horticultural plants. Sugar supply constitutes the basis of VP in bulb flowers, but the underlying molecular basis remains elusive. By performing a combined sequencing technologies coupled with ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry approach for metabolic analyses, we compared two Lycoris species with contrasting regeneration rates: high-regeneration Lycoris sprengeri and low-regeneration Lycoris aurea. A comprehensive multi-omics analyses identified both expected processes involving carbohydrate metabolism and transcription factor networks, as well as the metabolic characteristics for each developmental stage. A higher abundance of the differentially expressed genes including those encoding ethylene responsive factors was detected at bulblet initiation stage compared to the late stage of bulblet development. High hexose-to-sucrose ratio correlated to bulblet formation across all the species examined, indicating its role in the VP process in Lycoris bulb. Importantly, a clear difference between cell wall invertase (CWIN)-catalyzed sucrose unloading in high-regeneration species and the sucrose synthase-catalyzed pathway in low-regeneration species was observed at the bulblet initiation stage, which was supported by findings from carboxyfluorescein tracing and quantitative real-time PCR analyses. Collectively, the findings indicate a sugar-mediated model of the regulation of VP in which high CWIN expression or activity may promote bulblet initiation via enhancing apoplasmic unloading of sucrose or sugar signals, whereas the subsequent high ratio of hexose-to-sucrose likely supports cell division characterized in the next phase of bulblet formation.

WUSCHEL-related homeobox genes cooperate with cytokinin to promote bulbil formation in Lilium lancifolium

The bulbil is an important vegetative reproductive organ in triploid tiger lily (Lilium lancifolium). Based on our previously obtained transcriptome data, we screened two WUSCHEL-related homeobox (WOX) genes closely related to bulbil formation, LlWOX9 and LlWOX11. However, the biological functions and regulatory mechanisms of LlWOX9 and LlWOX11 are unclear. In this study, we cloned the full-length coding sequences of LlWOX9 and LlWOX11. Transgenic Arabidopsis (Arabidopsis thaliana) showed increased branch numbers, and the overexpression of LlWOX9 and LlWOX11 in stem segments promoted bulbil formation, while the silencing of LlWOX9 and LlWOX11 inhibited bulbil formation, indicating that LlWOX9 and LlWOX11 are positive regulators of bulbil formation. Cytokinin type-B response regulators could bind to the promoters of LlWOX9 and LlWOX11 and promote their transcription. LlWOX11 could enhance cytokinin pathway signaling by inhibiting the transcription of type-A LlRR9. Our study enriches the understanding of the regulation of plant development by the WOX gene family and lays a foundation for further research on the molecular mechanism of bulbil formation in lily.

Full-length transcriptome-referenced analysis reveals crucial roles of hormone and wounding during induction of aerial bulbils in lily

Aerial bulbils are important vegetative reproductive organs in Lilium. They are often perpetually dormant in most Lilium species, and little is known about the induction of these vegetative structures. The world-famous Oriental hybrid lily cultivar 'Sorbonne', which blooms naturally devoid of aerial bulbils, is known for its lovely appearance and sweet fragrance. We found that decapitation stimulated the outgrowth of aerial bulbils at lower stems (LSs) and then application of low and high concentrations of IAA promoted aerial bulbils emergence around the wound at upper stems (USs) of 'Sorbonne'. However, the genetic basis of aerial bulbil induction is still unclear. Herein, 'Sorbonne' transcriptome has been sequenced for the first time using the combination of third-generation long-read and next-generation short-read technology. A total of 46,557 high-quality non-redundant full-length transcripts were generated. Transcriptomic profiling was performed on seven tissues and stems with treatments of decapitation and application of low and high concentrations of IAA, respectively. Functional annotation of 1918 DEGs within stem samples of different treatments showed that hormone signaling, sugar metabolism and wound-induced genes were crucial to bulbils outgrowth. The expression pattern of auxin-, shoot branching hormone-, plant defense hormone- and wound-inducing-related genes indicated their crucial roles in bulbil induction. Then we established five hormone- and wounding-regulated co-expression modules and identified some candidate transcriptional factors, such as MYB, bZIP, and bHLH, that may function in inducing bulbils. High connectivity was observed among hormone signaling genes, wound-induced genes, and some transcriptional factors, suggesting wound- and hormone-invoked signals exhibit extensive cross-talk and regulate bulbil initiation-associated genes via multilayered regulatory cascades. We propose that the induction of aerial bulbils at LSs after decapitation can be explained as the release of apical dominance. In contrast, the induction of aerial bulbils at the cut surface of USs after IAA application occurs via a process similar to callus formation. This study provides abundant candidate genes that will deepen our understanding of the regulation of bulbil outgrowth, paving the way for further molecular breeding of lily.

Changes of starch and sucrose content and related gene expression during the growth and development of Lanzhou lily bulb

As the main forms of carbohydrates, starch and sucrose play a vital role in the balance and coordination of various carbohydrates. Lanzhou lily is the most popular edible lily in China, mainly distributed in the central region of Gansu. To clarify the relationship between carbohydrate metabolism and bulb development of Lanzhou lily, so as to provide a basis for the promotion of the growth and development in Lanzhou lily and its important economic value, we studied lily bulbs in the squaring stage, flowering stage, half withering stage and withering stage. The plant height, fresh weight of mother and daughter bulbs continued to increase during the whole growth period and fresh weight of stem and leaf began to decrease in the half withering stage. The content of starch, sucrose and total soluble sugar in the lily mother bulb accumulated mostly in the flowering, withering and half withering stages, respectively. Starch, sucrose and total soluble sugar accumulated in the daughter bulb with the highest concentration during the withering stage. In the transcription level, sucrose synthase (SuSy1) and sucrose invertase (INV2) expressed the highest in squaring stage, and the expression was significantly higher in the mother bulb than in the daughter bulb. In flowering stage, the expression levels of soluble starch synthase (SSS1), starch-branching enzyme (SBE) and adenosine diphosphate-glucose pyrophosphorylase (AGP1) genes were higher in the mother bulb than in the daughter bulb. Altogether, our results indicate that starch and sucrose are important for the bulb growth and development of Lanzhou lily.

EST-SSR Marker Development and Full-Length Transcriptome Sequence Analysis of Tiger Lily (Lilium lancifolium Thunb)

The fast advancement and deployment of sequencing technologies after the Human Genome Project have greatly increased our knowledge of the eukaryotic genome sequences. However, due to technological concerns, high-quality genomic data has been confined to a few key organisms. Moreover, our understanding of which portions of genomes make up genes and which transcript isoforms synthesize these genes is scarce. Therefore, the current study has been designed to explore the reliability of the tiger lily (Lilium lancifolium Thunb) transcriptome. The PacBio-SMRT was used for attaining the complete transcriptomic profile. We obtained a total of 815,624 CCS (Circular Consensus Sequence) reads with an average length of 1295 bp. The tiger lily transcriptome has been sequenced for the first time using third-generation long-read technology. Furthermore, unigenes (38,707), lncRNAs (6852), and TF members (768) were determined based on the transcriptome data, followed by evaluating SSRs (3319). It has also been revealed that 105 out of 128 primer pairs effectively amplified PCR products. Around 15,608 transcripts were allocated to 25 distinct KOG Clusters, and 10,706 unigenes were grouped into 52 functional categories in the annotated transcripts. Until now, no tiger lily lncRNAs have been discovered. Results of this study may serve as an extensive set of reference transcripts and help us learn more about the transcriptomes of tiger lilies and pave the path for further research.

Hormonal Regulatory Patterns of LaKNOXs and LaBEL1 Transcription Factors Reveal Their Potential Role in Stem Bulblet Formation in LA Hybrid Lily

In lily reproduction, the mechanism of formation of bulbs has been a hot topic. However, studies on stem bulblet formation are limited. Stem bulblets, formed in the leaf axils of under- and above-ground stems, provide lilies with a strong capacity for self-propagation. First, we showed that above-ground stem bulblets can be induced by spraying 100 mg/L 6-BA on the LA hybrid lily 'Aladdin', with reduced endogenous IAA and GA₄ and a higher relative content of cytokinins. Then, expression patterns of three potential genes (two KNOTTED1-like homeobox (KNOX) and one partial BEL1-like homeobox (BELL)), during stem bulblet formation from our previous study, were determined by RT-qPCR, presenting a down-up trend in KNOXs and a rising tendency in BELL. The partial BELL gene was cloned by RACE from L. 'Aladdin' and denoted LaBEL1. Physical interactions of LaKNOX1-LaBEL1 and LaKNOX1-LaKNOX2 were confirmed by yeast two-hybrid and bimolecular fluorescence complementation assays. Furthermore, hormonal regulatory patterns of single LaKNOX1, LaKNOX2, LaBEL1, and their heterodimers, were revealed in transgenic Arabidopsis, suggesting that the massive mRNA accumulations of LaKNOX1, LaKNOX2 and LaBEL1 genes during stem bulblet formation could cause the dramatic relative increase of cytokinins and the decline of GAs and IAA. Taken together, a putative model was proposed that LaKNOX1 interacts with LaKNOX2 and LaBEL1 to regulate multiple phytohormones simultaneously for an appropriate hormonal homeostasis, which suggests their potential role in stem bulblet formation in L. 'Aladdin'.

Methane-rich water induces bulblet formation of scale cuttings in Lilium davidii var. unicolor by regulating the signal transduction of phytohormones and their levels

Previous studies have shown that methane (CH₄ ) has promoting roles in the adventitious root (AR) and lateral root (LR) formation in plants. However, whether CH₄ could trigger the bulblet formation in scale cutting of Lilium davidii var. unicolor has not been elucidated. To gain insight into the effect of CH₄ on the bulblet formation, different concentrations (1, 10, 50, and 100%) of methane-rich water (MRW) and distilled water were applied to treat the scale cuttings of Lilium. We observed that treatment with 100% MRW obviously induced the bulblet formation in scale cuttings. To explore the mechanism of CH₄ -induced bulblet formation, the transcriptome of scales was analyzed. A total of 2078 differentially expressed genes (DEGs) were identified. The DEGs were classified into different metabolic pathways, especially phenylpropanoid biosynthesis, starch and sucrose metabolism, and plant signal transduction. Of these, approximately 38 candidate DEGs involved in the plant signal transduction were further studied. In addition, the expression of AP2-ERF/ERF, WRKY, GRAS, ARF, and NAC transcription factors (TFs) was changed by MRW treatment, suggesting their potential involvement in bulblet formation. As for hormones, exogenous IAA, GA, and ABA could induce the bulblet formation. Additional experiments suggested that MRW could increase the endogenous IAA, GA, and JA levels, but decrease the levels of ABA during bulblet formation, which showed that higher IAA, GA, JA levels and lower ABA content might facilitate bulblet formation. In addition, the levels of endogenous hormones were consistent with the expression level of genes involved in phytohormone signal transduction. Overall, this study has revealed that CH₄ might improve the bulblet formation of cutting scales in Lilium by regulating the expression of genes related to phytohormone signal transduction and TFs, as well as by changing the endogenous hormone levels.

Antagonism between abscisic acid and gibberellin regulates starch synthesis and corm development in Gladiolus hybridus

Understanding corm development in flower bulbs is of importance for securing the quality of cut flowers and propagation of commercial stocks. Gladiolus is one of the most popular bulb plants worldwide. Its corm development is characterized by starch accumulation. Previous research has shown that phytohormones (especially gibberellin (GA)) are involved in tuber development. However, the relationship between abscisic acid (ABA)/GA and starch during corm development remains unclear. To gain deeper insights into the biological process of corm development, we performed a detailed anatomical characterization of different stages of corm development and analyzed phytohormone levels. Our study showed that corm development is linked to hormones (ABA and GA) and carbohydrates (sucrose and starch). Exogenous hormone treatment and silencing of endogenous hormone biosynthesis genes indicated that ABA positively regulates corm development, while GA acts as an antagonist of ABA function. A sucrose synthase gene (GhSUS2) was shown to be involved in the antagonism between ABA and GA. GhSUS2 was upregulated by ABA and downregulated by GA. The increase in the transcript level of GhSUS2 coincided with the development of corm/cormels. Silencing of GhSUS2 repressed corm development and starch accumulation. In conclusion, we propose that GhSUS2, an essential enzyme in sucrose degradation, is differentially regulated by ABA and GA and controls corm development in Gladiolus.

Molecular cloning, characterization and expression analysis of three key starch synthesis-related genes from the bulb of a rare lily germplasm, Lilium brownii var. giganteum

Starch is the predominant compound in bulb scales, and previous studies have shown that bulblet development is closely associated with starch enrichment. However, how starch synthesis affects bulbification at the molecular level is unclear. In this study, we demonstrate that Lilium brownii var. giganteum, a wild lily with a giant bulb in nature, and L. brownii, the native species, have different starch levels and characteristics according to cytological and ultra-structural observations. We cloned the complete sequence of three key gene-encoding enzymes (LbgAGPS, LbgGBSS, andLbgSSIII) during starch synthesis by rapid amplification of 5' and 3' complementary DNA (cDNA) ends (RACE) technology. Bioinformatics analysis revealed that the proteins deduced by these genes contain the canonical conserved domains. Constructed phylogenetic trees confirmed the evolutionary relationships with proteins from other species, including monocotyledons and dicotyledons. The transcript levels of various tissues and time course samples obtained during bulblet development uncovered relatively high expression levels in bulblets and gradual increase expression accompanying bulblet growth. Moreover, a set of single nucleotide polymorphisms (SNPs) was discovered in the AGPS genes of four lily genotypes, and a purifying selection fashion was predicted according to the non-synonymous/synonymous (Ka/Ks) values. Taken together, our results suggested that key starch-synthesizing genes might play important roles in bulblet development and lead to distinctive phenotypes in bulblet size.

Cytokinin Type-B Response Regulators Promote Bulbil Initiation in Lilium lancifolium

The bulbil is an important vegetative reproductive organ in triploid Lilium lancifolium whose development is promoted by cytokinins. Type-B response regulators (RRs) are critical regulators that mediate primary cytokinin responses and promote cytokinin-induced gene expression. However, the function of cytokinin type-B Arabidopsis RRs (ARRs) in regulating bulbil formation is unclear. In this study, we identified five type-B LlRRs, LlRR1, LlRR2, LlRR10, LlRR11 and LlRR12, in L. lancifolium for the first time. The five LlRRs encode proteins of 715, 675, 573, 582 and 647 amino acids. All of the regulators belong to the B-I subfamily, whose members typically contain a conserved CheY-homologous receiver (REC) domain and an Myb DNA-binding (MYB) domain at the N-terminus. As transcription factors, all five type-B LlRRs localize at the nucleus and are widely expressed in plant tissues, especially during axillary meristem (AM) formation. Functional analysis showed that type-B LlRRs are involved in bulbil formation in a functionally redundant manner and can activate LlRR9 expression. In summary, our study elucidates the process by which cytokinins regulate bulbil initiation in L. lancifolium through type-B LlRRs and lays a foundation for research on the molecular mechanism of bulbil formation in the lily.

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.

Comparative transcriptome analyses provide novel insights into etiolated shoot development of walnut (Juglans regia L.)

In general, genes promoting IAA, CTK GA and ethylene biosynthesis were upregulated, while genes participating in ABA, chlorophyll and starch biosynthesis pathways performed opposite tendency during etiolation. Etiolation as a method for rejuvenation plays an important role in the vegetative propagation of woody plants. However, the molecular mechanism of etiolated shoot development remains unclear. In this study, we investigated changes at different etiolation stages of Juglans regia. The histology and transcriptome of J. regia were analysed using etiolated stems, which were treated in darkness for 30, 60, 90 days. The results showed that the ratios of pith (Pi) diameter/stem diameter (D), cortex (Co) width/D, and phloem (Ph) width/D increased, while the ratio of xylem (Xy) width/D decreased after etiolation, and the difference in these ratios between etiolated stems and the control was more significant at 60 days than 90 days. Differentially expressed genes (DEGs) were significantly enriched in pathways such as plant hormone biosynthesis and signal transduction, chlorophyll biosynthesis and degradation, and starch and sucrose metabolism. The difference in the contents of indole-3-acetic acid (IAA), abscisic acid (ABA), sugar and chlorophyll between etiolated stems and the control increased with increasing treatment duration; in contrast, the concentrations of gibberellin (GA), zeatin (ZT), and starch, as well as the difference between the etiolated stems and control were lowest at 60 days among the three stages. On the whole, the positive effect of etiolation on the rejuvenation of walnut stems changed as the treatment period increased. The present investigation lays a foundation for future studies on the effect of etiolation on rejuvenation and for promoting the efficiency of vegetative propagation.

Mechanism of exogenous cytokinins inducing bulbil formation in Lilium lancifolium in vitro

The cytokinin pathway promotes the initiation of bulbil formation, and iPA may an important type of cytokinin during bulbil formation in Lilium lancifolium. Bulbils are important vegetative reproductive organs in triploid Lilium lancifolium. We previously showed that cytokinins are involved in bulbil formation, but how cytokinins participate in bulbil formation is not clear. In this study, bulbil formation was divided into three stages on the basis of anatomical and histological observations: the bulbil initiation stage, bulbil primordium-formation stage and bulbil structure-formation stage. The results indicated that iPA was the most critical cytokinin during the bulbil initiation. qRT-PCR revealed that increased iPA content during bulbil initiation was mainly due to increased expression of cytokinin synthesis genes (IPT1/5) and cytokinin activation genes (LOG1/3/5/7) and significantly decreased expression of the cytokinin degradation gene CKX4. Exogenous 6-BA and lovastatin affected the cytokinin pathway and promoted or inhibited bulbil initiation by increasing or decreasing the content of endogenous iPA, respectively. In summary, we demonstrate that cytokinins positively regulate bulbil formation and provide preliminary insight into the regulatory mechanisms by which the cytokinin pathway promotes bulbil initiation.

Changes in carbohydrate metabolism and endogenous hormone regulation during bulblet initiation and development in Lycoris radiata

BACKGROUND

Lycoris species have great ornamental and medicinal values; however, their low regeneration efficiency seriously restricts their commercial production. Understanding the mechanism of bulblet propagation in this genus, which has remained underexplored to date, could provide a theoretical basis for improving the reproductive efficiency. Therefore, we studied the bulblet initiation and developmental processes in Lycoris radiata.

RESULTS

We found that bulblets are formed on the junctions of the innermost layers of scales and the basal plate, and initially present as an axillary bud and gradually develop into a bulblet. We also determined the changes in carbohydrate and endogenous hormone contents during bulblet initiation and development, as well as the expression patterns of genes involved in carbohydrate metabolism and hormone biosynthesis and signaling through transcriptome analysis. Soluble sugars derived from starch degradation in the outer scales are transported to and promote bulblet initiation and development through starch synthesis in the inner scales. This process is mediated by several genes involved in carbohydrate metabolism, especially genes encoding ADP glucose pyrophosphorylase, a crucial starch synthesis enzyme. As for hormones, endogenous IAA, GA, and ABA content showed an increase and decrease during bulblet initiation and development, respectively, which were consistent with the expression patterns of genes involved in IAA, GA, and ABA synthesis and signal transduction. In addition, a decrease in ZR content may be down- and up-regulated by CK biosynthesis and degradation related genes, respectively, with increasing auxin content. Furthermore, expression levels of genes related to BR, JA, and SA biosynthesis were increased, while that of ethylene biosynthesis genes was decreased, which was also consistent with the expression patterns of their signal transduction genes.

CONCLUSIONS

The present study provides insights into the effect of carbohydrate metabolism and endogenous hormone regulation on control of L. radiata bulblet initiation and development. Based on the results, we propose several suggestions to improve L. radiata propagation efficiency in production, which will provide directions for future research.

Morphological and stage-specific transcriptome analyses reveal distinct regulatory programs underlying yam (Dioscorea alata L.) bulbil growth

In yam (Dioscorea spp) species, bulbils at leaf axils are the most striking species-specific axillary structure and exhibit important ecological niches. Genetic regulation underlying bulbil growth remains largely unclear so far. Here, we characterize yam (Dioscorea alata L.) bulbil development using histological analysis, and perform full transcriptional profiling on key developmental stages together with phytohormone analyses. Using the stage-specific scoring algorithm, we have identified 3451 stage-specifically expressed genes that exhibit a tight link between major transcriptional changes and stages. Co-expressed gene clusters revealed an obvious over-representation of genes associated with cell division and expansion at the initiation stage of bulbils (T1). Transcriptional changes of hormone-related genes highly coincided with hormone levels, indicating that bulbil initiation and growth are coordinately controlled by multiple phytohormones. In particular, localized auxin is transiently required to trigger bulbil initiation, and be further depleted or exported from bulbils to promote growth by up-regulation of genes involved in auxinconjugation and efflux. The sharp increase in supply of sucrose and an enhanced trehalose-6-phophate pathway at T1 were observed, suggesting that sucrose probably functions as a key signal and promotes bulbil initiation. Analysis of the expression of transcription factors (TFs) predicated 149 TFs as stage-specifically expressed; several T1-specific TFs (from Aux/IAA, E2F, MYB, and bHLH families) have been shown to play key roles in triggering bulbil formation. Together, our work provides a crucial angle for in-depth understanding of the molecular programs underlying yam's unique bulbil development processes. Stage-specific gene sets can be queried to obtain key candidates regulating bulbil growth, serving as valuable resources for further functional research.

WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update

WEGO (Web Gene Ontology Annotation Plot), created in 2006, is a simple but useful tool for visualizing, comparing and plotting GO (Gene Ontology) annotation results. Owing largely to the rapid development of high-throughput sequencing and the increasing acceptance of GO, WEGO has benefitted from outstanding performance regarding the number of users and citations in recent years, which motivated us to update to version 2.0. WEGO uses the GO annotation results as input. Based on GO's standardized DAG (Directed Acyclic Graph) structured vocabulary system, the number of genes corresponding to each GO ID is calculated and shown in a graphical format. WEGO 2.0 updates have targeted four aspects, aiming to provide a more efficient and up-to-date approach for comparative genomic analyses. First, the number of input files, previously limited to three, is now unlimited, allowing WEGO to analyze multiple datasets. Also added in this version are the reference datasets of nine model species that can be adopted as baselines in genomic comparative analyses. Furthermore, in the analyzing processes each Chi-square test is carried out for multiple datasets instead of every two samples. At last, WEGO 2.0 provides an additional output graph along with the traditional WEGO histogram, displaying the sorted P-values of GO terms and indicating their significant differences. At the same time, WEGO 2.0 features an entirely new user interface. WEGO is available for free at http://wego.genomics.org.cn.

Transcriptome analysis reveals the different compatibility between LAAA × AA and LAAA × LL in Lilium

To unveil the mechanism of the compatibility of odd-allotetraploid lily (LAAA) as female with diploid male lily, the differences of expressed unigenes in the ovaries and leaves between LAAA × AA and LAAA × LL were investigated using transcriptome analysis. The results showed the fruits of LAAA × AA well developed, while those of LAAA × LL aborted. The number of differentially expressed genes was less in the ovaries of LAAA × AA than those of LAAA × LL, but it showed opposite trend in those of leaves. The unigenes related with auxins, cytokinins, gibberellins, antioxidants, expansins, chlorophylls, carbohydrates, transport proteins were usually up-expressed in the ovaries and leaves of LAAA × AA but not in LAAA × LL; while those of abscisic acid, ethylene, jasmonic acid, and salicylic acid were increased in the ovaries or leaves of LAAA × LL but not in LAAA × AA. The up-expressed unigenes in the ovaries and leaves of LAAA × AA played positive roles in its fruit development because the products of the genes, like phytohormones and antioxidants, had functions protecting leaves from senescence or scavenging ROS, and thus LAAA was compatible with AA, while those of LAAA × LL played negative roles and caused its fruits aborted, and hence LAAA was incompatible with LL.

Comparison of Structural and Functional Properties of Starches from the Rhizome and Bulbil of Chinese Yam

Chinese yam is an important edible starch plant and widely cultivated in China. Its rhizome and bulbil are starch storage tissues below and above ground, respectively. In this paper, starches were isolated from the rhizome and bulbil of Chinese yam, and their structural and functional properties were compared. Both starches had an oval shape with an eccentric hilum and a CA-type crystalline structure. Their short-range ordered structure and lamellar structure had no significant difference. However, the rhizome starch had a significantly bigger granule size and lower amylose content than the bulbil starch. The swelling power and water solubility were significantly lower in the rhizome starch than in the bulbil starch. The onset and peak gelatinization temperatures were significantly higher in the rhizome starch than in the bulbil starch. The rhizome starch had a significantly higher breakdown viscosity and a lower setback viscosity than the bulbil starch. The thermal stability was lower in the rhizome starch than in the bulbil starch. The rhizome starch had a significantly lower resistance to hydrolysis and in vitro digestion than the bulbil starch. The above results provide important information for the utilization of rhizome and bulbil starches of Chinese yam.