OrchidBase 5.0: updates of the orchid genome knowledgebase
BMC Plant Biology volume 22, Article number: 557 (2022)
Containing the largest number of species, the orchid family provides not only materials for studying plant evolution and environmental adaptation, but economically and culturally important ornamental plants for human society. Previously, we collected genome and transcriptome information of Dendrobium catenatum, Phalaenopsis equestris, and Apostasia shenzhenica which belong to two different subfamilies of Orchidaceae, and developed user-friendly tools to explore the orchid genetic sequences in the OrchidBase 4.0. The OrchidBase 4.0 offers the opportunity for plant science community to compare orchid genomes and transcriptomes and retrieve orchid sequences for further study.
In the year 2022, two whole-genome sequences of Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis, were de novo sequenced, assembled and analyzed. In addition, systemic transcriptomes from these two species were also established. Therefore, we included these datasets to develop the new version of OrchidBase 5.0. In addition, three new functions including synteny, gene order, and miRNA information were also developed for orchid genome comparisons and miRNA characterization.
OrchidBase 5.0 extended the genetic information to three orchid subfamilies (including five orchid species) and provided new tools for orchid researchers to analyze orchid genomes and transcriptomes. The online resources can be accessed at https://cosbi.ee.ncku.edu.tw/orchidbase5/
Orchids have been considered mysterious and fascinating to scientists, collectors, growers, and general public. Orchidaceae represents approximately 10 percent of angiosperms and is one of the largest angiosperm family containing more than 25,000 species categorized in about 900 genera . The basic blueprint of orchid flower including three sepals at outermost floral whorl, two lateral petals and one elaborated medial labellum at second whorl, and gynostemium (fused male and female reproductive organ at central whorl) with diversified shapes and color patterns has captured the vision of people worldwide for hundreds of years . The floral features of orchid, especially the delicate and complex structures of labellum and gynostemium, attract and interact with pollinators to achieve reproductive assurance. In addition to floral diversification, orchids also evolved distinct lifestyles, such as terrestrial growth, epiphytism, and mycoheterotrophy, for adaptation and radiation of heterogeneously ecological niches. It is well established that all orchid species could be monophylletically gathered into Orchidaceae composed of five subfamilies including Apostasioideae (most basal), Epidendroideae, Vanilloideae, Cypripedioideae, and Orchidoideae.
The previous version of OrchidBase 4.0 accommodates the Sanger-sequenced transcriptomes from 11 various tissues of three Phalaenopsis species, Illumina-sequenced floral transcriptomes of 10 orchid species across five Orchidaceae subfamilies, as well as whole-genome sequence of Phalaenopsis equestris (Epidendroideae) , Dendrobium catenatum (Epidendroideae) , and Apostasia shenzhenica (Apostasioideae)  and their relative transcriptomes [6,7,8,9]. However, the complete genomic sequences from Vanilloideae, Cypripedioideae, and Orchidoideae are still absent in the OrchidBase, hindering the genetic comparisons among orchid species distributed in more than two different orchid subfamilies.
Taxonomic diversity of orchids was obviously observed in two recently expanded sister-subfamilies: Orchidoideae (occupying about 14% species in Orchidaceae) and especially Epidendroideae (about 84%). Floral morphological differences between Orchidoideae and Epidendroideae includes loss of stamen vasculature in Orchidoideae and development of hard pollinia in Epidendroideae . In addition, species in Orchidoideae have a shared terrestrial habitat while those in Epidendroideae are epiphytic growth.
Recent whole-genome sequenced orchids just concentrated on the species in Epidendroideae, Vanilloideae, and Apostasioideae [3,4,5, 11,12,13,14,15,16]. Thus, increasing the genome information of species in Orchidoideae as well as in Cypripedioideae is very important for understanding the orchid genome evolution and the gene function. In the OrchidBase 5.0, we added two whole-genome sequenced Orchidoideae species, Platanthera zijinensis and Platanthera guangdongensis (Fig. 1), and their transcriptomes derived from various tissues. Furthermore, three tools were newly developed for synteny, gene order and miRNA comparison among five orchid genomes. The data and tools launched in the OrchidBased 5.0 make it an excellent resource for orchid biology research.
Construction and content
Implementation and architecture
Figure 3 shows the main page of each orchid species in OrchidBase 5.0 which includes genome (described in OrchidBase 4.0) and transcriptome (described in OrchidBase and OrchidBase 2.0) information. The OrchidBase 5.0 simplifies the workflow for huge and complicated biological data analysis and visualization. OrchidBase 5.0 is an open-access, web-available portal that integrates the available data for the genomes of the five orchid species (P. equestris, D. catenatum, A. shenzhenica, Pl. zijinensis, and Pl. guangdongensis) and the related transcriptome information.
Expanded database content
Pl. zijinensis and Pl. guangdongensis both have a karyotype of 2 N = 2X = 42 chromosomes. A total of 441.16 Gb and 414.20 Gb data were generated respectively for Pl. zijinensis and Pl. guangdongensis using PacBio technologies . The total length of the genome assembly was 4.19 Gb with a contig N50 value of 1.77 Mb for the Pl. zijinensis genome and 4.20 Gb with a contig N50 value of 1.57 Mb for the Pl. guangdongensis genome (Table 2) . Hi-C libraries were built and sequenced by Illumina technology to reconstruct physical maps by ordering and clustering the assembled scaffolds into 21 pseudomolecules for each species. The raw data and whole genome-assembled scaffold sequences of the two Platanthera species (BioProject PRJNA739531) were downloaded from the National Center for Biotechnology Information (NCBI) database. Statistics of the added orchid genomes is presented in Table 2. Based on these datasets, protein-coding genes were predicted by combination of homology-based prediction, de novo gene prediction, RNA sequence-aided prediction  for each species. Twenty-four thousands five hundreds and thirteen and 22,559 protein-coding genes in the respective genome of Pl. zijinensis and Pl. guangdongensis were predicted. Furthermore, 31 and 33 miRNAs were identified in the Pl. zijinensis and Pl. guangdongensis genomes, respectively . Each predicted gene is assigned to a specific Gene ID as the identifier of the predicted gene. The specific genes could be selected to investigate their annotated functions involved in biological processes of interest.
The transcriptomics data derived from the two Platanthera species were also downloaded from BioProject PRJNA739531. All RNA sequencing (RNA-seq) reads were mapped to the predicted genes and the fragments per kilobase of transcript per million mapped reads (FPKM) values for each gene in various tissues and different developmental stages were calculated. All the gene expression values have been integrated into the updated version of OrchidBase 5.0.
Utility and discussion
Searching the genome information of the two Orchidoideae species in the database
Through the web interface, the genome information of the two Platanthera species contained in OrchidBase 5.0 could be freely obtained. The information can be linked via the “Orchid Genome" icon (Fig. 3, step 1). With the interface, a page allows users to select one of the three already existed orchid genomes (P. equestris, D. catenatum, A. shenzhenica) or two newly added orchid genomes (Pl. zijinensis, and Pl. guangdongensis) (Fig. 3, step 2). Users then could access the Genome browser (Fig. 3, step 3), Gene annotation (Fig. 3, step 4), Gene expression (Fig. 3, step 5), Metabolism pathway (Fig. 3, step 6), BLAST (Fig. 3, step 7), Synteny (Fig. 3, step 8), Gene order (Fig. 3, step 9), and miRNA (Fig. 3, step 10) for querying the genome and obtain the gene information and comparison analysis in the selected orchid genome. Owing to Genome browser, Gene annotation, Gene expression, Metabolism pathway, and BLAST have been introduced in OrchidBase 4.0 , here we only explained the new pages of Synteny, Gene order and miRNA in details.
Database user protocol
Synteny is a locally conserved gene order found in the compared genomes. Two species that have recently diverged from a common ancestor might be expected to share a similar set of genes positioned along the DNA strand in the same order. Genomic synteny comparisons among species have provided an additional opportunity to study evolutionary trajectories that lead to diversity of chromosomal structure in many lineages . The “Synteny” page is a graphical interface for displaying syntenic relationships (predicted using MCScanX ) between the selected genomic region and potential orthologous region of the other genome (Fig. 4). For visiting the “Synteny” page, users could click the Orchid Genome (Fig. 3, step 1), and then choose one of the orchid species (Fig. 3, step 2). Users will be directed to the main page of the selected genome and then enter the “Synteny page” (Fig. 3, step 8). In the “Synteny page”, users could select again one of the orchid species (Fig. 4, step 1). Users then could anchor one of the scaffolds (Fig. 4, step 2) in the selected orchid species and tick the other orchid genome for further comparison (Fig. 4, step 3). When users click the “search” icon at the left corner of the page (Fig. 4, step 4), they will be lead to the next page for showing the result indicating the orthologous blocks between the selected genomic scaffolds. In the page, users could pick one of the syntenic blocks (Fig. 4, step 5) to further understand the relative position of the genes in the syntenic block (Fig. 4, step 6) and get the detailed gene information (Fig. 4, step 7 and 8) and also gene annotation linked to the Gene Annotation page (Fig. 4, step 9).
“Gene Order” page
We provided another way to display permutation of gene arrangement based on the BLASTP and genomic location of genes, and the system we named here is “Gene Order” (Fig. 5). Users would not necessary to specify a scaffold, but select a unique “gene ID” of the genome. Based on the best hit of the BLASTP, the system provides the location of the orthologous genes in the other two selected genomes, and further indicates the position of the other homologous genes among the three compared species. To help users execute “Gene Order” analysis, an intuitive graphic interface was developed (Fig. 5). By accessing the “Gene Order” page, orchid genome could be selected (Fig. 3, step 1), and one of the orchid species could be opted (Fig. 3, step 2). In the main page of the selected orchid species, users could click the Gene order (Fig. 3, step 9). In the “Gene Order” page, users have to select one of the orchid species (Fig. 5, step 1) and offer the Gene ID (Fig. 5, step 2) for showing the best BLASTP hits among the other orchid genomes. The Gene ID marked by blue color could be clicked (Fig. 5, step 3), then users would be directed to the next page presenting the detailed comparison results among orchid orthologous proteins (Fig. 5, step 4). In this page, users could further select two orchid genomes (Fig. 5, step 5 and 6). After clicking the “Enter” icon (Fig. 5, step 7), the comparison of gene order among the three chosen genomes will be displayed (Fig. 5, step 8). In this page, the red rectangle at the middle row is the specified gene in the orchid genome selected in (Fig. 5, step 1), the red one at the top row is the orthologous gene in the genome selected in (Fig. 5, step 5), and the red one at the bottom row is another orthologous hit in the genome selected in (Fig. 5, step 6). User could point each of the colored rectangle to have the annotated information of the gene (Fig. 5, step 9 and 10) or each line connected between two orthologs to get the detailed comparison information (Fig. 5, step 11 and 12). The page also provides “Zoom in” and “Zoom out” function (Fig. 5, step 13) to enlarge or narrow specific region in the scaffold of the selected genome and offers a movable gray region to navigate the selected zone (Fig. 5, step 14).
miRNAs are small noncoding RNAs involved in regulating many biological processes at the post-transcriptional and translational levels. The molecular function of miRNAs is suppression the expression of their targets. For providing pre-miRNAs of orchid genomes, we identified them by using INFERNAL  to search the Rfam database. We further used miRLocator  to predict mature miRNAs. For identification of orchid miRNA targets, we uploaded the mature miRNAs sequences and orchid coding sequences (CDSs) to TAPIR , and adopted RNAhybrid search engine equipped in the server to get the target sequences. For going to the miRNA page, users could visit the Orchid Genome page (Fig. 3, step 1), then select one of the orchid species (Fig. 3, step 2). Users will step into the main page of the selected orchid species and then could click the miRNA page (Fig. 3, step 10). In the miRNA page, users have the further choice to pick one of the orchid species (Fig. 6, step 1) and click the View/Search item to browse the miRNAs in the selected orchid species (Fig. 6, step 2). Users also could type “miRNA ID” or keywords to query the miRNA information (Fig. 6, step 3). In this page, users could get more detailed information by accessing one of the “mature miRNAs” to obtain mature miRNA sequence, its annotation, and secondary structure of its pre-miRNA (Fig. 6. step 4). Users could also click one of the “number of targets” to visualize the “target mRNA” information (Fig. 6, step 5), and find the target gene names (Fig. 6, step 6). Users could click the target gene name (Fig. 6, step 7) and go to the “Gene Annotation” page (Fig. 6, step 8). Finally, the miRNA targeting site in the target mRNA could be shown by clicking one of choice at the “Link” (Fig. 6, step 9 and step 10).
A case study
NO APICAL MERISTEM (NAC), ARABIDOPSIS THALIANA ACTIVATING FACTOR (AtAF), and CUP-SHAPED COTYLEDON (CUC) encode distinct members of the plant NAC family . This family is suggested to be specific for land plants . Previous researches concluded that CUCs have conserved function in gynoecium and ovule development across angiosperms and some of CUC transcripts (for example Arabidopsis CUC1 and CUC2) were regulated by miR-164 . Because the orchid gynostemium (fused by androecium and genoecium floral organs) is distinct from those in most angiosperm species, we are interested in the function of CUC homologs played in orchid gynostemium development. Arabidopsis CUC1 (NCBI accession number: BAB20598) was used as query to BLASTP against Apostasia protein database, and the Apostasia Ash001107 located at fragScaff_scaffold_67: 6,149,945–6,151,832 was hit. We then went to “Synteny” page (Fig. 7) and selected library as “Apostasia”, Scaffold/Chromosome as “fragScaff_scaffold_67”, and chose any one of the orchid genomes (in this case: Phalaenopsis), then clicked search (Fig. 7A). The page was changed to the synteny blocks screened between Apostasia and Phalaenopsis genomes (Fig. 7B). We chose one of the lines connected Apostasia and Phalaenopsis blocks and searched the gene location to match Phalaenopsis scaffold harboring CUC gene (Fig. 7C). We selected the synteny block containing Apostasia Ash001107 to show the detailed information of this region embedding order and direction of genes and orthologous relationship among each gene (Fig. 7D). We found that the CUC gene in Phalaenopsis is Peq010053 and the CUC genes in these two orchid species formed a synteny block. Finally, the Ash001107 and Peq010053 could be clicked to reveal their gene annotation information (Fig. 7E).
“Gene Order” also could be applied to reveal the permutation of gene arrangement among orchid genomes. We filled in the gene ID (in this case: Ash001107) into the “Search” blank at the “Gene Order” page and the putative orthologs in other orchid genomes were shown (Fig. 8A). The Ash001107 marked by blue color was clicked. We then selected two orchid genomes for comparison (in this case: Phalaenopsis and Pl. guangdongensis, Fig. 8B), and clicked “Enter”. We clicked “ ± ” to zoom in/out the region around the indicated gene ID (marked by red color) to adjust the comparison region among the selected orchid genomes (Fig. 8B). We found that the regions around the CUC genes in Apostasia and Phalaenopsis may have synteny relationship since many genes in the compared regions are orthologs. On the contrary, the regions around the CUC genes in Apostasia and Pl. guangdongensis may not have synteny relationship since many neighboring genes of the CUC gene in Apostasia could not find orthologs nearby the CUC gene in Pl. guangdongensis. We also could click the red or blue rectangle to reveal the specific gene annotation (Fig. 8C). One of the advantages of “Gene order” tool compared to the “Synteny” is that if the putative orthologs do not locate in the corresponding synteny block, the putative orthologs and its neighboring genes could still be displayed and compared. The other advantage is that the “Gene order” tool could simultaneously compare three orchid genomes.
Arabidopsis CUC genes were known to be regulated by miR-164 . To check whether Ash001107 (CUC gene in Apostasia) is also regulated by miR-164, we entered “miRNA” page, and filled Ash001107 in the “Target mRNA”, then clicked “View/Search”. The results showed that Ash001107 is regulated by 24 different miRNAs, one of which is Ash_miRNA_000001_5p (Fig. 9A). Note that Ash_miRNA_000001_5p is the miR-164 in Apostasia and it has eight targets (Fig. 9B and C). Therefore, we concluded that CUC gene in Apostasia is also regulated by miR-164. The detailed information of the miR-164 binding sites in Ash001107 mRNA could be seen in Fig. 9D. The miRNA binding site information is useful for biologists to design experiments for exploring the regulatory mechanism of CUC-like gene controlled by miR-164 in orchid.
We added two whole genome sequences of Orchidoideae species (Pl. zijinensis and Pl. guangdongensis) and their transcriptomes into OrchidBase 5.0. In addition, two functions for genome comparison (synteny and gene order) and miRNA characterization were developed in this version. Both expanding the Orchidoideae genomes and creating new tools for exploring knowledge embedded in the nucleotide sequences provide the opportunity for the users to get novel insights into the conservation and diversification of orchid genomes. Users could precisely and efficiently adopt “synteny” and “gene order” tools to find specific collinear region among the five orchid genomes (P. equestris, D. catenatum, A. shenzhenica, Pl. zijinensis and Pl. guangdongensis). Although studies of miRNA function involved in orchid growth and development are still scarce [27, 28], we believe that the researches will be increased after the miRNA characterization tool in the OrchidBase 5.0. has been released. In the future, we will continue to include new orchid whole genome sequences in the OrchidBase and also develop new -omic analysis tools for plant scientists.
Availability of data and materials
The raw data and whole genome-assembled scaffold sequences for the two Platanthera species (BioProject PRJNA739531) were downloaded from the National Center for Biotechnology Information (NCBI) database. The transcriptome data derived from the two Platanthera species were also downloaded from BioProject PRJNA739531 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA739531).
- A. shenzhenica :
- D. catenatum :
Kyoto Encyclopedia of Genes and Genomes
Fragments Per Kilobase of transcript per Millions mapped reads
- P. equestris :
- Pl. zijinensis :
- Pl. guangdongensis :
Structured Query Language
Chase MW, Cameron KM, Freudenstein JV, Pridgeon AM, Salazar G, Van Den Berg C, et al. An updated classification of Orchidaceae. Bot J Linn Soc. 2015;177(2):151–74.
Albert VA, Carretero-Paulet L. A genome to unveil the mysteries of orchids. Nat Genet. 2015;47(1):3–4.
Cai J, Liu X, Vanneste K, Proost S, Tsai WC, Liu KW, et al. The genome sequence of the orchid Phalaenopsis equestris. Nat Genet. 2015;47(1):65–72.
Zhang GQ, Xu Q, Bian C, Tsai WC, Yeh CM, Liu KW, et al. The Dendrobium catenatum Lindl. genome sequence provides insights into polysaccharide synthase, flower development and adaptive evolution. Sci Rep. 2016;6:19029.
Zhang GQ, Liu KW, Li Z, Lohaus R, Hsiao YY, Niu SC, et al. The Apostasia genome and the evolution of orchids. Nature. 2017;579(7672):379–83.
Fu CH, Chen YW, Hsiao YY, Pan ZJ, Liu ZJ, Huang YM, et al. OrchidBase: a collection of sequences of the transcriptome derived from orchids. Plant Cell Physiol. 2011;52(2):238–43.
Tsai WC, Fu CH, Hsiao YY, Huang YM, Chen LJ, Liu ZJ, et al. OrchidBase 2.0: comprehensive collection of Orchidaceae floral transcriptomes. Plant Cell Physiol. 2013;54(2):e7.
Tsai WC, Fu CH, Hsiao YY, Wu WL, Zhang DY, Lan SR, et al. OrchidBase 3.0: a resource for studying gene function and genome evolution in orchids. J Fujian Agric Univ. 2019;48:440–6.
Hsiao YY, Fu CH, Ho SY, Li CI, Chen YY, Wu WL, et al. OrchidBase 4.0: a database for orchid genomics and molecular biology. BMC Plant Biol. 2021;21:371.
Rudall PJ, Bateman RM. Roles of synorganisation, zygomorphy and heterotopy in floral evolution: the gynostemium and labellum of orchids and other lilioid monocots. Biol Rev Camb Philos Soc. 2002;77(3):403–41.
Yan L, Wang X, Liu H, Tian Y, Lian J, Yang R, et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese orchid herb. Mol Plant. 2015;8(6):922–34.
Chao YT, Chen WC, Chen CY, Ho HY, Yeh CH, Kuo YT, et al. Chromosome-level assembly, genetic and physical mapping of Phalaenopsis aphrodite genome provides new insights into species adaptation and resources for orchid breeding. Plant Biotechnol J. 2018;16(12):2027–41.
Yuan Y, Jin X, Liu J, Zhao X, Zhou J, Wang X, et al. The Gastrodia elata genome provides insights into plant adaptation to heterotrophy. Nature Comm. 2018;9(1):1615.
Hasing T, Tang H, Brym M, Khazi F, Huang T, Chambers AH. A phased Vanilla planifolia genome enables genetic improvement of flavour and production. Nature Food. 2020;1:811–9.
Yang FX, Gao J, Wei YL, Ren R, Zhang GQ, Lu CQ, et al. The genome of Cymbidium sinense revealed the evolution of orchid traits. Plant Biotechnol J. 2021;19:2501–16.
Ai Y, Li Z, Sun WH, Chen J, Zhang D, Ma L, et al. The Cymbidium genome reveals the evolution of unique morphological traits. Hortic Res. 2021;8(1):255.
Li MH, Liu KW, Li Z, Lu HC, Ye QL, Zhang D, et al. Genomes of leafy and leafless Platanthera orchids provide insights into the evolution of mycoheterotrophy. Nature Plants. 2022;8(4):373–88.
Keilwagen J, Hartung F, Paulini M, Twardziok SO, Grau J. Combining RNA-seq data and homology-based gene prediction for plants, animals and fungi. BMC Bioinformatics. 2018;19(1):189.
Lyons E, Pedersen B, Kane J, Alam M, Ming R, Tang H, et al. Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with rosids. Plant Physiol. 2008;148(4):1772–81.
Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012;40(7):e49.
Nawrocki EP, Kolbe DL, Eddy SR. Infernal 1.0: inference of RNA alignments. Bioinformatics. 2009;25(10):1335–7.
Cui H, Zhai J, Ma C. miRLocator: machine learning-based prediction of mature microRNAs within Plant pre-miRNA sequences. PLoS ONE. 2015;10(11):e0142753.
Bonnet E, He Y, Billiau K, Van de Peer Y. TAPIR, a web server for the prediction of plant microRNA targets, including target mimics. Bioinformatics. 2010;26(12):1566–8.
Ooka H, Satoh K, Doi K, Nagata T, Otomo Y, Murakami K, et al. Comprehensive analysis of NAC family genes in Oryza sativa and Arabidopsis thaliana. DNA Res. 2003;10(6):239–47.
Finet C, Timme RE, Delwiche CF, Marletaz F. Multigene phylogeny of the green lineage reveals the origin and diversification of land plants. Curr Biol. 2010;20(24):2217–22.
Vialette-Guiraud ACM, Adam H, Finet C, Jasinski S, Jouannic S, Scutt CP. Insights from ANA-grade angiosperms into the early evolution of cup-shaped cotyledon genes. Ann Bot. 2011;107(9):1511–9.
Aceto S, Sica M, De Paolo S, D’Argenio V, Cantiello P, Salvatore F, et al. The analysis of the inflorescence miRNome of the orchid Orchis italica reveals a DEF-like MADS-box gene as a new miRNA target. PLoS ONE. 2014;9(5):e97839.
Xu Z, Liu Q, Chen Y, He Y, Hu F. miR390 family of Cymbidium goeringii is involved in the development of reproductive organs in transgenic Arabidopsis. BMC Plant Biol. 2022;22(1):149.
We thank National Center for High-performance Computing (NCHC, https://www.nchc.org.tw) to kindly provide the cloud space for storage and assessment of transcriptomic data.
Publication costs was funded by the Ministry of Science and Technology, Taiwan (MOST 110–2313-B-006–002-MY3, MOST 110–2221-E-006–198-MY3 and MOST-111–2221-E-006–151-MY3), National Key R&D Program of China (Grant No.2019YFD1000400), and National Natural Science Foundation of China [no. 31870199]. The funders have no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, YY., Li, C., Hsiao, YY. et al. OrchidBase 5.0: updates of the orchid genome knowledgebase. BMC Plant Biol 22, 557 (2022). https://doi.org/10.1186/s12870-022-03955-5