Integrated multi-omic data and analyses reveal the pathways underlying key ornamental traits in carnation flowers

A high-quality haplotype genome of Michelia alba DC reveals differences in methylation patterns and flower characteristics

Michelia alba DC is a highly valuable ornamental plant of the Magnoliaceae family. This evergreen tropical tree commonly grows in Southeast Asia and is adored for its delightful fragrance. Our study assembled the M. alba haplotype genome MC and MM by utilizing Nanopore ultralong reads, Pacbio Hifi long reads and parental second-generation data. Moreover, the first methylation map of Magnoliaceae was constructed based on the methylation site data obtained using Nanopore data. Metabolomic datasets were generated from the flowers of three different species to assess variations in pigment and volatile compound accumulation. Finally, transcriptome data were generated to link genomic, methylation, and morphological patterns to reveal the reasons underlying the differences between M. alba and its parental lines in petal color, flower shape, and fragrance. We found that the AP1 and AP2 genes are crucial in M. alba petal formation, while the 4CL, PAL, and C4H genes control petal color. The data generated in this study serve as a foundation for future physiological and biochemical research on M. alba, facilitate the targeted improvement of M. alba varieties, and offer a theoretical basis for molecular research on Michelia L.

Haplotype-resolved genome assembly of the diploid Rosa chinensis provides insight into the mechanisms underlying key ornamental traits

Roses are consistently ranked at the forefront in cut flower production. Increasing demands of market and changing climate conditions have resulted in the need to further improve the diversity and quality of traits. However, frequent hybridization leads to highly heterozygous nature, including the allelic variants. Therefore, the absence of comprehensive genomic information leads to them making it challenging to molecular breeding. Here, two haplotype-resolved chromosome genomes for Rosa chinensis 'Chilong Hanzhu' (2n = 14) which is high heterozygous diploid old Chinese rose are generated. An amount of genetic variation (1,605,616 SNPs, 209,575 indels) is identified. 13,971 allelic genes show differential expression patterns between two haplotypes. Importantly, these differences hold valuable insights into regulatory mechanisms of traits. RcMYB114b can influence cyanidin-3-glucoside accumulation and the allelic variation in its promoter leads to differences in promoter activity, which as a factor control petal color. Moreover, gene family expansion may contribute to the abundance of terpenes in floral scents. Additionally, RcANT1, RcDA1, RcAG1 and RcSVP1 genes are involved in regulation of petal number and size under heat stress treatment. This study provides a foundation for molecular breeding to improve important characteristics of roses.

Integration of Metabolomic and Transcriptomic Analyses Reveals the Molecular Mechanisms of Flower Color Formation in Prunus mume

Flower color is an important trait that affects the economic value of Prunus mume, a famous ornamental plant in the Rosaceae family. P. mume with purple-red flowers is uniquely charming and highly favored in landscape applications. However, little is known about its flower coloring mechanism, which stands as a critical obstacle on the path to innovative breeding for P. mume flower color. In this study, transcriptomic and targeted metabolomic analyses of purple-red P. mume and white P. mume were performed to elucidate the mechanism of flower color formation. In addition, the expression patterns of key genes were analyzed using an RT-qPCR experiment. The results showed that the differential metabolites were significantly enriched in the flavonoid synthesis pathway. A total of 14 anthocyanins emerged as the pivotal metabolites responsible for the differences in flower color between the two P. mume cultivars, comprising seven cyanidin derivatives, five pelargonium derivatives, and two paeoniflorin derivatives. Moreover, the results clarified that the metabolic pathway determining flower color in purple-red P. mume encompasses two distinct branches: cyanidin and pelargonidin, excluding the delphinidin branch. Additionally, through the integrated analysis of transcriptomic and metabolomic data, we identified 18 key genes responsible for anthocyanin regulation, thereby constructing the gene regulatory network for P. mume anthocyanin synthesis. Among them, ten genes (PmCHI, PmGT2, PmGT5, PmGST3, PmMYB17, PmMYB22, PmMYB23, PmbHLH4, PmbHLH10, and PmbHLH20) related to anthocyanin synthesis were significantly positively correlated with anthocyanin contents, indicating that they may be the key contributors to anthocyanin accumulation. Our investigation contributes a novel perspective to understanding the mechanisms responsible for flower color formation in P. mume. The findings of this study introduce novel strategies for molecular design breeding aimed at manipulating flower color in P. mume.

The haplotype-resolved telomere-to-telomere carnation (Dianthus caryophyllus) genome reveals the correlation between genome architecture and gene expression

Carnation (Dianthus caryophyllus) is one of the most valuable commercial flowers, due to its richness of color and form, and its excellent storage and vase life. The diverse demands of the market require faster breeding in carnations. A full understanding of carnations is therefore required to guide the direction of breeding. Hence, we assembled the haplotype-resolved gap-free carnation genome of the variety 'Baltico', which is the most common white standard variety worldwide. Based on high-depth HiFi, ultra-long nanopore, and Hi-C sequencing data, we assembled the telomere-to-telomere (T2T) genome to be 564 479 117 and 568 266 215 bp for the two haplotypes Hap1 and Hap2, respectively. This T2T genome exhibited great improvement in genome assembly and annotation results compared with the former version. The improvements were seen when different approaches to evaluation were used. Our T2T genome first informs the analysis of the telomere and centromere region, enabling us to speculate about specific centromere characteristics that cannot be identified by high-order repeats in carnations. We analyzed allele-specific expression in three tissues and the relationship between genome architecture and gene expression in the haplotypes. This demonstrated that the length of the genes, coding sequences, and introns, the exon numbers and the transposable element insertions correlate with gene expression ratios and levels. The insertions of transposable elements repress expression in gene regulatory networks in carnation. This gap-free finished T2T carnation genome provides a valuable resource to illustrate the genome characteristics and for functional genomics analysis in further studies and molecular breeding.

Multi-integrated genomic data for Passiflora foetida provides insights into genome size evolution and floral development in Passiflora

Passiflora is a plant genus known for its extremely distinctive and colorful flowers and a wide range of genome size variation. However, how genome characteristics are related to flower traits among Passiflora species remains poorly understood. Here, we assembled a chromosome-scale genome of P. foetida, which belongs to the same subgenus as the commercial passionfruit P. edulis. The genome of P. foetida is smaller (424.16 Mb) and contains fewer copies of long terminal repeat retrotransposons (LTR-RTs). The disparity in LTR-RTs is one of the main contributors to the differences in genome sizes between these two species and possibly in floral traits. Additionally, we observed variation in insertion times and copy numbers of LTR-RTs across different transposable element (TE) lineages. Then, by integrating transcriptomic data from 33 samples (eight floral organs and flower buds at three developmental stages) with phylogenomic and metabolomic data, we conducted an in-depth analysis of the expression, phylogeny, and copy number of MIKC-type MADS-box genes and identified essential biosynthetic genes responsible for flower color and scent from glandular bracts and other floral organs. Our study pinpoints LRT-RTs as an important player in genome size variation in Passiflora species and provides insights into future genetic improvement.

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

Functional Characterization of AP2/ERF Transcription Factors during Flower Development and Anthocyanin Biosynthesis Related Candidate Genes in Lycoris

The APETALA2/ethylene-responsive transcription factor (AP2/ERF) family has been extensively investigated because of its significant involvement in plant development, growth, fruit ripening, metabolism, and plant stress responses. To date, there has been little investigation into how the AP2/ERF genes influence flower formation and anthocyanin biosynthesis in Lycoris. Herein, 80 putative LrAP2/ERF transcription factors (TFs) with complete open reading frames (ORFs) were retrieved from the Lycoris transcriptome sequence data, which could be divided into five subfamilies dependent on their complete protein sequences. Furthermore, our findings demonstrated that genes belonging to the same subfamily had structural similarities and conserved motifs. LrAP2/ERF genes were analyzed for playing an important role in plant growth, water deprivation, and flower formation by means of gene ontology (GO) enrichment analysis. The expression pattern of the LrAP2/ERF genes differed across tissues and might be important for Lycoris growth and flower development. In response to methyl jasmonate (MeJA) exposure and drought stress, the expression of each LrAP2/ERF gene varied across tissues and time. Moreover, a total of 20 anthocyanin components were characterized using ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis, and pelargonidin-3-O-glucoside-5-O-arabinoside was identified as the major anthocyanin aglycone responsible for the coloration of the red petals in Lycoris. In addition, we mapped the relationships between genes and metabolites and found that LrAP2/ERF16 is strongly linked to pelargonidin accumulation in Lycoris petals. These findings provide the basic conceptual groundwork for future research into the molecular underpinnings and regulation mechanisms of AP2/ERF TFs in anthocyanin accumulation and Lycoris floral development.

A chromosome-scale and haplotype-resolved genome assembly of carnation (Dianthus caryophyllus) based on high-fidelity sequencing

Dianthus caryophyllus is an economic species often considered excellent cut flowers and is suitable for bouquets and gardens. Here, we assembled the haplotype-resolved genome of D. caryophyllus 'Aili' at the chromosome level for the first time. The total lengths of the two assembled haplotypes of carnation were 584.88 Mb for haplotype genome 1 (hap1) and 578.78 Mb for haplotype genome 2 (hap2), respectively. We predicted a total of 44,098 and 42,425 protein-coding genes, respectively. The remarkable structure variation was identified between two haplotypes. Moreover, we identified 403.80 Mb of transposable elements (TEs) in hap1, which accounted for 69.34% of the genome. In contrast, hap2 had 402.70 Mb of TEs, representing 69.61% of the genome. Long terminal repeats were the predominant transposable elements. Phylogenetic analysis showed that the species differentiation time between carnation and gypsophila was estimated to be ~54.43 MYA. The unique gene families of carnation genomes were identified in 'Aili' and previously published 'Francesco' and 'Scarlet Queen'. The assembled and annotated haplotype-resolved D. caryophyllus genome not only promises to facilitate molecular biology studies but also contributes to genome-level evolutionary studies.

A high-quality Bougainvillea genome provides new insights into evolutionary history and pigment biosynthetic pathways in the Caryophyllales

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana 'Mrs Butt', a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Genome-Wide Identification of the MYB and bHLH Families in Carnations and Expression Analysis at Different Floral Development Stages

Carnations are one of the most popular ornamental flowers in the world with varied flower colors that have long attracted breeders and consumers alike. The differences in carnation flower color are mainly the result of the accumulation of flavonoid compounds in the petals. Anthocyanins are a type of flavonoid compound that produce richer colors. The expression of anthocyanin biosynthetic genes is mainly regulated by MYB and bHLH transcription factors. However, these TFs have not been comprehensively reported in popular carnation cultivars. Herein, 106 MYB and 125 bHLH genes were identified in the carnation genome. Gene structure and protein motif analyses show that members of the same subgroup have similar exon/intron and motif organization. Phylogenetic analysis combining the MYB and bHLH TFs from Arabidopsis thaliana separates the carnation DcaMYBs and DcabHLHs into 20 subgroups each. Gene expression (RNAseq) and phylogenetic analysis shows that DcaMYB13 in subgroup S4 and DcabHLH125 in subgroup IIIf have similar expression patterns to those of DFR, ANS, and GT/AT, which regulate anthocyanin accumulation, in the coloring of carnations, and in red-flowered and white-flowered carnations, DcaMYB13 and DcabHLH125 are likely the key genes responsible for the formation of red petals in carnations. These results lay a foundation for the study of MYB and bHLH TFs in carnations and provide valuable information for the functional verification of these genes in studies of tissue-specific regulation of anthocyanin biosynthesis.

Widely targeted metabolomics reveals the antioxidant and anticancer activities of different colors of Dianthus caryophyllus

Carnation is edible flower that has potent antioxidant properties and is used in traditional Chinese medicinal system and food industry. The phytochemicals responsible for these various proprieties, however, are not fully understood. Thus, in order to recognize metabolite diversity and variability in carnation flowers of different colors and to discover key metabolites that contribute to the differences in antioxidant and anticancer activities, widely targeted LC-MS/MS-based metabolomics analysis was conducted on purple, green, yellow, and white carnation flowers. We identified and chemically categorized 932 metabolites. Metabolic compounds varied significantly with flower color. Several flavonoids, organic acids, phenolic acids, and nucleotides and their derivatives were found to be specific differential metabolites in purple flowers. A total of 128 key differential metabolites were screened. The purple flowers were found to have the highest antioxidant and anticancer activities compared to the other colored flowers. Correlation analysis revealed that the 6-hydroxykaempferol-3,6-O-diglucoside, 6-hydroxykaempferol-7-O-glucoside, quercetin-3-O-sophoroside, and 2'-deoxyguanosine were found to be the major constituents of the antioxidant and anticancer activities. 2'-Deoxyguanosine has effective antiproliferative activity against A549 and U2OS cells for the first report. At the same time, the combination of 2'-deoxyguanosine with 6-hydroxykaempferol-3, 6-O-diglucoside, or quercetin-3-O-sophoroside have also been found to increase the antitumor activity of 2'-deoxyguanosine. These discoveries enrich information on the phytochemical composition of carnation of different colors and provide resources for the overall use and improvement of carnation flowers quality.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

Charting plant gene functions in the multi-omics and single-cell era

Despite the increased access to high-quality plant genome sequences, the set of genes with a known function remains far from complete. With the advent of novel bulk and single-cell omics profiling methods, we are entering a new era where advanced and highly integrative functional annotation strategies are being developed to elucidate the functions of all plant genes. Here, we review different multi-omics approaches to improve functional and regulatory gene characterization and highlight the power of machine learning and network biology to fully exploit the complementary information embedded in different omics layers. Finally, we discuss the potential of emerging single-cell methods and algorithms to further increase the resolution, allowing generation of functional insights about plant biology.

Multiomics analysis reveals the mechanisms underlying the different floral colors and fragrances of Rosa hybrida cultivars

The color and fragrance of rose flowers affect their commercial value. However, several rose varieties with new floral colors developed by the bud mutation method lost their fragrance during the breeding process, raising the question: Is there a relationship between floral color and aroma traits? Rose cultivar 'Yellow Island' (YI) with intensely aroma and yellow petals, while its bud mutant 'Past Feeling' (PF) with light aroma and pink petals mixing some yellow, two cultivars were used to explore this question using multiomics approaches. We investigated the genomic polymorphisms between PF and YI by whole-genome resequencing. 71 differentially abundant metabolites and 155 related differentially expressed genes identified in petals between PF and YI. From this, we constructed a model of metabolic changes affecting floral color and fragrance integrating shikimate, terpenoid, carotenoid, and green leaf volatile metabolites and predicted the associated key genes and transcription factors. This study provides a reference for understanding the molecular mechanism of variation in rose floral color and aroma traits.

Characterization of the WRKY Gene Family Related to Anthocyanin Biosynthesis and the Regulation Mechanism under Drought Stress and Methyl Jasmonate Treatment in Lycoris radiata

Lycoris radiata, belonging to the Amaryllidaceae family, is a well-known Chinese traditional medicinal plant and susceptible to many stresses. WRKY proteins are one of the largest families of transcription factors (TFs) in plants and play significant functions in regulating physiological metabolisms and abiotic stress responses. The WRKY TF family has been identified and investigated in many medicinal plants, but its members and functions are not identified in L. radiata. In this study, a total of 31 L. radiata WRKY (LrWRKY) genes were identified based on the transcriptome-sequencing data. Next, the LrWRKYs were divided into three major clades (Group I-III) based on the WRKY domains. A motif analysis showed the members within same group shared a similar motif component, indicating a conservational function. Furthermore, subcellular localization analysis exhibited that most LrWRKYs were localized in the nucleus. The expression pattern of the LrWRKY genes differed across tissues and might be important for Lycoris growth and flower development. There were large differences among the LrWRKYs based on the transcriptional levels under drought stress and MeJA treatments. Moreover, a total of 18 anthocyanin components were characterized using an ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis and pelargonidin-3-O-glucoside-5-O-arabinoside as well as cyanidin-3-O-sambubioside were identified as the major anthocyanin aglycones responsible for the coloration of the red petals in L. radiata. We further established a gene-to-metabolite correlation network and identified LrWRKY3 and LrWRKY27 significant association with the accumulation of pelargonidin-3-O-glucoside-5-O-arabinoside in the Lycoris red petals. These results provide an important theoretical basis for further exploring the molecular basis and regulatory mechanism of WRKY TFs in anthocyanin biosynthesis and in response to drought stress and MeJA treatment.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Comprehensive physiological, transcriptomic, and metabolomic analyses reveal the synergistic mechanism of Bacillus pumilus G5 combined with silicon alleviate oxidative stress in drought-stressed Glycyrrhiza uralensis Fisch

Glycyrrhiza uralensis Fisch. is often cultivated in arid, semi-arid, and salt-affected regions that suffer from drought stress, which leads to the accumulation of reactive oxygen species (ROS), thus causing oxidative stress. Plant growth-promoting bacteria (PGPB) and silicon (Si) have been widely reported to be beneficial in improving the tolerance of plants to drought stress by maintaining plant ROS homeostasis. Herein, combining physiological, transcriptomic, and metabolomic analyses, we investigated the response of the antioxidant system of G. uralensis seedlings under drought stress to Bacillus pumilus (G5) and/or Si treatment. The results showed that drought stress caused the overproduction of ROS, accompanied by the low efficiency of antioxidants [i.e., superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), the ascorbate (AsA)-glutathione (GSH) pool, total carotenoids, and total flavonoids]. Inversely, supplementation with G5 and/or Si enhanced the antioxidant defense system in drought-stressed G. uralensis seedlings, and the complex regulation of the combination of G5 and Si differed from that of G5 or Si alone. The combination of G5 and Si enhanced the antioxidant enzyme system, accelerated the AsA-GSH cycle, and triggered the carotenoid and flavonoid metabolism, which acted in combination via different pathways to eliminate the excess ROS induced by drought stress, thereby alleviating oxidative stress. These findings provide new insights into the comparative and synergistic roles of PGPB and Si in the antioxidant system of plants exposed to drought and a guide for the application of PGPB combined with Si to modulate the tolerance of plants to stress.