RhMYB108, an R2R3-MYB transcription factor, is involved in ethylene- and JA-induced petal senescence in rose plants

The transcription factor RhMYB17 regulates the homeotic transformation of floral organs in rose (Rosa hybrida) under cold stress

Low temperatures affect flower development in rose (Rosa hybrida), increasing petaloid stamen number and reducing normal stamen number. We identified the low-temperature-responsive R2R3-MYB transcription factor RhMYB17, which is homologous to Arabidopsis MYB17 by similarity of protein sequences. RhMYB17 was up-regulated at low temperatures, and RhMYB17 transcripts accumulated in floral buds. Transient silencing of RhMYB17 by virus-induced gene silencing decreased petaloid stamen number and increased normal stamen number. According to the ABCDE model of floral organ identity, class A genes APETALA 1 (AP1) and AP2 contribute to sepal and petal formation. Transcription factor binding analysis identified RhMYB17 binding sites in the promoters of rose APETALA 2 (RhAP2) and APETALA 2-LIKE (RhAP2L). Yeast one-hybrid assays, dual-luciferase reporter assays, and electrophoretic mobility shift assays confirmed that RhMYB17 directly binds to the promoters of RhAP2 and RhAP2L, thereby activating their expression. RNA sequencing further demonstrated that RhMYB17 plays a pivotal role in regulating the expression of class A genes, and indirectly influences the expression of the class C gene. This study reveals a novel mechanism for the homeotic transformation of floral organs in response to low temperatures.

Tulip transcription factor TgWRKY75 activates salicylic acid and abscisic acid biosynthesis to synergistically promote petal senescence

WRKY transcription factors play a central role in controlling plant organ senescence; however, it is unclear whether and how they regulate petal senescence in the widely grown ornamental plant tulip (Tulipa gesneriana). In this study, we report that TgWRKY75 promotes petal senescence by enhancing the synthesis of both abscisic acid (ABA) and salicylic acid (SA) in tulip and in transgenic Arabidopsis. The expression level of TgWRKY75 was up-regulated in senescent petals, and exogenous ABA or SA treatment induced its expression. The endogenous contents of ABA and SA significantly increased during petal senescence and in response to TgWRKY75 overexpression. Two SA synthesis-related genes, TgICS1 and TgPAL1, were identified as direct targets of TgWRKY75, which binds to their promoters. In parallel, TgWRKY75 activated the expression of the ABA biosynthesis-related gene TgNCED3 via directly binding to its promoter region. Site mutation of the W-box core motif located in the promoters of TgICS1, TgPAL1, and TgNCED3 eliminated their interactions with TgWRKY75. In summary, our study demonstrates a dual regulation of ABA and SA biosynthesis by TgWRKY75, revealing a synergistic process of tulip petal senescence through feedback regulation between TgWRKY75 and the accumulation of ABA and SA.

The MYB Transcription Factor GmMYB78 Negatively Regulates Phytophthora sojae Resistance in Soybean

Phytophthora root rot is a devastating disease of soybean caused by Phytophthora sojae. However, the resistance mechanism is not yet clear. Our previous studies have shown that GmAP2 enhances sensitivity to P. sojae in soybean, and GmMYB78 is downregulated in the transcriptome analysis of GmAP2-overexpressing transgenic hairy roots. Here, GmMYB78 was significantly induced by P. sojae in susceptible soybean, and the overexpressing of GmMYB78 enhanced sensitivity to the pathogen, while silencing GmMYB78 enhances resistance to P. sojae, indicating that GmMYB78 is a negative regulator of P. sojae. Moreover, the jasmonic acid (JA) content and JA synthesis gene GmAOS1 was highly upregulated in GmMYB78-silencing roots and highly downregulated in overexpressing ones, suggesting that GmMYB78 could respond to P. sojae through the JA signaling pathway. Furthermore, the expression of several pathogenesis-related genes was significantly lower in GmMYB78-overexpressing roots and higher in GmMYB78-silencing ones. Additionally, we screened and identified the upstream regulator GmbHLH122 and downstream target gene GmbZIP25 of GmMYB78. GmbHLH122 was highly induced by P. sojae and could inhibit GmMYB78 expression in resistant soybean, and GmMYB78 was highly expressed to activate downstream target gene GmbZIP25 transcription in susceptible soybean. In conclusion, our data reveal that GmMYB78 triggers soybean sensitivity to P. sojae by inhibiting the JA signaling pathway and the expression of pathogenesis-related genes or through the effects of the GmbHLH122-GmMYB78-GmbZIP25 cascade pathway.

Thidiazuron combined with cyclanilide modulates hormone pathways and ROS systems in cotton, increasing defoliation at low temperatures

Low temperatures decrease the thidiazuron (TDZ) defoliation efficiency in cotton, while cyclanilide (CYC) combined with TDZ can improve the defoliation efficiency at low temperatures, but the mechanism is unknown. This study analyzed the effect of exogenous TDZ and CYC application on cotton leaf abscissions at low temperatures (daily mean temperature: 15°C) using physiology and transcriptomic analysis. The results showed that compared with the TDZ treatment, TDZ combined with CYC accelerated cotton leaf abscission and increased the defoliation rate at low temperatures. The differentially expressed genes (DEGs) in cotton abscission zones (AZs) were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses to compare the enriched GO terms and KEGG pathways between the TDZ treatment and TDZ combined with CYC treatment. TDZ combined with CYC could induce more DEGs in cotton leaf AZs at low temperatures, and these DEGs were related to plant hormone and reactive oxygen species (ROS) pathways. CYC is an auxin transport inhibitor. TDZ combined with CYC not only downregulated more auxin response related genes but also upregulated more ethylene and jasmonic acid (JA) response related genes at low temperatures, and it decreased the indole-3-acetic acid (IAA) content and increased the JA and 1-aminocyclopropane-1-carboxylic acid (ACC) contents, which enhanced cotton defoliation. In addition, compared with the TDZ treatment alone, TDZ combined with CYC upregulated the expression of respiratory burst oxidase homologs (RBOH) genes and the hydrogen peroxide content in cotton AZs at low temperatures, which accelerated cotton defoliation. These results indicated that CYC enhanced the TDZ defoliation efficiency in cotton by adjusting hormone synthesis and response related pathways (including auxin, ethylene, and JA) and ROS production at low temperatures.

Abscisic acid-induced transcription factor PsMYB306 negatively regulates tree peony bud dormancy release

Bud dormancy is a crucial strategy for perennial plants to withstand adverse winter conditions. However, the regulatory mechanism of bud dormancy in tree peony (Paeonia suffruticosa) remains largely unknown. Here, we observed dramatically reduced and increased accumulation of abscisic acid (ABA) and bioactive gibberellins (GAs) GA1 and GA3, respectively, during bud endodormancy release of tree peony under prolonged chilling treatment. An Illumina RNA sequencing study was performed to identify potential genes involved in the bud endodormancy regulation in tree peony. Correlation matrix, principal component, and interaction network analyses identified a downregulated MYB transcription factor gene, PsMYB306, the expression of which positively correlated with 9-CIS-EPOXYCAROTENOID DIOXYGENASE 3 (PsNCED3) expression. Protein modeling analysis revealed 4 residues within the R2R3 domain of PsMYB306 to possess DNA binding capability. Transcription of PsMYB306 was increased by ABA treatment. Overexpression of PsMYB306 in petunia (Petunia hybrida) inhibited seed germination and plant growth, concomitant with elevated ABA and decreased GA contents. Silencing of PsMYB306 accelerated cold-triggered tree peony bud burst and influenced the production of ABA and GAs and the expression of their biosynthetic genes. ABA application reduced bud dormancy release and transcription of ENT-KAURENOIC ACID OXIDASE 1 (PsKAO1), GA20-OXIDASE 1 (PsGA20ox1), and GA3-OXIDASE 1 (PsGA3ox1) associated with GA biosynthesis in PsMYB306-silenced buds. In vivo and in vitro binding assays confirmed that PsMYB306 specifically transactivated the promoter of PsNCED3. Silencing of PsNCED3 also promoted bud break and growth. Altogether, our findings suggest that PsMYB306 negatively modulates cold-induced bud endodormancy release by regulating ABA production.

Epigenetics governs senescence

Petal is one of the most esthetic and essential parts of a flower that fascinates the pollinators to enhance pollination. Petal senescence is a highly controlled and organized natural phenomenon assisted by phytohormones and gene regulation. It is an inelastically programmed event preceding to which petals give rise to color and scent that captivate pollinators, representing a flower's maturity for sexual reproduction. Till today, many genes involved in the petal senescence through genetic as well as epigenetic changes in response to hormones have been identified. In most of the species, petal senescence is controlled by ethylene, whereas others are independent of this hormone. It has also been proved that the increase in the carbohydrate contents like mannitol, inositol and trehalose delayed the senescence in tulips and Gladiolus. An increased sugar content prevents the biosynthesis of EIN3-like mRNA and further upregulates several senescence correlated genes. A wide range of different transcription factors as well as regulators are disparately expressed in ethylene insensitive and ethylene sensitive petal senescence. DcHB30, a downregulating factor, which upon linking physically to DcWRKY75 leads to the upregulation of ethylene promoting petal senescence. Here we describe the role of ethylene in petal senescence through epigenetic changes. Studies show that ethylene causes petal senescence through epigenetic changes. Feng et al. (Plant Physiol 192:546-564, 2023) observed that ARABIDOPSIS HOMOLOG OF TRITHORAX1 (DcATX1) promotes trimethylation of histone 3 (H3) at 4th lysine (H3K4me3) in Carnation. H3K4me3 further stimulates the expression of genes of ethylene biosynthesis and senescence, leading to senescence in Carnation.

5-Aminolevulinic acid improves cold resistance through regulation of SlMYB4/SlMYB88-SlGSTU43 module to scavenge reactive oxygen species in tomato

Cold stress severely affects the growth and quality of tomato. 5-Aminolevulinic acid (ALA) can effectively improve tomato's cold stress tolerance. In this study, a tomato glutathione S-transferase gene, SlGSTU43, was identified. Results showed that ALA strongly induced the expression of SlGSTU43 under cold stress. SlGSTU43-overexpressing lines showed increased resistance to cold stress through an enhanced ability to scavenge reactive oxygen species. On the contrary, slgstu43 mutant lines were sensitive to cold stress, and ALA did not improve their cold stress tolerance. Thus, SlGSTU43 is a key gene in the process of ALA improving tomato cold tolerance. Through yeast library screening, SlMYB4 and SlMYB88 were preliminarily identified as transcription factors that bind to the SlGSTU43 promoter. Electrophoretic mobility shift, yeast one-hybrid, dual luciferase, and chromatin immunoprecipitation assays experiments verified that SlMYB4 and SlMYB88 can bind to the SlGSTU43 promoter. Further experiments showed that SlMYB4 and SlMYB88 are involved in the process of ALA-improving tomato's cold stress tolerance and they positively regulate the expression of SlGSTU43. The findings provide new insights into the mechanism by which ALA improves cold stress tolerance. SlGSTU43, as a valuable gene, could be added to the cold-responsive gene repository. Subsequently, it could be used in genetic engineering to enhance the cold tolerance of tomato.

Genome-Wide Identification and Expression Pattern of MYB Family Transcription Factors in Erianthus fulvus

MYB family genes have many functions and are widely involved in plant abiotic-stress responses. Erianthus fulvus is an important donor material for stress-resistance genes in sugarcane breeding. However, the MYB family genes in E. fulvus have not been systematically investigated. In this study, 133 EfMYB genes, including 48 Ef1R-MYB, 84 EfR2R3-MYB and 1 Ef3R-MYB genes, were identified in the E. fulvus genome. Among them, the EfR2R3-MYB genes were classified into 20 subgroups. In addition, these EfMYB genes were unevenly distributed across 10 chromosomes. A total of 4 pairs of tandemly duplicated EfMYB genes and 21 pairs of segmentally duplicated EfMYB genes were identified in the E. fulvus genome. Protein-interaction analysis predicted that 24 EfMYB proteins had potential interactions with 14 other family proteins. The EfMYB promoter mainly contains cis-acting elements related to the hormone response, stress response, and light response. Expression analysis showed that EfMYB39, EfMYB84, and EfMYB124 could be significantly induced using low-temperature stress. EfMYB30, EfMYB70, EfMYB81, and EfMYB101 responded positively to drought stress. ABA treatment significantly induced EfMYB1, EfMYB30, EfMYB39, EfMYB84, and EfMYB130. All nine genes were induced using MeJA treatment. These results provide comprehensive information on EfMYB genes and can serve as a reference for further studies of gene function.

Arabidopsis MYB21 Negatively Regulates KTN1 to Fine-Tune the Filament Elongation

The growth process of the stamen filament is crucial for plant reproduction. However, the molecular mechanisms underlying the regulation of filament growth remain largely unclear. Our study has identified that MYB21 is involved in the regulation of filament growth in Arabidopsis. In comparison to the wild type, the cell length of the filaments is notably reduced in the myb21 mutant. Moreover, we found that KTN1, which encodes a microtubule-severing enzyme, is significantly upregulated in the myb21 mutant. Additionally, yeast one-hybrid assays demonstrated that MYB21 can bind to the promoter region of KTN1, suggesting that MYB21 might directly regulate the expression of KTN1. Finally, transcriptional activity experiments showed that MYB21 is capable of suppressing the driving activity of the KTN1 promoter. This study indicates that the MYB21-KTN1 module may play a precise regulatory role in the growth of Arabidopsis filament cells.

Petal size is controlled by the MYB73/TPL/HDA19-miR159-CKX6 module regulating cytokinin catabolism in Rosa hybrida

The size of plant lateral organs is determined by well-coordinated cell proliferation and cell expansion. Here, we report that miR159, an evolutionarily conserved microRNA, plays an essential role in regulating cell division in rose (Rosa hybrida) petals by modulating cytokinin catabolism. We uncover that Cytokinin Oxidase/Dehydrogenase6 (CKX6) is a target of miR159 in petals. Knocking down miR159 levels results in the accumulation of CKX6 transcripts and earlier cytokinin clearance, leading to a shortened cell division period and smaller petals. Conversely, knocking down CKX6 causes cytokinin accumulation and a prolonged developmental cell division period, mimicking the effects of exogenous cytokinin application. MYB73, a R2R3-type MYB transcription repressor, recruits a co-repressor (TOPLESS) and a histone deacetylase (HDA19) to form a suppression complex, which regulates MIR159 expression by modulating histone H3 lysine 9 acetylation levels at the MIR159 promoter. Our work sheds light on mechanisms for ensuring the correct timing of the exit from the cell division phase and thus organ size regulation by controlling cytokinin catabolism.

Evolution of ethylene as an abiotic stress hormone in streptophytes

All land plants modulate their growth and physiology through intricate signaling cascades. The majority of these are at least modulated-and often triggered-by phytohormones. Over the past decade, it has become apparent that some phytohormones have an evolutionary origin that runs deeper than plant terrestrialization-many emerged in the streptophyte algal progenitors of land plants. Ethylene is such a case. Here we synthesize the current knowledge on the evolution of the phytohormone ethylene and speculate about its deeply conserved role in adjusting stress responses of streptophytes for more than half a billion years of evolution.

The CALCINEURIN B-LIKE 4/CBL-INTERACTING PROTEIN 3 module degrades repressor JAZ5 during rose petal senescence

Flower senescence is genetically regulated and developmentally controlled. The phytohormone ethylene induces flower senescence in rose (Rosa hybrida), but the underlying signaling network is not well understood. Given that calcium regulates senescence in animals and plants, we explored the role of calcium in petal senescence. Here, we report that the expression of calcineurin B-like protein 4 (RhCBL4), which encodes a calcium receptor, is induced by senescence and ethylene signaling in rose petals. RhCBL4 interacts with CBL-interacting protein kinase 3 (RhCIPK3), and both positively regulate petal senescence. Furthermore, we determined that RhCIPK3 interacts with the jasmonic acid response repressor jasmonate ZIM-domain 5 (RhJAZ5). RhCIPK3 phosphorylates RhJAZ5 and promotes its degradation in the presence of ethylene. Our results reveal that the RhCBL4-RhCIPK3-RhJAZ5 module mediates ethylene-regulated petal senescence. These findings provide insights into flower senescence, which may facilitate innovations in postharvest technology for extending rose flower longevity.

The ARF2-MYB6 module mediates auxin-regulated petal expansion in rose

In cut rose (Rosa hybrida), the flower-opening process is closely associated with vase life. Auxin induces the expression of transcription factor genes that function in petal growth via cell expansion. However, the molecular mechanisms underlying the auxin effect during flower opening are not well understood. Here, we identified the auxin-inducible transcription factor gene RhMYB6, whose expression level is high during the early stages of flower opening. Silencing of RhMYB6 delayed flower opening by controlling petal cell expansion through down-regulation of cell expansion-related genes. Furthermore, we demonstrated that the auxin response factor RhARF2 directly interacts with the promoter of RhMYB6 and represses its transcription. Silencing of RhARF2 resulted in larger petal size and delayed petal movement. We also showed that the expression of genes related to ethylene and petal movement showed substantial differences in RhARF2-silenced petals. Our results indicate that auxin-regulated RhARF2 is a critical player that controls flower opening by governing RhMYB6 expression and mediating the crosstalk between auxin and ethylene signaling.

Comparative transcriptome analysis reveals the involvement of an MYB transcriptional activator, SmMYB108, in anther dehiscence in eggplant

Male sterility is a highly attractive agronomic trait as it effectively prevents self-fertilization and facilitates the production of high-quality hybrid seeds in plants. Timely release of mature pollen following anther dehiscence is essential for stamen development in flowering plants. Although several theories have been proposed regarding this, the specific mechanism of anther development in eggplant remains elusive. In this study, we selected an R2R3-MYB transcription factor gene, SmMYB108, that encodes a protein localized primarily in the nucleus by comparing the transcriptomics of different floral bud developmental stages of the eggplant fertile line, F142. Quantitative reverse transcription polymerase chain reaction revealed that SmMYB108 was preferentially expressed in flowers, and its expression increased significantly on the day of flowering. Overexpression of SmMYB108 in tobacco caused anther dehiscence. In addition, we found that SmMYB108 primarily functions as a transcriptional activator via C-terminal activation (amino acid 262-317). Yeast one-hybrid and dual-luciferase reporter assays revealed that genes (SmMYB21, SmARF6, and SmARF8) related to anther development targeted the SmMYB108 promoter. Overall, our results provide insights into the molecular mechanisms involved in the regulation of anther development by SmMYB108.

Polygalacturonase gene family analysis identifies FcPG12 as a key player in fig (Ficus carica L.) fruit softening

BACKGROUND

The fig (Ficus carica L.) tree has high economic value. However, its fruit have a short shelf life due to rapid softening. Polygalacturonases (PGs) are essential hydrolases, responsible for the pectin degradation that plays a key role in fruit softening. However, fig PG genes and their regulators have not yet been characterized.

RESULTS

In this study, 43 FcPGs were identified in the fig genome. They were non-uniformly distributed on 13 chromosomes, and tandem repeat PG gene clusters were found on chromosomes 4 and 5. Ka/Ks calculation and collinear analysis indicated negative selection as the main driver of FcPG family expansion. Fourteen FcPGs were found expressed in fig fruit with FPKM values > 10, of which seven were positively correlated, and three, negatively correlated with fruit softening. Eleven FcPGs were upregulated and two downregulated in response to ethephon treatment. FcPG12, a member of the tandem repeat cluster on chromosome 4, was selected for further analyses due to its sharp increment in transcript abundance during fruit softening and its response to ethephon treatment. Transient overexpression of FcPG12 led to decreased fig fruit firmness and increased PG enzyme activity in the tissue. Two ethylene response factor (ERF)-binding GCC-box sites were found on the FcPG12 promoter. Yeast one-hybrid and dual luciferase assays showed that FcERF5 binds directly to the FcPG12 promoter and upregulates its expression. Transient overexpression of FcERF5 upregulated FcPG12 expression, thereby increasing PG activity and fruit softening.

CONCLUSIONS

Our study identified FcPG12 as a key PG gene in fig fruit softening, and its direct positive regulation by FcERF5. The results provide new information on the molecular regulation of fig fruit softening.

BASIC PENTACYSTEINE2 fine-tunes corm dormancy release in Gladiolus

Bud dormancy is an important trait in geophytes that largely affects their flowering process and vegetative growth after dormancy release. Compared with seed dormancy, the regulation of bud dormancy is still largely unclear. Abscisic acid (ABA) acts as the predominant hormone that regulates the whole dormancy process. In Gladiolus (Gladiolus hybridus), cold storage promotes corm dormancy release (CDR) by repressing ABA biosynthesis and signaling. However, the mechanisms governing ABA-related processes during CDR via epigenetics are poorly understood. Here, we show that class I BASIC PENTACYSTEINE2, (GhBPC2) directly binds to 9-CIS-EPOXYCAROTENOID DIOXYGENASE (GhNCED) and ABA INSENSITIVE5 (GhABI5) loci and down-regulates their expression to accelerate CDR. During CDR, histone modifications change dramatically at the GhBPC2-binding loci of GhABI5 with an increase in H3K27me3 and a decrease in H3K4me3. GhBPC2 is involved in both H3K27me3 and H3K4me3 and fine-tunes GhABI5 expression by recruiting polycomb repressive complex 2 (PRC2) and the chromatin remodeling factor EARLY BOLTING IN SHORT DAYS (GhEBS). These results show GhBPC2 epigenetically regulates CDR in Gladiolus by mediating GhABI5 expression with PRC2 and GhEBS.

Recent Advances in Research into Jasmonate Biosynthesis and Signaling Pathways in Agricultural Crops and Products

Jasmonates (JAs) are phospholipid-derived hormones that regulate plant development and responses to environmental stress. The synthesis of JAs and the transduction of their signaling pathways are precisely regulated at multiple levels within and outside the nucleus as a result of a combination of genetic and epigenetic regulation. In this review, we focus on recent advances in the regulation of JA biosynthesis and their signaling pathways. The biosynthesis of JAs was found to be regulated with an autocatalytic amplification mechanism via the MYC2 regulation pathway and inhibited by an autonomous braking mechanism via the MYC2-targeting bHLH1 protein to terminate JA signals in a highly ordered manner. The biological functions of JAs mainly include the promotion of fruit ripening at the initial stage via ethylene-dependent and independent ways, the regulation of mature coloring via regulating the degradation of chlorophyll and the metabolism of anthocyanin, and the improvement of aroma components via the regulation of fatty acid and aldehyde alcohol metabolism in agricultural crops. JA signaling pathways also function in the enhancement of biotic and abiotic stress resistance via the regulation of secondary metabolism and the redox system, and they relieve cold damage to crops through improving the stability of the cell membrane. These recently published findings indicate that JAs are an important class of plant hormones necessary for regulating plant growth and development, ripening, and the resistance to stress in agricultural crops and products.

Senescence in dahlia flowers is regulated by a complex interplay between flower age and floret position

Mechanisms regulating flower senescence are not fully understood in any species and are particularly complex in composite flowers. Dahlia (Dahlia pinnata Cav.) florets develop sequentially, hence each composite flower head includes florets of different developmental stages as the whole flower head ages. Moreover, the wide range of available cultivars enables assessment of intraspecific variation. Transcriptomes were compared amongst inner (younger) and outer (older) florets of two flower head ages to assess the effect of floret vs. flower head ageing. More gene expression, including ethylene and cytokinin pathway expression changed between inner and outer florets of older flower heads than between inner florets of younger and older flower heads. Additionally, based on Arabidopsis network analysis, different patterns of co-expressed ethylene response genes were elicited. This suggests that changes occur in young inner florets as the whole flower head ages that are different to ageing florets within a flower head. In some species floral senescence is orchestrated by the plant growth regulator ethylene. However, there is both inter and intra-species variation in its importance. There is a lack of conclusive data regarding ethylene sensitivity in dahlia. Speed of senescence progression, effects of ethylene signalling perturbation, and patterns of ethylene biosynthesis gene expression differed across three dahlia cultivars ('Sylvia', 'Karma Prospero' and 'Onesta') suggesting differences in the role of ethylene in their floral senescence, while effects of exogenous cytokinin were less cultivar-specific.

Genome-wide identification and expression analysis of MIKCC genes in rose provide insight into their effects on flower development

The MIKC^C-type gene family plays important roles in plant growth, development, and tolerance of biotic and abiotic stress, especially during floral organ differentiation. However, there have been no studies of MIKC^C-type genes in rose, and functional differentiation of family members has not been explored. In this study, we identified 42 MIKC^C-type genes in rose, classified the genes into 12 subfamilies, and constructed a phylogenetic tree. We performed expression analysis of these genes, and found that expression patterns correlated with the predicted subfamily, indicating that the features of MIKC^C-type genes were broadly conserved during evolution. Collinear analysis of MIKC^C genes among Rosaceae species confirmed the occurrence of whole genome duplications (WGD) and revealed some species-specific MIKC^C genes. Transcriptome analysis showed that the expression of some MIKC^C-type genes responded to low temperatures (4°C, 24 h) during flower organ differentiation. These conserved, duplicated, and novel expression patterns of MIKC^C-type genes may have facilitated the adaptation of rose to various internal and external environmental changes. The results of this study provide a theoretical basis for future functional analysis of the MIKC^C genes in rose and investigation of the evolutionary pattern of the MIKC^C gene family in the Rosaceae genome.

Protein Kinase RhCIPK6 Promotes Petal Senescence in Response to Ethylene in Rose (Rosa Hybrida)

Cultivated roses have the largest global market share among ornamental crops. Postharvest release of ethylene is the main cause of accelerated senescence and decline in rose flower quality. To understand the molecular mechanism of ethylene-induced rose petal senescence, we analyzed the transcriptome of rose petals during natural senescence as well as with ethylene treatment. A large number of differentially expressed genes (DEGs) were observed between developmental senescence and the ethylene-induced process. We identified 1207 upregulated genes in the ethylene-induced senescence process, including 82 transcription factors and 48 protein kinases. Gene Ontology enrichment analysis showed that ethylene-induced senescence was closely related to stress, dehydration, and redox reactions. We identified a calcineurin B-like protein (CBL) interacting protein kinase (CIPK) family gene in Rosa hybrida, RhCIPK6, that was regulated by age and ethylene induction. Reducing RhCIPK6 expression through virus-induced gene silencing significantly delayed petal senescence, indicating that RhCIPK6 mediates petal senescence. In the RhCIPK6-silenced petals, several senescence associated genes (SAGs) and transcription factor genes were downregulated compared with controls. We also determined that RhCIPK6 directly binds calcineurin B-like protein 3 (RhCBL3). Our work thus offers new insights into the function of CIPKs in petal senescence and provides a genetic resource for extending rose vase life.

The chromosome-level genome of Gypsophila paniculata reveals the molecular mechanism of floral development and ethylene insensitivity

Gypsophila paniculata, belonging to the Caryophyllaceae of the Caryophyllales, is one of the most famous worldwide cut flowers. It is commonly used as dried flowers, whereas the underlying mechanism of flower senescence has not yet been addressed. Here, we present a chromosome-scale genome assembly for G. paniculata with a total size of 749.58 Mb. Whole-genome duplication signatures unveil two major duplication events in its evolutionary history: an ancient one occurring before the divergence of Caryophyllaceae and a more recent one shared with Dianthus caryophyllus. The integrative analyses combining genomic and transcriptomic data reveal the mechanisms regulating floral development and ethylene response of G. paniculata. The reduction of AGAMOUS expression probably caused by sequence polymorphism and the mutation in miR172 binding site of PETALOSA are associated with the double flower formation in G. paniculata. The low expression of ETHYLENE RESPONSE SENSOR (ERS) and the reduction of downstream ETHYLENE RESPONSE FACTOR (ERF) gene copy number collectively lead to the ethylene insensitivity of G. paniculata, affecting flower senescence and making it capable of making dried flowers. This study provides a cornerstone for understanding the underlying principles governing floral development and flower senescence, which could accelerate the molecular breeding of the Caryophyllaceae species.

Transcriptome Level Analysis of Genes of Exogenous Ethylene Applied under Phosphorus Stress in Chinese Fir

Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) is a widely grown gymnosperm in China. Phosphorus (P) is an indispensable nutrient for the growth of Chinese fir. Inorganic phosphate (Pi) deficiency exists in soils of many Chinese fir planting area regions, and the trees themselves have limited efficiency in utilizing P from the soil. Ethylene is important in regulation responses to nutrient deficiencies. However, little is known about how ethylene signals participate in Pi stress in Chinese fir. A total of six different treatments were performed to reveal the transcript levels of Chinese fir under Pi, ethephon (an ethylene-releasing compound), and CoCl₂ (cobalt chloride, an ethylene biosynthesis inhibitor) treatments. We assembled a full-length reference transcriptome containing 22,243 unigenes as a reference for UMI RNA-seq (Digital RNA-seq). There were 586 Differentially Expressed Genes (DEGs) in the Pi starvation (NP) group, while DEGs from additional ethephon or CoCl₂ in NP were 708 and 292, respectively. Among the DEGs in each treatment, there were 83 TFs in these treatment groups. MYB (v-myb avian myeloblastosis viral oncogene homolog) family was the most abundant transcription factors (TFs). Three ERF (Ethylene response factor) family genes were identified when only ethylene content was imposed as a variable. Enrichment analysis indicated that the ascorbate and aldarate metabolism pathway plays a key role in resistance to Pi deficiency. This study provides insights for further elucidating the regulatory mechanism of Pi deficiency in Chinese fir.

Natural variation MeMYB108 associated with tolerance to stress-induced leaf abscission linked to enhanced protection against reactive oxygen species in cassava

Natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers led to alleviate leaf abscission under drought in cassava. The reactive oxygen scavengers play important roles in regulating the cassava (Manihot esculenta Crantz) leaf abscission induced by stresses. To date, the relationship between natural variations of MYB genes and reactive oxygen scavengers under drought in cassava genotypes remains unclear. Here, we reported the transcription factor MeMYB108 played an important role in regulating leaf abscission exposed to drought in cassava. The expression levels of MeMYB108 in abscission zones of cassava leaf pulvinus were higher in cassava genotype SC124, which were less easy to shed leaves under stress than cassava genotype SC8 when the leaf abscission induced by the same drought condition. Compared with wild type and interference expression plants, overexpression of MeMYB108 significantly reduced the drought-induced leaf abscission rate under drought. The consecutively 2-year analysis of reactive oxygen scavengers showed significant differences among different cassava genotypes under drought-induced leaf abscission, indicating the relevance between reactive oxygen scavengers and leaf abscission. Correlation analysis revealed the natural variation of the MeMYB108 exon was associated with reactive oxygen scavengers during drought-induced leaf abscission. Association analysis between pairwise LD of DNA polymorphism indicated the MeMYB108 allele enhanced the tolerance of cassava to drought-induced leaf abscission. Complementation transgenic lines containing the elite allele of MeMYB108 ^SC124 decreased the leaf abscission rate induced by drought conditions, demonstrating natural variation in MeMYB108 contributed to leaf abscission tolerance induced by drought in cassava. Further studies showed MeMYB108 played an active role in the tolerance of cassava to drought-induced leaf abscission by inducing scavenging of reactive oxygen species.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Thidiazuron Promotes Leaf Abscission by Regulating the Crosstalk Complexities between Ethylene, Auxin, and Cytokinin in Cotton

Thidiazuron (TDZ) is widely used as a defoliant to induce leaf abscission in cotton. However, the underlying molecular mechanism is still unclear. In this study, RNA-seq and enzyme-linked immunosorbent assays (ELISA) were performed to reveal the dynamic transcriptome profiling and the change of endogenous phytohormones upon TDZ treatment in leaf, petiole, and abscission zone (AZ). We found that TDZ induced the gene expression of ethylene biosynthesis and signal, and promoted ethylene accumulation earlier in leaf than that in AZ. While TDZ down-regulated indole-3-acetic acid (IAA) biosynthesis genes mainly in leaf and IAA signal and transport genes. Furthermore, the IAA content reduced more sharply in the leaf than that in AZ to change the auxin gradient for abscission. TDZ suppressed CTK biosynthesis genes and induced CTK metabolic genes to reduce the IPA accumulation for the reduction of ethylene sensitivity. Furthermore, TDZ regulated the gene expression of abscisic acid (ABA) biosynthesis and signal and induced ABA accumulation between 12-48 h, which could up-regulate ABA response factor genes and inhibit IAA transporter genes. Our data suggest that TDZ orchestrates metabolism and signal of ethylene, auxin, and cytokinin, and also the transport of auxin in leaf, petiole, and AZ, to control leaf abscission.

A combined transcriptomic and proteomic analysis of chrysanthemum provides new insights into petal senescence

Numerous transcription factor genes and methylation-related genes were differentially expressed in senescent petals compared with control petals. Studying petal senescence is crucial for extending the postharvest longevity of cut flowers, but petal senescence remains relatively unexplored compared to well-studied leaf senescence. In this study, a combined transcriptomic and proteomic analysis of senescent (22 days after cutting) and control (0 day after cutting) petals was performed to investigate the molecular processes underlying petal senescence of chrysanthemum (Chrysanthemum morifolium Ramat.), an important cut flower crop worldwide. A total of 11,324 differentially expressed genes (DEGs), including 4888 up-regulated and 6436 down-regulated genes, and 403 differentially expressed proteins (DEPs), including 210 up-regulated and 193 down-regulated proteins, were identified at transcript and protein levels, respectively. A cross-comparison of transcriptomic and proteomic data identified 257 consistent DEGs/DEPs, including 122 up-regulated and 135 down-regulated DEGs/DEPs. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that "cutin, suberine and wax biosynthesis" is a main pathway for both DEGs and DEPs, especially for down-regulated DEGs/DEPs. Functional analysis indicated that chrysanthemum genes mainly encoding putative cytochrome P450s, non-specific lipid-transfer proteins, subtilisin-like proteases, AAA-ATPases, proteins essential for cuticular wax biosynthesis, and proteins in hormone signal transduction or ubiquitination were differentially expressed at both transcript and protein levels. In addition, numerous transcription factor genes and methylation-related genes were also differentially expressed, inferring an involvement of transcriptional and epigenetic regulation in petal senescence. These results provide a valuable resource of studying chrysanthemum senescence and significant insights into petal senescence.

Transcriptome and Metabolome Analyses Provide Insights into the Stomium Degeneration Mechanism in Lily

Lily (Lilium spp.) is a widely cultivated horticultural crop that has high ornamental and commercial value but also the serious problem of pollen pollution. However, mechanisms of anther dehiscence in lily remain largely unknown. In this study, the morphological characteristics of the stomium zone (SZ) from different developmental stages of ‘Siberia’ lily anthers were investigated. In addition, transcriptomic and metabolomic data were analyzed to identify the differentially expressed genes (DEGs) and secondary metabolites involved in stomium degeneration. According to morphological observations, SZ lysis occurred when flower buds were 6–8 cm in length and was completed in 9 cm. Transcriptomic analysis identified the genes involved in SZ degeneration, including those associated with hormone signal transduction, cell structure, reactive oxygen species (ROS), and transcription factors. A weighted co-expression network showed strong correlations between transcription factors. In addition, TUNEL (TdT-mediated dUTP nick-end labeling) assays showed that programmed cell death was important during anther SZ degeneration. Jasmonates might also have key roles in anther dehiscence by affecting the expression of the genes involved in pectin lysis, water transport, and cysteine protease. Collectively, the results of this study improve our understanding of anther dehiscence in lily and provide a data platform from which the molecular mechanisms of SZ degeneration can be revealed.

RhWRKY33 Positively Regulates Onset of Floral Senescence by Responding to Wounding- and Ethylene-Signaling in Rose Plants

Rose plants are one of the most important horticultural crops, whose commercial value mainly depends on long-distance transportation, and wounding and ethylene are the main factors leading to their quality decline and accelerated senescence in the process. However, underlying molecular mechanisms of crosstalk between wounding and ethylene in the regulation of flower senescence remain poorly understood. In relation to this, transcriptome analysis was performed on rose flowers subjected to various treatments, including control, wounding, ethylene, and wounding- and ethylene- (EW) dual treatment. A large number of differentially expressed genes (DEGs) were identified, ranging from 2,442 between the ethylene- and control-treated groups to 4,055 between the EW- and control-treated groups. Using weighted gene co-expression network analysis (WGCNA), we identified a hub gene RhWRKY33 (rchiobhmchr5g0071811), accumulated in the nucleus, where it may function as a transcription factor. Moreover, quantitative reverse transcription PCR (RT-qPCR) results showed that the expression of RhWRKY33 was higher in the wounding-, ethylene, and EW-treated petals than in the control-treated petals. We also functionally characterized the RhWRKY33 gene through virus-induced gene silencing (VIGS). The silencing of RhWRKY33 significantly delayed the senescence process in the different treatments (control, wounding, ethylene, and EW). Meanwhile, we found that the effect of RhWRKY33-silenced petals under ethylene and EW dual-treatment were stronger than those under wounding treatment in delaying the petal senescence process, implying that RhWRKY33 is closely involved with ethylene and wounding mediated petal senescence. Overall, the results indicate that RhWRKY33 positively regulates the onset of floral senescence mediated by both ethylene and wounding signaling, but relies heavily on ethylene signaling.

Genome and transcriptome-based characterization of high energy carbon-ion beam irradiation induced delayed flower senescence mutant in Lotus japonicus

BACKGROUND

Flower longevity is closely related to pollen dispersal and reproductive success in all plants, as well as the commercial value of ornamental plants. Mutants that display variation in flower longevity are useful tools for understanding the mechanisms underlying this trait. Heavy-ion beam irradiation has great potential to improve flower shapes and colors; however, few studies are available on the mutation of flower senescence in leguminous plants.

RESULTS

A mutant (C416) exhibiting blossom duration eight times longer than that of the wild type (WT) was isolated in Lotus japonicus derived from carbon ion beam irradiation. Genetic assays supported that the delayed flower senescence of C416 was a dominant trait controlled by a single gene, which was located between 4,616,611 Mb and 5,331,876 Mb on chromosome III. By using a sorting strategy of multi-sample parallel genome sequencing, candidate genes were narrowed to the gene CUFF.40834, which exhibited high identity to ethylene receptor 1 in other model plants. A physiological assay demonstrated that C416 was insensitive to ethylene precursor. Furthermore, the dynamic changes of phytohormone regulatory network in petals at different developmental stages was compared by using RNA-seq. In brief, the ethylene, jasmonic acid (JA), and salicylic acid (SA) signaling pathways were negatively regulated in C416, whereas the brassinosteroid (BR) and cytokinin signaling pathways were positively regulated, and auxin exhibited dual effects on flower senescence in Lotus japonicus. The abscisic acid (ABA) signaling pathway is positively regulated in C416.

CONCLUSION

So far, C416 might be the first reported mutant carrying a mutation in an endogenous ethylene-related gene in Lotus japonicus, rather than through the introduction of exogenous genes by transgenic techniques. A schematic of the flower senescence of Lotus japonicus from the perspective of the phytohormone regulatory network was provided based on transcriptome profiling of petals at different developmental stages. This study is informative for elucidating the molecular mechanism of delayed flower senescence in C416, and lays a foundation for candidate flower senescence gene identification in Lotus japonicus. It also provides another perspective for the improvement of flower longevity in legume plants by heavy-ion beam.

Transcriptome Analysis Reveals the Senescence Process Controlling the Flower Opening and Closure Rhythm in the Waterlilies (Nymphaea L.)

Most waterlily flowers open at dawn and close after noon usually for three to four days, and thereafter wilt. The short lifespan of flowers restricts the development of the flower postharvest industry. The termination of flower movements is a key event during flower aging process. However, it is still unclear when the senescence process initiates and how it terminates the movement rhythm. In this study, we observed that the opening diameter of flowers was the smallest on the fourth (last) flowering day. Subsequent transcriptome profiles generated from petals at different flowering stages showed that the multiple signaling pathways were activated at the last closure stage (Time 3, T3) of the flowers, including Ca²⁺, reactive oxygen species and far red light signaling pathways, as well as auxin, ethylene and jasmonic acid signaling pathways. Moreover, In terms of cell metabolism regulation, the genes related to hydrolase (protease, phospholipase, nuclease) were upregulated at T3 stage, indicating that petals entered the senescence stage at that time; and the genes related to water transport and cell wall modification were also differentially regulated at T3 stage, which would affect the ability of cell expand and contract, and eventually lead to petal not open after the fourth day. Collectively, our data provided a new insight into the termination of flower opening in the waterlilies, and a global understanding of the senescence process of those opening-closure rhythm flowers.

SMRT and Illumina RNA Sequencing and Characterization of a Key NAC Gene LoNAC29 during the Flower Senescence in Lilium oriental 'Siberia'

Lily (Lilium spp.) is an important cut flower around the world. Flower senescence in lilies is characterized by the wilting and abscission of tepals, which results in a decrease in flower quality and huge economic loss. However, the mechanism underlying flower senescence in lilies is largely unknown. In this study, single-molecule, real-time (SMRT) and Illumina sequencing were carried out in L. oriental ‘Siberia’. Sequencing yielded 73,218 non-redundant transcripts, with an N50 of 3792 bp. These data were further integrated with three published transcriptomes through cogent analysis, which yielded 62,960 transcripts, with an increase in N50 of 3935 bp. Analysis of differentially expressed genes showed that 319 transcription factors were highly upregulated during flower senescence. The expression of twelve NAC genes and eleven senescence-associated genes (SAGs) showed that LoNAC29 and LoSAG39 were highly expressed in senescent flowers. Transient overexpression of LoNAC29 and LoSAG39 in tepals of lily notably accelerated flower senescence, and the promoter activity of LoSAG39 was strongly induced by LoNAC29. This work supported new evidence for the molecular mechanism of flower senescence and provided better sequence data for further study in lilies.

Mediation of JA signalling in glandular trichomes by the woolly/SlMYC1 regulatory module improves pest resistance in tomato

Almost all plants form trichomes, which protect them against insect herbivores by forming a physical barrier and releasing chemical repellents. Glandular trichomes produce a variety of specialized defensive metabolites, including volatile terpenes. Previous studies have shown that the defence hormone jasmonic acid (JA) affects trichome development and induces terpene synthases (TPSs) but the underlying molecular mechanisms remain unclear. Here, we characterized a loss-of-function allele of the HD-ZIP IV transcription factor woolly (wo) and analysed its role in mediating JA signalling in tomato. We showed that knockout of wo led to extensive trichome defects, including structural and functional changes in type VI glandular trichomes, and a dramatic reduction in terpene levels. We further found that wo directly binds to TPS gene promoters to recruit SlMYC1, a JA signalling modulator, and that together these transcription factors promote terpene biosynthesis in tomato trichomes. The wo/SlMYC1 regulatory module is inhibited by SlJAZ2 through a competitive binding mechanism, resulting in a fine-tuned JA response in tomato trichomes. Enhanced expression of SlMYC1 substantially increased terpene levels and improved tomato resistance to spider mites. Interestingly, we also found that SlMYC1 plays an additional role in glandular cell division and expansion in type VI trichomes, independent of JA. Together, our results reveal a novel, JA-mediated regulatory mechanism that promotes insect resistance in tomato.

The regulatory framework of developmentally programmed cell death in floral organs: A review

Developmentally programmed cell death (dPCD) is a tightly controlled biological process. In recent years, vital roles of dPCD on regulating floral organ growth and development have been reported. It is well known that flower is an essential organ for reproduction and a turning point of plants' life cycle. Hence, uncovering the complex molecular networks which regulates dPCD processes in floral organs is utmost important. So far, our understanding of dPCD on floral organ growth and development is just starting. Herein, we summarize the important factors that involved in the tapetal degeneration, pollen tube rupture, receptive synergid cell death, nucellar degradation, and antipodal cell degradation. Meanwhile, the known factors that involved in transmitting tract formation and self-incompatibility-induced PCD were also introduced. Furthermore, the genes that associated with anther dehiscence and petal senescence and abscission were reviewed as well. The functions of various types of factors involved in floral dPCD processes are highlighted principally. The regulatory panorama described here can provide us some insights about flower-specific dPCD process.

An R2R3-MYB Transcription Factor RmMYB108 Responds to Chilling Stress of Rosa multiflora and Conferred Cold Tolerance of Arabidopsis

A sudden cooling in the early spring or late autumn negatively impacts the plant growth and development. Although a number of studies have characterized the role of the transcription factors (TFs) of plant R2R3-myeloblastosis (R2R3-MYB) in response to biotic and abiotic stress, plant growth, and primary and specific metabolisms, much less is known about their role in Rosa multiflora under chilling stress. In the present study, RmMYB108, which encodes a nuclear-localized R2R3-MYB TF with a self-activation activity, was identified based on the earlier published RNA-seq data of R. multiflora plants exposed to short-term low-temperature stress and also on the results of prediction of the gene function referring Arabidopsis. The RmMYB108 gene was induced by stress due to chilling, salt, and drought and was expressed in higher levels in the roots than in the leaves. The heterologous expression of RmMYB108 in Arabidopsis thaliana significantly enhanced the tolerance of transgenic plants to freezing, water deficit, and high salinity, enabling higher survival and growth rates, earlier flowering and silique formation, and better seed quantity and quality compared with the wild-type (WT) plants. When exposed to a continuous low-temperature stress at 4°C, transgenic Arabidopsis lines-overexpressing RmMYB108 showed higher activities of superoxide dismutase and peroxidase, lower relative conductivity, and lower malondialdehyde content than the WT. Moreover, the initial fluorescence (F _o) and maximum photosynthetic efficiency of photosystem II (F _v/F _m) changed more dramatically in the WT than in transgenic plants. Furthermore, the expression levels of cold-related genes involved in the ICE1 (Inducer of CBF expression 1)-CBFs (C-repeat binding factors)-CORs (Cold regulated genes) cascade were higher in the overexpression lines than in the WT. These results suggest that RmMYB108 was positively involved in the tolerance responses when R. multiflora was exposed to challenges against cold, freeze, salt, or drought and improved the cold tolerance of transgenic Arabidopsis by reducing plant damage and promoting plant growth.

The RhHB1/RhLOX4 module affects the dehydration tolerance of rose flowers (Rosa hybrida) by fine-tuning jasmonic acid levels

Phytohormones are key factors in plant responsiveness to abiotic and biotic stresses, and maintaining hormone homeostasis is critically important during stress responses. Cut rose (Rosa hybrida) flowers experience dehydration stress during postharvest handling, and jasmonic acid (JA) levels change as a result of this stress. However, how JA is involved in dehydration tolerance remains unclear. We investigated the functions of the JA- and dehydration-induced RhHB1 gene, which encodes a homeodomain-leucine zipper I γ-clade transcription factor, in rose flowers. Silencing RhHB1 decreased petal dehydration tolerance and resulted in a persistent increase in JA-Ile content and reduced dehydration tolerance. An elevated JA-Ile level had a detrimental effect on rose petal dehydration tolerance. RhHB1 was shown to lower the transient induction of JA-Ile accumulation in response to dehydration. In addition to transcriptomic data, we obtained evidence that RhHB1 suppresses the expression of the lipoxygenase 4 (RhLOX4) gene by directly binding to its promoter both in vivo and in vitro. We propose that increased JA-Ile levels weaken the capacity for osmotic adjustment in petal cells, resulting in reduced dehydration tolerance. In conclusion, a JA feedback loop mediated by an RhHB1/RhLOX4 regulatory module provides dehydration tolerance by fine-tuning bioactive JA levels in dehydrated flowers.

Molecular Evidences for the Interactions of Auxin, Gibberellin, and Cytokinin in Bent Peduncle Phenomenon in Rose (Rosa sp.)

In roses (Rosa sp.), peduncle morphology is an important ornamental feature. The common physiological abnormality known as the bent peduncle phenomenon (BPP) seriously decreases the quality of rose flowers and thus the commercial value. Because the molecular mechanisms underlying this condition are poorly understood, we analysed the transcriptional profiles and cellular structures of bent rose peduncles. Numerous differentially expressed genes involved in the auxin, cytokinin, and gibberellin signaling pathways were shown to be associated with bent peduncle. Paraffin sections showed that the cell number on the upper sides of bent peduncles was increased, while the cells on the lower sides were larger than those in normal peduncles. We also investigated the large, deformed sepals that usually accompany BPP and found increased expression level of some auxin-responsive genes and decreased expression level of genes that are involved in cytokinin and gibberellin synthesis in these sepals. Furthermore, removal of the deformed sepals partially relieved BPP. In summary, our findings suggest that auxin, cytokinin, and gibberellin all influence the development of BPP by regulating cell division and expansion. To effectively reduce BPP in roses, more efforts need to be devoted to the molecular regulation of gibberellins and cytokinins in addition to that of auxin.

Erratum: Author Correction: RhMYB108, an R2R3-MYB transcription factor, is involved in ethylene- and JA-induced petal senescence in rose plants

[This corrects the article DOI: 10.1038/s41438-019-0221-8.].