ORIGINAL RESEARCH article

Front. Plant Sci., 12 May 2025

Sec. Plant Breeding

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1586225

The mechanism of tetraploidization in tree peony, and its implications for speciation and evolution of genus Paeonia L.

  • 1State Key Laboratory of Efficient Production of Forest Resources, National Engineering Research Center for Floriculture, Peony International Institute, School of Landscape Architecture, Beijing Forestry University, Beijing, China
  • 2Beijing Key Laboratory of Ornamental Plants Germplasm Innovation & Molecular Breeding, National Engineering Research Center for Floriculture, Peony International Institute, School of Landscape Architecture, Beijing Forestry University, Beijing, China

Introduction: Polyploidization is not only an important driving force for plant speciation and evolution, but also an effective approach for plant domestication and improvement. Polyploid taxa are quite common in Paeonia section Paeonia (herbaceous peonies), but very rare in section Moutan (tree peonies), which are known as the ‘king of flowers’.

Methods: In this paper, we studied the origination of a tetraploid tree peony, P. × lemoinei ‘Golden Era’ (‘GE’, AABB), by comparing its karyotype with its seed parent, P. × lemoinei ‘Golden Isles’ (‘GI’, AB), based on genomic in situ hybridization (GISH) and rDNA fluorescence in situ hybridization (FISH). The meiosis behaviors of ‘GI’ were observed to make clear the cytogenetic process of unreduced (2n) gamete generation.

Results: Four chromosomes with inter-genome translocations were identified in ‘GE’, two of which might be reciprocal translocations. The 2n female gametes via first division restitution (FDR) from ‘GI’ might play an important role in tetraploidization of ‘GE’.

Disucssion: The distant hybridization between intersterile species from different subsections of section Moutan probably promotes the tetraploidization of tree peonies by facilitating 2n gamete production. The mechanism of tetraploidization in section Moutan is highly consistent with that in section Paeonia, but is inseparable with the assist from mankind. The divergence of life history between tree peonies and herbaceous peonies is speculated to contribute to the different level of polyploidization, and distinct tempo of speciation and evolution, between section Moutan and section Paeonia in genus Paeonia. These findings bring new insights to polyploid breeding, speciation and evolution in genus Paeonia.

1 Introduction

Polyploidization is one of the most important driving forces for plant speciation and evolution (Alix et al., 2017; De Storme and Mason, 2014; Jiao et al., 2011; Rieseberg and Willis, 2007; Soltis et al., 2009). The effects of polyploidy such as enlarged organs, vigorous growth and enhanced stress resistance are beneficial for plants to adapt to harsher environments in nature, and are also used to improve yield, quality, stress resistance and other important traits in plant breeding. So far, polyploid breeding has been widely carried out in the cultivation of new plant varieties in cereals, vegetables, fruits, trees, ornamental plants, and industrial crops (Chen, 2010; Fang and Morrell, 2016; Kang and Wei, 2022; Salman-Minkov et al., 2016; Zhang et al., 2019).

Tree peonies have been well known ornamental plants in China for more than 1000 years, and are popular in many other countries now (Cheng, 2007; Page, 2005). There are about 34 species in genus Paeonia (Paeoniaceae), which are classified into three sections (Cheng, 2007; Hong, 2010) or two subgenera (Hong, 2021). Most species of genus Paeonia are diploids (2n=2x=10), therein, only section Paeonia (belonging to subgenus Paeonia in the new taxonomic system of Hong (2021)) has natural tetraploid species, subspecies or populations (Hong, 2021). While all species of section Moutan (subgenus Moutan of Hong (2021)), known as tree peonies, are diploids, in addition, only two triploids (Zhong et al., 2024a; Li and Zhang, 1982) and several tetraploids (Zhong et al., 2024b; Zhong et al., 2023) have been found in the past decades in cultivated varieties. The tetraploid taxa of section Paeonia are mostly allotetraploids, which mainly originated from natural hybridization and polyploidization among diploid species. Moreover, there are some naturally formed autotetraploid populations in some species of section Paeonia (Ferguson and Sang, 2001; Sang et al., 2004; Zhou et al., 2021). Different from section Paeonia, the plants in section Moutan do not form polyploids under natural conditions (Hong, 2021), but only under artificial cultivation. There are great differences in the occurrence frequency and generation mode of polyploids between the two sections, implying that their polyploidization mechanisms might be different. Therefore, revealing the origination mechanism of polyploid tree peonies has important theoretical and practical significance for understanding the polyploidization mechanism of genus Paeonia, and promoting polyploid breeding of tree peonies.

The vast majority of the polyploid tree peonies are derived from P. × lemoinei, the artificially bred inter-subsectional hybrid tree peonies (Zhong et al., 2024a, 2024b, 2023). The inter-subsectional distant hybridization began around 1900s, with P. delavayi (B-genome) from subsection Delavayanae as the seed parents, and P. suffruticosa (A-genome) from subsection Vaginatae as the pollen parents (Hong, 2010; Page, 2005; Wister, 1995). Karyotype comparison based on GISH and rDNA FISH has shown that 2n gametes produced by the diploid seed parent, P. × lemoinei ‘High Noon’ (2n=2x=10, AB), play a key role in the formation of the allotriploid tree peony, P. × lemoinei ‘Oukan’ (2n=3x=15, ABC) (Zhong et al., 2024a), but the origination mechanism of allotetraploid tree peonies remains unclear. Although a few allotetraploid tree peonies have been found, the parents of them are mostly unknown. P. × lemoinei ‘Golden Era’ (‘GE’) is one of the earliest found allotetraploid tree peonies, which has higher fertility and better compatibility than many other cultivars in P. × lemoinei when crossing with herbaceous peonies, and is one of the most important pollen parents in intersectional distant hybridization (Page, 2005; Zhong et al., 2024b). The seed parent of ‘GE’ is P. × lemoinei ‘Golden Isles’ (‘GI’) (Page, 2005), a diploid intersubsectional hybrid tree peony. Although its pollen parent is unknown, ‘GE’ is still the best material available to study the origin of allotetraploid tree peony. The GISH study showed that the genome component of ‘GE’ is AABB (Zhong et al., 2024b), but its detailed karyotype information is lacking. Meanwhile, the karyotype of ‘GI’ has not been reported.

Unreduced gametogenesis is a key step in sexual polyploidization (Ramsey and Schemske, 1998). Different types of meiosis behaviors may produce unreduced (2n) gametes with diverse genetic compositions, thus affect the fertility and traits of polyploids (De Storme and Geelen, 2013; De Storme and Mason, 2014; Kreiner et al., 2017; Younis et al., 2014). Hence, the characteristics of meiosis behavior of diploid parents can provide important information to reveal the mechanism of polyploidization. Some studies have shown that ‘High Noon’ (2n=2x=10, AB) produces FDR (first division restitution) type 2n gametes through abnormal spindle localization, and transmits 10 chromosomes to the progeny intactly, resulting in the formation of the allotriploid, ‘Oukan’ (2n=3x=15, ABC) (Zhong et al., 2024a, 2019). However, the producing pathway and genetic type of 2n gametes from ‘GI’ are unclear.

Therefore, in this study, we used ‘GE’, a tetraploid tree peony, and its diploid seed parent, ‘GI’, as materials to conduct a karyotype comparison study by GISH and FISH, and observed the meiosis behaviors of ‘GI’. Based on the above results, the origination pathway of ‘GE’ was revealed, and the similarities and differences of polyploidization mechanisms between section Moutan and section Paeonia were discussed. The findings provide a theoretical basis for understanding the polyploidization mechanism of species in genus Paeonia, and for promoting polyploid breeding in section Moutan.

2 Materials and methods

2.1 Plant materials

Germinating buds of the tetraploid tree peony cultivar, P. × lemoinei ‘Golden Era’ (‘GE’), and its diploid seed parent, P. × lemoinei ‘Golden Isles’ (‘GI’), were collected for miotic chromosome preparations and karyotype analysis by GISH and FISH. Anthers of ‘GI’ were used for meiotic chromosome preparations and observation on meiosis behaviors. Fresh leaves of P. delavayi, one parent of P. × lemoinei, were employed for genomic DNA extraction and probe labeling for GISH in karyotype analysis of ‘GE’ and ‘GI’, and meiosis observation on ‘GI’, to distinguish the chromosomes of B-genome (P. delavayi), from those of A-genome (P. suffruticosa), the pollen parent of P. × lemoinei. All the materials were collected from living plants cultivated in the Peony Study Base of Beijing Forestry University, in the Jiufeng Forestry Park, Beijing.

2.2 Slides preparation

Germinating buds of ‘GE’ and ‘GI’ were taken in spring to prepare the metaphase chromosomes for karyotype analysis, according to the protocols described by Zhong et al. (2024a). Slides with well spread metaphase chromosomes were selected for the following in situ hybridization.

Young anthers of ‘GI’ in yellow green color were collected and fixed with Carnot fixative (anhydrous ethanol: glacial acetic acid = 3:1) at 4°C for 24 h, rinsed several times with distilled water, and then transferred to 70% ethanol, stored at -20°C in a refrigerator. Some of the fixed anthers were squashed on slides after rinsing with distilled water, and stained with modified Carbol-Fuchsin solution (G1165, Solarbio). The slides were observed under an optical microscope (Leica DM500, Germany) and photographed for meiosis observation.

Some other fixed anthers were processed as described by Dang et al. (2015) with some modifications for meiosis observation based on GISH. Pollen mother cells (PMCs) were suspended in mixed enzyme solution (3% cellulose + 0.3% pectinase + 1% snailase, W/V). Bottles were capped and vertically immersed in a 37°C water bath for 2.5 h. The suspension was then transferred into polyethylene centrifuge tubes and centrifuged at 2000 × g for 3 min. The supernatant was removed and the precipitate was resuspended in 200 μL Carnot fixative. Centrifugation was performed again as in the previous step. The precipitate was resuspended in 100 μL Carnot fixative. Then, the suspension was dropped on greaseless slides and dried rapidly using an alcohol flame.

2.3 Probe labeling and in situ hybridization

Genomic DNA extracted from young leaves of P. delavayi by plant genomic DNA kit (DP305, Tiangen) was used as GISH probe. 45S rDNA and 5S rDNA from maize were used as FISH probes. All the probes were labeled by nick translation kits, Dig-Nick Translation Mix (No. 11745816910, Roche) or Biotin-Nick Translation Mix (No. 11745824910, Roche), at 15°C for 90 min. The in situ hybridization and the following photographing were carried out according to the protocols described by Zhong et al. (2024a).

2.4 Karyotype analysis and meiosis observation

Karyotype analysis was performed based on the GISH and FISH signals integrated with relative length and arm ratio, according to the criteria summarized by Li and Chen (1985). For each plant material, 10 cells were selected for the measurement of karyotype data with Auto CAD 2019. Microsoft Office Excel 2016 and Adobe Photoshop 2020 were used for data processing and image postprocessing (Zhong et al., 2024a). The difference on arm ratios of homologous chromosomes between ‘GI’ and ‘GE’ were tested based on one-way ANOVA and LSD test using SPSS 18.0.

Meiosis behaviors were first observed and classified using the PMCs stained with Carbol-Fuchsin. Meanwhile, the behaviors of A- and B-genome sourced chromosomes in meiosis were observed using the PMCs marked with GISH. The number of monads, dyads, triads, and tetrads were counted with PMCs at the tetrad stage. The predicted frequency of 2n gametes were calculated using the following formula by Kondo et al. (2022): Frequency of 2n gamete (%) = (2×dyad + triad)/(monad + 2×dyad + 3×triad + 4×tetrad) × 100%.

3 Results

3.1 Karyotype of ‘GE’ and ‘GI’ revealed by GISH and FISH

The chromosomes of ‘GE’ and ‘GI’ were clearly distinguished into A- and B-genomes, according to the GISH signals with B-genome DNA as probes (Figure 1). FISH signals of 45S and 5S rDNA showed similar distribution patterns in ‘GE’ and ‘GI’. 45S rDNA signals are located on the end of short arm of 3A~5A and 2B~5B chromosomes. The 5S rDNA signals are located near the middle of the short arm of chromosome 3 of both A- and B-genomes (Figures 1, 2).

Figure 1
www.frontiersin.org

Figure 1. Chromosome identification based on genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH) in ‘Golden Era’ and ‘Golden Isles’. (a, b): Composite images with GISH and FISH signals for chromosome identification in ‘Golden Era’ and ‘Golden Isles’, white arrows show the inter-genome translocated fragments; (c, d): DAPI stained chromosomes of ‘Golden Era’ and ‘Golden Isles’; (e, f) 45S rDNA sites (green) recognized by FISH and B genome (red) identified by GISH; (g, h): 5S rDNA sites (red) recognized by FISH.

Figure 2
www.frontiersin.org

Figure 2. Karyograms (a, b) and ideograms (c, d) prepared with results from the genomic in situ hybridization (GISH) and fluorescence in situ hybridization (FISH) in ‘Golden Era’ (a, c) and ‘Golden Isles’ (b, d).

Karyotype analysis based on GISH and FISH showed that ‘GE’ is an allotetraploid (2n=4x=20 = 15m+2sm+3st, AABB) (Figure 2; Supplementary Table S1). Chromosome 4 and 5 are heterozygous in arm ratio. The 41A, 42Aand 42B chromosomes are all median (m) type, which are similar to chromosome 1~3, while 41B chromosome is submedian (sm) type. Chromosome 5 of the two A-genomes and one of the B-genomes are subterminal (st) type, while only the 52B chromosome is sm type. Among the 20 chromosomes of ‘GE’, 4 ones have inter-genome translocation segments, on the long arms of 42A, 41B and 21B chromosomes, and on both long and short arms of the 32B chromosome (Figures 2a, c; Supplementary Table S1). ‘GI’ is a diploid (2n=2x=10 = 8m+1sm+1st, AB), 5B chromosome of which is sm type, while 5A chromosome is st type. Both 1A~4A and 1B~4B chromosomes of ‘GI’ are m type (Figures 2b, d; Supplementary Table S2).

Although the 41B chromosome in ‘GE’ was classified as sm type according to its arm ratio (1.71 ± 0.27), different from other 4B or 4A chromosomes in ‘GE’ (Supplementary Table S1). There was no significant difference on arm ratio between it and the 4B chromosome in ‘GI’ (Supplementary Tables S2, S3). The relative length of inter-genome translocation segments were 2.14 ± 0.28% and 2.15 ± 0.20% in 41B and 42A chromosomes (Supplementary Table S1), respectively. Therefore, the inter-genome translocations on the long arm of 41B and 42A chromosomes in ‘GE’ might be reciprocally translocated from the 4A and 4B chromosomes in ‘GI’, without changing their arm ratios significantly. The arm ratio of 52B chromosome in ‘GE’ is similar to that of 5B chromosome in ‘GI’ (Supplementary Tables S1, S2), which is significantly smaller than that of 51B chromosome in ‘GE’ (Supplementary Table S3), implying that the 52B chromosome might come from the 5B chromosome in ‘GI’. In addition, the other seven chromosomes of ‘GI’ are also highly similar to their homologous chromosomes in ‘GE’, suggesting that 10 of ‘GE’’s 20 chromosomes are likely come from its seed parent, ‘GI’.

3.2 The meiosis behaviors of ‘GI’

The observation on 471 pollen mother cells (PMCs) of ‘GI’ in metaphase I by Carbol-Fuchsin staining indicated that the mean meiosis configuration was 2n=5.8I+1.84II+0.04III+0.01IV, the pairing index was 15.79%, and about 11.06% of the PMCs showed no chromosome pairing (Figure 3). The PMCs bearing polyvalents accounted for about 9.98% of the total number. Monovalents were found in 98.09% of the PMCs (Supplementary Table S4). GISH was used to observe the 198 PMCs in metaphase I. The results showed that the proportion of three types of bivalents (Figure 3), IIAB, IIAA and IIBB, were 92.83%, 1.79% and 5.38%, respectively. Inter-genome translocations were observed in IIAB bivalents. The proportion of rod bivalents with single chiasma was 89.00%, which was much higher than the ring bivalents with double chiasmata, indicating that the probability of simultaneous translocations in both long and short arms was much lower than that of single translocation in only one arm (Figure 4). The proportion of abnormal PMCs in anaphase I was about 75.82%, and the abnormal chromosome behaviors included unequal segregation, lagging chromosomes, chromosome bridges and fragments, and premature separation of sister chromatids (Figure 5). In 15.38% of the PMCs, the chromosomes did not migrate to the poles, resulted in monocytes in telophase I, indicating that a considerable proportion of PMCs did not form reduced nuclei at the end of the meiotic first division.

Figure 3
www.frontiersin.org

Figure 3. Chromosome pairing in metaphase I of ‘Golden Isles’ (2n=2x=10, AB). (a) absence of pairing; (b) 1IIAB+8I; (c) 2IIAB+6I; (d) 3IIAB+4I; (e) 4IIAB+2I; (f) 5IIAB; (g) 1IIIABB+7I; (h) 1IIIAAB+1IIAB+5I; (i) 1IIIBAB+2IIAB+3I; (j) 1IIIABB+2IIAB+1IIAA+1I (red arrow: inter-genome pairing and translocation, yellow arrow head: pairing in A-genome, white arrow head: pairing in B-genome).

Figure 4
www.frontiersin.org

Figure 4. Morphology and constitute of bivalents in metaphase I of ‘Golden Isles’ (2n=2x=10, AB). (a–e): rod bivalents with one chiasma on each, (f, g): ring bivalents with two chiasmata on each; 1: bivalents stained with Carbol-Fuchsin, 2: the chromosomes from B-genome (red) were identified by GISH; inter-genome translocated segments were indicated by arrows (white for A-genome, and red for B-genome).

Figure 5
www.frontiersin.org

Figure 5. Chromosome behaviors in anaphase I and telophase I of ‘Golden Isles’ (2n=2x=10, AB). (a–f): anaphase I, (g–l): telophase I; (a, g): equational separation, (b, h): unequal separation, (c, i): lagging chromosome, (d, j, k): chromosome bridge and fragments, e: premature separation of sister chromatids, (f, l): chromosomes failed in separation; 1: stained with Carbol-Fuchsin solution; 2: the chromosomes from A-genome (blue) were stained with DAPI, and those from B-genome (red) were identified by GISH; blue arrow head: lagging chromosome, yellow arrow head: chromosome bridge, red arrow head: fragment, green arrow head: premature separated sister chromatid.

Among the 655 PMCs in tetrad phase, the proportion of monads, dyads and triads was 1.83%, 29.77% and 12.98%, respectively (Figure 6; Supplementary Table S4), while the tetrads without extra microspores or micronuclei accounted for 22.90%, which meant that the proportion of unreduced gametes was up to 37.79% ((2×29.77 + 12.98)/(1.83 + 2×29.77 + 3×12.98 + 4×22.90) × 100%). In the various spindle orientations of metaphase II, fused spindles accounted for 32.06% and tripolar spindles accounted for 15.68%, indicating that most of the unreduced gametes produced by ‘GI’ might be caused by these two spindle orientation anomalies, and belonging to the first division restitution (FDR) type. Since about 15.38% of the PMCs in the first division of meiosis did not form reduced nuclei, these PMCs would form fused spindles and dyads in the second division of meiosis. Therefore, the unreduced gametes formed through this pathway account for about 16.03% of the total gametes and 42.42% of all unreduced gametes.

Figure 6
www.frontiersin.org

Figure 6. Chromosome behaviors in metaphase II to tetrad stage in ‘Golden Isles’ (2n=2x=10, AB). (a): fused spindles, (b): tripolar spindles, (c): parallel spindles, (d): linear spindles, (e): multipolar spindles, (f): failure to form spindles; (g–l): separation of sister chromatids in telophase II from different types of spindles; (m–r): different number of nuclei resulted from different types of spindles; (s): dyad, (t): triad, (u, v): tetrad, (w): polyad, (x): monad; 1: stained with Carbol-Fuchsin solution, 2: the chromosomes from A-genome (blue) were stained with DAPI, and those from B-genome (red) were identified by GISH.

4 Discussion

4.1 The origin of allotetraploid tree peony

Karyotype analysis based on GISH and FISH showed that 10 of the 20 chromosomes of ‘GE’ (2n=4x=20, AABB) were probably derived from the FDR type unreduced gamete produced by its diploid parent ‘GI’ (2n=2x=10, AB), in which the reciprocal inter-genome translocations were speculated to take place between 4A and 4B chromosomes in their long arms. Meiosis studies verified that ‘GI’ can indeed produce a high proportion (37.97%) of unreduced gametes. These results indicate that 2n gametes produced by diploid inter-subsectional hybrid tree peony (P. × lemoinei) played an important role in the origination of allotetraploid tree peony, similar to that of the allotriploid tree peony (Zhong et al., 2024a). The rest 10 chromosomes of ‘GE’ also contain both A and B genomes, with inter-genome translocations in 2B and 3B chromosomes, implying that the unknown pollen parent of ‘GE’ is also P. × lemoinei.

Meiosis studies have shown that the diploids of P. × lemoinei are prone to gengerate translocations in the long arms, and occasionally in both the long and short arms (Figure 4, unpublished data). Therefore, the pollen parent of ‘GE’ might also be another diploid of P. × lemoinei. However, the fertile 2n (2x) gametes from diploid hybrids are rare in most cases, as a result, the probability of producing tetraploids by the cross of female and male 2n gametes is extremely low (De Storme and Mason, 2014; Ramsey and Schemske, 1998). In contrast, allotetraploids are more fertile and easier to produce 2x gametes (Alix et al., 2017). The earliest fertile individuals in P. × lemoinei can be traced back to the first two F2 plants produced before 1950s, followed by a number of fertile ones obtained in advanced generation hybridizations (Page, 2005; Wister, 1995). The fertility restoration of hybrids is often accompanied by ploidy increasement (Abbott et al., 2013; Alix et al., 2017). Hence, it is speculated that most of these fertile offsprings might be tetraploids (Zhong et al., 2023). So, it is more likely that the unknown pollen parent of ‘GE’ was another allotetraploid tree peony produced earlier than ‘GE’, which was generated in 1980s.

Taken together, ‘GE’ might not be the first allotetraploid tree peony, but they may share the similar genome constitutes. The earliest tetraploid tree peonies are probably those two F2 plants, which might be produced by the hybridization of female and male 2n gametes produced by F1 diploids in P. × lemoinei. The earliest formation of allotetraploid tree peonies is a minimal probability event, which indeed took place only a few decades ago. Once the first allotetraploids are formed, they can hybridize with each other or with diploid F1 hybrids, and it is possible to produce more allotetraploids (Ramsey and Schemske, 1998).

4.2 Inter-subsectional allotetraploid formation vs. intra-subsectional homoploid hybrid speciation in section Moutan

Hybridization is not only regarded as evolutionary noise, but also recognized as a creative evolutionary force in plant speciation and evolution (Mallet, 2007; Soltis and Soltis, 2009). Compared with the frequent allopolyploid speciation, homoploid hybrid speciation is much less common (Long and Rieseberg, 2024). However, in Paeonia section Moutan, allotetraploid formation has only taken place in inter-subsectional hybridization, between P. delavayi from subsection Delavayanae and P. suffruticosa from subsection Vaginatae (Zhong et al., 2024b, 2023). While the hybridization among close relative species within subsection Vaginatae always produce homoploid hybrids or even hybrid species. Homoploid hybridization involving five wild species in subsection Vaginatae that people transplanted to their gardens could have given rise to the diverse cultivars of P. suffruticosa, the cultivated tree peony, which are now grown all over the world (Zhou et al., 2014). In addition, P. × yananensis is documented as a naturally formed homoploid hybrid species, with P. jishanensis as the maternal parent and P. rockii as the paternal parent (Yuan et al., 2010).

Unreduced (2n) gametes are crucial for sexual polyploidization in plants (De Storme and Geelen, 2013). The probability of 2n gamete production is much higher in hybrid plants than in non-hybrid ones (Ramsey and Schemske, 1998). All the species in section Moutan is very stable in ploidy (diploid, 2n=2x=10) (Hong, 2021), suggesting that the extremely rare natural production of 2n gametes might be common in plants of section Moutan. Due to the triploid block, it is difficult to obtain triploid hybrid offsprings by crossing 2n gametes with normal gametes, but easier to obtain tetraploids by crossing between 2n gametes (Ramsey and Schemske, 1998). But when the female 2n gametes artificially induced in P. ostii by high temperature treatment were open pollinated, only a small number of triploid plants, but no tetraploid ones, were obtained (Liu et al., 2023), indicating that it is definitely difficult for tree peonies to produce 2n gametes naturally. Nevertheless, 2n gamete generation have been observable in inter-subsectional hybrid tree peonies (Zhong et al., 2024a, 2019), sometimes in a high frequency (this research).

The genomic differentiation between progenitor taxa influences the likelihood of diploid (homoploid) versus polyploid hybrid speciation, because genetic divergence between parents of allopolyploids is found to be significantly greater than in the case of homoploid hybrid species (Paun et al., 2009). Only when there is a suitable (neither too low nor too high) level of divergence between the parental taxa, they can produce the diploid F1 hybrids which are capable to produce 2n gametes in a high frequency (Sang et al., 2004). Most F1 hybrids observed to produce polyploids were highly sterile (Ramsey and Schemske, 1998). The crossing between inter- or intra-subsectional species of section Moutan both generates homoploid hybrids in F1 generation, but with diverse difficulty. The intra-subsectional hybridization is easy to get viable seeds and fertile hybrids, while the inter-subsectional hybridization only generates few viable seeds and fewer sterile hybrids (Wang et al., 2013), which implying that severe incompatibility exists between inter-subsectional species, instead of intra-subsectional species.

Therefore, the suitable level of divergence between the parents of inter-subsectional hybridization, P. delavayi and P. suffruticosa, which belong to subsection Delavayanae and subsection Vaginatae, respectively, probably promotes the polyploidization of tree peonies by facilitating 2n gamete production.

4.3 Comparison of polyploidization mechanisms between section Moutan and section Paeonia: similarities and differences

Different from the status in section Moutan, polyploid speciation is quite common in section Paeonia. In the 22 tetraploid taxa, including 8 species and 14 subspecies, of section Paeonia, only four taxa are autotetraploids, whereas 18 are allotetraploids. Therein, 8 allotetraploid taxa are formed by hybridization and polyploidization between diploid taxa, while 10 are derived from homoploid tetraploid hybridization (Zhou et al., 2021), implying that the hybridization among diploid taxa play an important role in polyploidization in section Paeonia. Unlike the different level of compatibility between inter- and intra-subsectional hybridization in section Moutan, almost all the diploid species in section Paeonia are intersterile, including those involved in hybridization and polyploidization (Hong, 2021; Saunders and Stebbins, 1938; Zhou et al., 2021). Hence, the mechanism of allotetraploidization in section Paeonia is consistent with that in section Moutan, in which the key point is the hybridization between divergent diploid parents with significant incompatibility. Additionally, there are totally 5 pairs of parents involving 6 diploid taxa, which contribute to the hybridization and polyploidization in section Paeonia (Zhou et al., 2021). The parents of each pair uniformly diverged 18.72~21.61 Mya, approximate with 22.40 Mya, the time of divergence between P. delavayi (subsection Delavayanae) and P. suffruticosa (subsection Vaginatae) in section Moutan (Table 1). Therefore, the divergence level (18.72~ 22.40 Mya) of diploid parents might be the universally fundamental element in hybridization and polyploidization in genus Paeonia.

Table 1
www.frontiersin.org

Table 1. Divergence time of the diploid parents involved in hybridization and polyploidization in section Paeonia (Zhou et al., 2021).

The main difference is that allotetraploidization occurred in nature in section Paeonia, but in gardens in section Moutan. The habitats of species within subsection Vaginatae are completely isolated from those of subsection Delavayanae by the severe ecological shifts of the Hengduan Mountains, Southwestern China (Hong, 2010; Zhou et al., 2014), which might eliminate the possibility of their natural hybridization. In contrast, the diploid parents of allotetraploids in section Paeonia are sympatric (e.g. P. veitchii vs. P. obovata, and P. daurica vs. P. tenuifolia) (Hong, 2010), or might be once sympatric (e.g. P. obovata vs. other species) (Zhou et al., 2021), with opportunities to hybridize in nature. Hence, the allotetraploidization in section Moutan shares the same mechanism with that in section Paeonia, but in a different way to make it realized, that is with the assist from mankind or not.

Furthermore, four autotetraploid taxa are found in section Paeonia, while none in section Moutan (Zhou et al., 2021), suggesting that the probability of 2n gametes production in non-hybrid plants of the two sections might be very different. Some tetraploid plants were obtained from the cross between induced 2n gametes and natural gametes from P. lactiflora, a diploid species of section Paeonia (Zhu et al., 2022), indicating that P. lactiflora can produce some 2n gametes naturally. However, the induced 2n gametes of P. ostii, a diploid species of section Moutan, could only obtain triploids, but not tetraploids, through open pollination (Liu et al., 2023), indicating that P. ostii can hardly produce 2n gametes naturally. This difference might be the main reason for the inability of section Moutan to produce autotetraploids naturally.

4.4 Life history divergence is likely to result in the different patterns of polyploidization, as well as speciation and evolution in genus Paeonia

As summarized above, interparent incompatibility is necessary for allopolyploidization in genus Paeonia. However, intersterility exists not only between species from different clades that diverged earlier than 18.72 Mya both in section Moutan and section Paeonia, but also exits between those from the clades that diverged 9.34~14.08 Mya in section Paeonia (Zhou et al., 2021). While in section Moutan, P. rockii is compatible with other species, including P. jishanensis and P. ostii, that diverged 12.06 Mya from it (Zhou et al., 2021), implying that the level of incompatibility is in a different proportion with genetic divergence in section Paeonia, from that in section Moutan. Moreover, there are much more clades and species (over twice as much), and much wider distribution ranges, in section Paeonia, than in section Moutan, indicating the different tempo of speciation and habitat expansion between them (Figure 7).

Figure 7
www.frontiersin.org

Figure 7. Summary of polyploidization mechanisms in genus Paeonia. The different life history between section Moutan and section Paeonia, which is woody and herbaceous, respectively, could be the main contributor to their diverse level of polyploidization, and different tempo of speciation and evolution.

One of the most divergent traits between section Moutan and section Paeonia is life history, which is woody and herbaceous, respectively. Rates of molecular evolution are linked to life history in flowering plants, which are consistently low in trees and shrubs, as compared with related herbaceous plants (Gaut et al., 2011; Smith and Donoghue, 2008). The positive relationship between species numbers and the rate of neutral molecular evolution have been demonstrated in flowering plants (Barraclough and Savolainen, 2001). Reproductive isolation is positively correlated with genetic distance, and the postzygotic incompatibilities could accumulate and even increase faster than a linear rate with time (Eric et al., 2015). In turn, reproductive isolation can also facilitate the accumulation of genetic differences, with the resulting feedback loop, given enough time, usually leads to complete genetic isolation (Bock et al., 2023; Rieseberg and Willis, 2007). The establishment of reproductive barriers is not only the prerequisite for species formation (Fernández-Mazuecos and Glover, 2017), but also important for allopolyploidization (Ramsey and Schemske, 1998).

In addition, the rate of climatic niche evolution is also growth-form-dependent, namely, woody lineages accumulate fewer changes per million years in climatic niche space, and explore smaller climate space than related herbaceous lineages (Smith and Beaulieu, 2009). Adaptive evolution in plants is largely connected with genomic variations, which could affect genome size variation, gene expression, phenotype variation, and adaptation (Hu et al., 2023). Moreover, ecological adaptation is the major driver of reproductive isolation (Sobel et al., 2010). The Pan-Himalaya is deemed to be the refugium of both woody and herbaceous peonies. There were five vicariance events and 21 dispersal events in the evolutionary history of genus Paeonia (Zhou et al., 2021). Only one vicariance event and five dispersal events are involved in the speciation and evolution of woody peonies. The woody and herbaceous peonies both dispersed into the adjacent areas, including East Asia, in the early dispersal events. But only the herbaceous ones dispersed latter to Central and West Asia, North America, Europe, and North Africa (Hong, 2010; Zhou et al., 2021), where the climatic oscillations and sea level eustacy might have caused isolation and secondary contacts of the previously isolated species, giving the opportunity of hybridization (Mao et al., 2021; Zhou et al., 2021). Repeated cycles of connectivity and isolation may act as drivers of species diversification (Mosbrugger et al., 2018), as well as allopolyploidization (Stebbins, 1985).

Therefore, the different life history between section Moutan and section Paeonia, which is woody and herbaceous, respectively, could be the main contributor to their diverse patterns of polyploidization, and even speciation and evolution (Figure 7).

5 Conclusion

The distant hybridization between intersterile species from different subsections of section Moutan probably promotes the tetraploidization of tree peonies by facilitating 2n gamete production. The mechanism of tetraploidization in section Moutan is highly consistent with that in section Paeonia, but needs the assist from mankind. The divergence of life history between tree peonies and herbaceous peonies is speculated to contribute to the different level of polyploidization, and the distinct tempo of speciation and evolution, between section Moutan and section Paeonia in genus Paeonia.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Author contributions

YZ: Conceptualization, Methodology, Writing – original draft, Writing – review & editing. MD: Methodology, Writing – original draft. RJ: Writing – original draft. FC: Conceptualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No.: 32471937, 31971701), National Key Research and Development Project of China (No.: 2023YFD1200105, 2020YFD1000503), and 5.5 Engineering Research & Innovation Team Project of Beijing Forestry University (No.: BLRC2023A06).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1586225/full#supplementary-material

References

Abbott, R., Albach, D., Ansell, S., Arntzen, J. W., Baird, S. J. E., Bierne, N., et al. (2013). Hybridization and speciation. J. Evolutionary Biol. 26, 229–246. doi: 10.1111/j.1420-9101.2012.02599.x

PubMed Abstract | Crossref Full Text | Google Scholar

Alix, K., Gerard, P. R., Schwarzacher, T., Heslop-Harrison, J. S. (2017). Polyploidy and interspecific hybridization: partners for adaptation, speciation and evolution in plants. Ann. Bot. 120, 183–194. doi: 10.1093/aob/mcx079

PubMed Abstract | Crossref Full Text | Google Scholar

Barraclough, T. G., Savolainen, V. (2001). Evolutionary rates and species diversity in flowering plants. Evolution 55, 677–683. doi: 10.1111/j.0014-3820.2001.tb00803.x

Crossref Full Text | Google Scholar

Bock, D. G., Cai, Z., Elphinstone, C., González-Segovia, E., Hirabayashi, K., Huang, K., et al. (2023). Genomics of plant speciation. Plant Commun. 4, 100599. doi: 10.1016/j.xplc.2023.100599

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Z. J. (2010). Molecular mechanisms of polyploidy and hybrid vigor. Trends Plant Sci. 15, 57–71. doi: 10.1016/j.tplants.2009.12.003

PubMed Abstract | Crossref Full Text | Google Scholar

Cheng, F. Y. (2007). Advances in the breeding of tree peonies and a cultivar system for the cultivar group. Int. J. Plant Breed. 1, 89–104.

Google Scholar

Dang, J., Zhao, Q., Yang, X., Chen, Z., Xiang, S., Liang, G. (2015). A modified method for preparing meiotic chromosomes based on digesting pollen mother cells in suspension. Mol. Cytogenetics 8, 80. doi: 10.1186/s13039-015-0184-x

PubMed Abstract | Crossref Full Text | Google Scholar

De Storme, N., Geelen, D. (2013). Sexual polyploidization in plants – cytological mechanisms and molecular regulation. New Phytol. 198, 670–684. doi: 10.1111/nph.12184

PubMed Abstract | Crossref Full Text | Google Scholar

De Storme, N., Mason, A. (2014). Plant speciation through chromosome instability and ploidy change: cellular mechanisms, molecular factors and evolutionary relevance. Curr. Plant Biol. 1, 10–33. doi: 10.1016/j.cpb.2014.09.002

Crossref Full Text | Google Scholar

Eric, B., Maria Clara, M., Rieseberg, L. H., Daniel, O. B. (2015). The origins of reproductive isolation in plants. New Phytol. 207, 968–984. doi: 10.1111/nph.13424

PubMed Abstract | Crossref Full Text | Google Scholar

Fang, Z., Morrell, P. L. (2016). Domestication: Polyploidy boosts domestication. Nat. Plants 2, 16116. doi: 10.1038/nplants.2016.116

PubMed Abstract | Crossref Full Text | Google Scholar

Ferguson, D., Sang, T. (2001). Speciation through homoploid hybridization between allotetraploids in peonies (Paeonia). Proc. Natl. Acad. Sci. 98, 3915–3919. doi: 10.1073/pnas.061288698

PubMed Abstract | Crossref Full Text | Google Scholar

Fernández-Mazuecos, M., Glover, B. J. (2017). The evo-devo of plant speciation. Nat. Ecol. Evol. 1, 110. doi: 10.1038/s41559-017-0110

PubMed Abstract | Crossref Full Text | Google Scholar

Gaut, B., Yang, L., Takuno, S., Eguiarte, L. E. (2011). The patterns and causes of variation in plant nucleotide substitution rates. Annu. Rev. Ecology Evolution Systematics 42, 245–266. doi: 10.1146/annurev-ecolsys-102710-145119

Crossref Full Text | Google Scholar

Hong, D. Y. (2010). Peonies of the World: Taxonomy and Phytogeography. (UK: Kew Publishing).

Google Scholar

Hong, D. Y. (2021). Peonies of the World: Phylogeny and Evolution. (UK: Kew Publishing).

Google Scholar

Hu, Y., Wang, X., Xu, Y., Yang, H., Tong, Z., Tian, R., et al. (2023). Molecular mechanisms of adaptive evolution in wild animals and plants. Sci. China. Life Sci. 66, 453–495. doi: 10.1007/s11427-022-2233-x

PubMed Abstract | Crossref Full Text | Google Scholar

Jiao, Y., Wickett, N. J., Ayyampalayam, S., Chanderbali, A. S., Landherr, L., Ralph, P. E., et al. (2011). Ancestral polyploidy in seed plants and angiosperms. Nature 473, 97–113. doi: 10.1038/nature09916

PubMed Abstract | Crossref Full Text | Google Scholar

Kang, X., Wei, H. (2022). Breeding polyploid Populus: progress and perspective. Forestry Res. 1, 4–10. doi: 10.48130/FR-2022-0004

PubMed Abstract | Crossref Full Text | Google Scholar

Kondo, H., Deguchi, A., Kikuchi, S., Miyoshi, K. (2022). Two pathways of 2n gamete formation and differences in the frequencies of 2n gametes between wild species and interspecific hybrids. Plant Cell Rep. 41, 2187–2200. doi: 10.1007/s00299-022-02915-5

PubMed Abstract | Crossref Full Text | Google Scholar

Kreiner, J. M., Kron, P., Husband, B. C. (2017). Evolutionary dynamics of unreduced gametes. Trends Genet. 33, 583–593. doi: 10.1016/j.tig.2017.06.009

PubMed Abstract | Crossref Full Text | Google Scholar

Li, M. X., Chen, R. Y. (1985). A suggestion on the standardization of karyotype analysis in plants. J. Wuhan Botanical Res. 3, 297–302.

Google Scholar

Li, M. X., Zhang, X. F. (1982). Cytogenetic observation of triploid tree peony. Hereditas 4, 19–21.

Google Scholar

Liu, C. Y., Peng, C. F., Cheng, S. P., Yao, P. Q., Geng, X. N., Xie, L. H. (2023). Creation of triploid germplasm in Paeonia ostii ‘Fengdan’ through 2n female gametes inducing with high temperature treatment. Acta Hortic. Sin. 50, 1455–1466. doi: 10.16420/j.issn.0513-353x.2022-0699

Crossref Full Text | Google Scholar

Long, Z., Rieseberg, L. H. (2024). Documenting homoploid hybrid speciation. Mol. Ecol., e17412. doi: 10.1111/mec.17412

PubMed Abstract | Crossref Full Text | Google Scholar

Mallet, J. (2007). Hybrid speciation. Nature 446, 279–283. doi: 10.1038/nature05706

PubMed Abstract | Crossref Full Text | Google Scholar

Mao, K. S., Wang, Y., Liu, J. Q. (2021). Evolutionary origin of species diversity on the Qinghai–Tibet Plateau. J. Systematics Evol. 59, 1142–1158. doi: 10.1111/jse.12809

Crossref Full Text | Google Scholar

Mosbrugger, V., Favre, A., Muellner-Riehl, A., Päckert, M., Mulch, A. (2018). “Cenozoic evolution of geo-biodiversity in the Tibeto-Himalayan region,” in Mountains, Climate and Biodiversity. Eds. Hoorn, C., Perrigo, A., Antonelli, A. (Chichester, UK: Wiley-Blackwell), 429–449.

Google Scholar

Page, M. (2005). The Gardener’s Peony: Herbaceous and Tree Peonies. (Portland • Cambridge: Timber Press).

Google Scholar

Paun, O., Forest, F., Fay, M. F., Chase, M. W. (2009). Hybrid speciation in angiosperms: parental divergence drives ploidy. New Phytol. 182, 507–518. doi: 10.1111/j.1469-8137.2009.02767.x

PubMed Abstract | Crossref Full Text | Google Scholar

Ramsey, J., Schemske, D. W. (1998). Pathways, mechanisms, and rates of polyploid formation in flowering plants. Annu. Rev. Ecol. Systematics 29, 467–501. doi: 10.1146/annurev.ecolsys.29.1.467

Crossref Full Text | Google Scholar

Rieseberg, L. H., Willis, J. H. (2007). Plant speciation. Science 317, 910–914. doi: 10.1126/science.1137729

PubMed Abstract | Crossref Full Text | Google Scholar

Salman-Minkov, A., Sabath, N., Mayrose, I. (2016). Whole-genome duplication as a key factor in crop domestication. Nat. Plants 2, 16115. doi: 10.1038/NPLANTS.2016.115

PubMed Abstract | Crossref Full Text | Google Scholar

Sang, T., Pan, J., Zhang, D. M., Ferguson, D., Wang, C., Pan, K. Y., et al. (2004). Origins of polyploids: an example from peonies (Paeonia) and a model for angiosperms. Biol. J. Linn. Soc. 82, 561–571. doi: 10.1111/j.1095-8312.2004.00341.x

Crossref Full Text | Google Scholar

Saunders, A. P., Stebbins, G. L. (1938). Cytogenetic studies in Paeonia I. the compatibility of the species and the appearance of the hybrids. Genetics 23, 65–82. doi: 10.1093/genetics/23.1.65

PubMed Abstract | Crossref Full Text | Google Scholar

Smith, S. A., Beaulieu, J. M. (2009). Life history influences rates of climatic niche evolution in flowering plants. Proc. R. Soc. B 276, 4345–4352. doi: 10.1098/rspb.2009.1176

PubMed Abstract | Crossref Full Text | Google Scholar

Smith, S. A., Donoghue, M. J. (2008). Rates of molecular evolution are linked to life history in flowering plants. Science 322, 86–89. doi: 10.1126/science.1163197

PubMed Abstract | Crossref Full Text | Google Scholar

Sobel, J. M., Chen, G. F., Watt, L. R., Schemske, D. W. (2010). The biology of speciation. Evolution 64, 295–315. doi: 10.1111/j.1558-5646.2009.00877.x

PubMed Abstract | Crossref Full Text | Google Scholar

Soltis, D. E., Albert, V. A., Leebensmack, J., Bell, C. D., Paterson, A. H., Zheng, C., et al. (2009). Polyploidy and angiosperm diversification. Am. J. Bot. 96, 336. doi: 10.3732/ajb.0800079

PubMed Abstract | Crossref Full Text | Google Scholar

Soltis, P. S., Soltis, D. E. (2009). The role of hybridization in plant speciation. Annu. Rev. Plant Biol. 60, 561–588. doi: 10.1146/annurev.arplant.043008.092039

PubMed Abstract | Crossref Full Text | Google Scholar

Stebbins, G. L. (1985). Polyploidy, hybridization, and the invasion of new habitats. Ann. Missouri Botanical Garden 72, 824–832. doi: 10.2307/2399224

Crossref Full Text | Google Scholar

Wang, L. Y., Yuan, T., Wang, F., Li, Q. D. (2013). The ex-situ conservation centre of Chinese Paeoniceae wild species and cultivation of new cultivars. (Beijing, China: China Forestry Publishing House).

Google Scholar

Wister, J. C. (1995). The Peonies. Hopkins: The American Peony Society.

Google Scholar

Younis, A., Hwang, Y., Lim, K. (2014). Exploitation of induced 2n-gametes for plant breeding. Plant Cell Rep. 33, 215–223. doi: 10.1007/s00299-013-1534-y

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, J. H., Cheng, F. Y., Zhou, S. L. (2010). Hybrid origin of Paeonia × yananensis revealed by microsatellite markers, chloroplast gene sequences, and morphological characteristics. Int. J. Plant Sci. 171, 409–420. doi: 10.1086/651228

Crossref Full Text | Google Scholar

Zhang, K., Wang, X., Cheng, F. (2019). Plant polyploidy: origin, evolution, and its influence on crop domestication. Hortic. Plant J. 5, 231–239. doi: 10.1016/j.hpj.2019.11.003

Crossref Full Text | Google Scholar

Zhong, Y., Du, M. J., Cheng, F. Y., Li, K. Y., Ji, R. Z., Cui, J. (2023). Karyotype and fertility analysis of allotetraploids in tree peony (Paeonia section Moutan). Acta Hortic. Sin. 50, 1444–1454. doi: 10.16420/j.issn.0513-353x.2022-0413

Crossref Full Text | Google Scholar

Zhong, Y., Du, M., Ji, R., Cheng, F. (2024a). Cytological origination of the first found allotriploid tree peony, Paeonia × lemoinei ‘Oukan’ (AAB), reveled by molecular karyotype comparation. Scientia Hortic. 324, 112563. doi: 10.1016/j.scienta.2023.112563

Crossref Full Text | Google Scholar

Zhong, Y., Du, M., Ji, R., Rong, Q., Guan, Z., Cheng, F. (2024b). Cytogenetic analysis reveals a mechanism of compatibility in distant hybridization between tree peony and herbaceous peony. Euphytica 220, 66. doi: 10.1007/s10681-024-03326-1

Crossref Full Text | Google Scholar

Zhong, Y., Du, M. J., Liu, Y. X., Cheng, F. Y. (2019). Karyotype analysis and meiotic chromosome behavior observation on Paeonia × lemoinei ‘High Noon’. J. Beijing Forestry Univ. 41, 68–73. doi: 10.13332/j.1000-1522.20180381

Crossref Full Text | Google Scholar

Zhou, S. L., Xu, C., Liu, J., Yu, Y., Wu, P., Cheng, T., et al. (2021). Out of the Pan-Himalaya: evolutionary history of the Paeoniaceae revealed by phylogenomics. J. Systematics Evol. 59, 1170–1182. doi: 10.1111/jse.12688

Crossref Full Text | Google Scholar

Zhou, S. L., Zou, X. H., Zhou, Z. Q., Liu, J., Xu, C., Yu, J., et al. (2014). Multiple species of wild tree peonies gave rise to the ‘king of flowers’. Paeonia suffruticosa Andrews. Proc. R. Soc. B: Biol. Sci. 281, 20141687. doi: 10.1098/rspb.2014.1687

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, W., Dong, Z., Chen, X., Cao, J., Zhang, W., Sun, R., et al. (2022). Induction of 2n pollen by colchicine during microsporogenesis to produce polyploids in herbaceous peony (Paeonia lactiflora Pall.). Scientia Hortic. 304, 111264. doi: 10.1016/j.scienta.2022.111264

Crossref Full Text | Google Scholar

Keywords: karyotype, meiosis, GISH, FISH, hybrid, polypoid

Citation: Zhong Y, Du M-j, Ji R-z and Cheng F-y (2025) The mechanism of tetraploidization in tree peony, and its implications for speciation and evolution of genus Paeonia L.. Front. Plant Sci. 16:1586225. doi: 10.3389/fpls.2025.1586225

Received: 02 March 2025; Accepted: 08 April 2025;
Published: 12 May 2025.

Edited by:

Rebecca Grumet, Michigan State University, United States

Reviewed by:

Kyong-Sook Chung, Jungwon University, Republic of Korea
Jian Yu Wang, North China University of Science and Technology, China

Copyright © 2025 Zhong, Du, Ji and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yuan Zhong, emhvbmd5dWFuYmpmdUAxMjYuY29t; Fang-yun Cheng, Y2hlbmdmeThAMjYzLm5ldA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.