Skip to main content

Does the evolution of micromorphology accompany chromosomal changes on dysploid and polyploid levels in the Barnardia japonica complex (Hyacinthaceae)?

Abstract

Background

Chromosome number and genome size changes via dysploidy and polyploidy accompany plant diversification and speciation. Such changes often impact also morphological characters. An excellent system to address the questions of how extensive and structured chromosomal changes within one species complex affect the phenotype is the monocot species complex of Barnardia japonica. This taxon contains two well established and distinct diploid cytotypes differing in base chromosome numbers (AA: x = 8, BB: x = 9) and their allopolyploid derivatives on several ploidy levels (from 3x to 6x). This extensive and structured genomic variation, however, is not mirrored by gross morphological differentiation.

Results

The current study aims to analyze the correlations between the changes of chromosome numbers and genome sizes with palynological and leaf micromorphological characters in diploids and selected allopolyploids of the B. japonica complex. The chromosome numbers varied from 2n = 16 and 18 (2n = 25 with the presence of supernumerary B chromosomes), and from 2n = 26 to 51 in polyploids on four different ploidy levels (3x, 4x, 5x, and 6x). Despite additive chromosome numbers compared to diploid parental cytotypes, all polyploid cytotypes have experienced genome downsizing. Analyses of leaf micromorphological characters did not reveal any diagnostic traits that could be specifically assigned to individual cytotypes. The variation of pollen grain sizes correlated positively with ploidy levels.

Conclusions

This study clearly demonstrates that karyotype and genome size differentiation does not have to be correlated with morphological differentiation of cytotypes.

Peer Review reports

Introduction

Chromosome number changes comprising polyploidy, and dysploidy, play an important role in plant genome diversification and have thus been frequently analyzed to better understand the angiosperm evolution [1, 2]. Polyploidy (i.e., the whole-genome duplication) including both auto- and allopolyploidy is one of the major processes involved in the diversification and speciation of plants [2,3,4,5,6]. Autopolyploidy has received relatively little attention due to high levels of morphological similarity between autopolyploids and their parental diploid taxa [7,8,9]. In contrast, allopolyploids have been studied extensively in several plant groups with focus on changes of genome sizes and karyotypes compared to their parental taxa (e.g., Brassica L. [10]; Melampodium L. [11, 12]; Nicotiana L. [13]; Prospero Salisb. [14]; Tragopogon L. [15]). Although polyploidy is recognized as process leading to plant diversification in natural populations, comparative analyses of the evolution of polyploid genomes are rather scarce and limited to several plant groups [11, 12, 14, 16,17,18].

Two mechanisms result in the stepwise change of the haploid chromosome number among related species [1, 19], aneuploidy and dysploidy. While aneuploidy is less relevant to diversification and evolution as it introduces immediate genetic imbalance by addition or removal of single chromosomes, dysploidy contributes to evolution by changing haploid chromosome numbers via chromosomal rearrangements and keeping changes of DNA amount to minimum [1, 2]. Dysploidy and auto- and allopolyploidy, with occasional indication of aneuploidy on higher ploidy levels within species complexes have been reported in many species in family Hyacinthaceae Batsch ex Borkh. (also treated as Asparagaceae Juss.; [14, 20,21,22,23]). The genus Barnardia Lindl. comprises two geographically disjunct species groups: B. numidica (Poir.) Speta and the species complex referred to as B. japonica (Thunb.) Schult. & Schult.f. [24,25,26]. Barnardia numidica is widespread in the Balearic Islands and North-West Africa (Algeria, Tunisia, and Libya) and is known only as diploid with 2n = 18 (x = 9; [27]). In contrast, the B. japonica complex is found in Eastern Asia and exhibits a spectacular array of chromosomal variation with two different base chromosome numbers (x = 8 and 9), an extensive range of polyploids (3x, 4x, 5x and 6x), various types of chromosomal polymorphisms including the presence of B-chromosomes [28, 29], as well as genome size variation [22]. Molecular phylogenetic studies of plastid and nuclear DNA sequence data have revealed that the B. japonica complex is the most basal clade in the family of Hyacinthaceae [30, 31].

Two distinct and stable diploid cytotypes with base chromosome numbers of either x = 8 (cytotype/genome A) and x = 9 (cytotype/genome B) have been described in the B. japonica complex [25, 28, 32]. Diploids of AA (2n = 2x = 16) and BB (2n = 2x = 18) cytotypes are known to hybridize and to form allopolyploids in natural populations. Diploid homoploid hybrids AB (2n = 2x = 17), as well as a myriad of polyploids of various genomic compositions have been reported (ABB, 2n = 3x = 26; BBB, 2n = 3x = 27; AAAA, 2n = 4x = 32; AABB, 2n = 4x = 34; ABBB, 2n = 4x = 35; BBBB, 2n = 4x = 36; AABBB, 2n = 5x = 43; AAABBB, 2n = 6x = 51; [24, 25, 33, 34]). The diploid AA cytotype is found in China and Korea covering nearly the entire distribution range of the B. japonica complex, with the exception of Japan [24, 25]. In contrast, the BB cytotype is geographically more restricted to a narrow area of eastern central China, Jeju Island in Korea, and Japan [35]. It is commonly accepted to refer to the hybrids between different, karyotypically well differentiated cytotypes (particularly in monocots), as allopolyploids, despite the fact that the diploid cytotypes are all part of the same species complex. We are, therefore, going to use the same approach. Allotetraploids of genomic constitution AABB are frequently found throughout the complex distribution range [25, 33, 36]. Cytogenomic analyses of AABB allotetraploids revealed that the parental subgenomes are retained without any rearrangements except for the loss of NOR (e.g., nucleolar organizer region) in the parental A-genome [25, 37]. Supernumerary genetic material, either as B-chromosomes (Bs) or supernumerary chromosomal segments (SCSs) that are physically integrated into the standard chromosome complement, is frequently found in various plant groups, particularly in monocots, with hot spots in Liliales and Commelinales [38,39,40,41,42,43]. The SCSs can be located either interstitially or terminally and easily identified in homologous chromosomes in the heterozygous condition due to chromosome length differences [2, 14]. Both SCSs and Bs are often, but not always, heterochromatic. B chromosomes are also dispensable components of the genomes and usually behave as selfish genetic elements [2, 40]. In plants, frequency of both SCSs and Bs in diploids and polyploids are similar [44, 45] but tends to be higher in taxa with large genome sizes [39, 46]. Although the presence of Bs have been well documented in the B. japonica complex [28, 29, 33], the occurrence of SCSs has yet to be shown [29, 33].

The plants of the genus Barnardia possess rather large bulbs (over 4 cm in diameter) with several brown scales, an elongated raceme with 30–70 purple flowers, white or pale pink tepals including six stamens, and a central pistil [27, 31, 47]. The micromorphology of leaf epidermis and characters of pollen grains are often used as a diagnostic character to establish relationships and to define major groups in family Hyacinthaceae [48,49,50] as well as other economically important plant families such as Lamiaceae Martinov [51], Polygonaceae Juss. [52], Melastomataceae Juss. [53], Asteraceae Giseke [54], Poaceae Barnhart [55], and Iridaceae Juss. [56,57,58]. Specifically, wax striation and size of stomatal complexes are considered to be important diagnostic value in Hyacinthaceae [48]. Few studies to date have addressed the question of evolution of those characters in relations to extensive dysploidy and polyploidy within species complexes. Similarly, no morphological diagnostic characters have so far been found to allow for identification of the two diploid cytotypes and their hybrid/polyploid descendants in the B. japonica complex [25]. Thus, the current species concept in the B. japonica complex is mostly based on classical karyotaxonomic studies with two diploid cytotypes differing in base chromosome number and a myriad of resulting polyploid cytotypes [22, 24,25,26, 28, 35, 59], while the underlying patterns of genetic and micromorphological differentiation are still unknown.

Thus, the main goals of this study are (1) to establish chromosome numbers, karyotype structures, and genome sizes of 131 diploid and polyploid individuals of the B. japonica complex collected from natural populations in South Korea; (2) to analyze the leaf and pollen micromorphology in diploid and polyploid cytotypes of various genomic constitution, and (3) to correlate the variation of the analyzed characters with structured karyotypic variation in an attempt to identify new diagnostic characters or character combinations for cytotype identification in the Barnardia japonica complex.

Materials and methods

Sampling, chromosome counts, and karyotype analysis

Bulbs of the Barnardia japonica complex were collected in the field between 2017 and 2022 in South Korea (Fig. 1) and cultivated at Chungnam National University. Sampling of diploid and polyploid cytotypes was guided by earlier studies [25] with the aim to collect most of the ploidy levels reported. Thirty-four diploid individuals representing two genomes differing in base chromosome numbers (A genome: x = 8 and B genome: x = 9) and 97 polyploid individuals representing four ploidy levels (3x, 4x, 5x and 6x) were collected from 34 natural populations (Table 1).

Fig. 1
figure 1

Populations of the Barnardia japonica complex analysed in the current study (Table 1). (A) Map with the locations of sampled diploid and polyploid populations; (B) Proportions of diploid and polyploid cytotypes among investigated 131 individuals

Table 1 Individuals of the Barnardia japonica complex analyzed in this study with detailed voucher information

Chromosome numbers and karyotypes were determined using standard Feulgen staining [60]. In total, chromosome numbers of 131 collected plants of the B. japonica complex were established. Twenty-two plants of representing all cytotype were karyotyped to aid identification of their genomic composition ( [25, 28, 34]; Table 1 and Table S1). For cytological investigations, actively growing root meristems were pretreated with 0.05% aqueous solution of colchicine for 4 h at room temperature, fixed in ethanol : acetic acid (3:1), and stored at − 20℃ until use. Root tips were hydrolyzed in 5 N HCl (VWR; Vienna, Austria) at room temperature for 30 min and stained with Schiff’s reagent (Sigma, Vienna, Austria) for 60 min in darkness. Squash preparations were made in a drop of 60% acetic acid. At least three well-spread and complete metaphase plates with a medium degree of chromosome condensation were chosen for karyotyping for each individual. Chromosomes were cut out and karyotypes arranged using Adobe Photoshop CS6. Chromosomes were measured using Micromeasure ver. 3.3 (< www.colostate.edu/Depts/Biology/MicroMeasure/>) following [61]. Chromosome arm lengths and total chromosome lengths were measured in at least three chromosomal spreads per individual (unless otherwise indicated; Table S1).

Genome size measurements

Genome sizes of 42 selected individuals of the B. japonica complex were measured using flow cytometry with Pisum sativum “Kleine Rheinländerin” (4.42 pg/1 C) for diploids and Vicia faba L. “Inovec” (13.45 pg/1 C) for polyploids as internal standards [62,63,64,65]. The methodology used for the measurement of the genome sizes followed [64, 66]. Measurements were performed using Sysmex CyFlow Cytometer (Sysmex Partec GmbH, Görlitz, Germany) and 1 C values were calculated following [63]. The CVs (coefficient values) of all measurements were usually lower than 3% and never exceeded 5% [66]. To test whether the genome sizes of each cytotype with and without B chromosomes was significantly different, the Mann-Whitney test was performed in GraphPad Prism version 9.

Leaf and pollen micromorphological traits

Fresh leaf and pollen samples of all analysed individuals representing all cytotypes of the Barnardia japonica complex were preserved in formalin-acetic-acid-alcohol (FAA). Leaf samples were inspected under Nikon SMZ1500 stereomicroscope (Nikon, Japan) and only fully mature leaves were selected for further leaf epidermis analyses (Table 1). Tissue samples were examined using a light microscopy (LM; BX53F, Olympus, Japan) and scanning electron microscopy (SEM; Hitachi E-1010, Japan) following [64]. The micromorphological variation of abaxial (AB) and adaxial (AD) epidermal surfaces of the leaf including epidermal cells and stomatal complexes were analyzed for 13 individuals (at least one individual per cytotype). The guard cells of leaf and pollen grain size were measured from 30 samples per each cytotype (Table S2). Pollen micromorphological characters (i.e., pollen shape and size, exine sculpture) and pollen viability of nine individuals were analyzed using pollen from fertile anthers of flowers collected in the field. Pollen grains were stained using aniline blue dye solution, which only stained fertile pollen grains [67, 68]. Pollen morphological description and terminology followed [69]. Pollen shape and size including the detailed pollen exine ornamentation and the sculpture of apertures were analyzed using SEM following [56, 70].

To test whether the pollen size along the long axis and ploidy level are correlated, Pearson correlation coefficients were computed, and a Mann-Whitney test was conducted to test if the pollen size measurements (long axis) in diploids significantly differ from those in polyploids (GraphPad Prism version 9).

Results

Chromosome numbers, karyotypes and genome size variation

Analyses of 131 individuals from 34 populations of the chromosomally hypervariable Barnardia japonica complex (Table 1) revealed the presence of both diploid and polyploid individuals. Chromosome numbers in the B. japonica complex varied from 2n = 16 to 25 in diploids due to the presence or regular A-chromosomes (2n = 16 or 2n = 18) accompanied by one to seven B-chromosomes and from 2n = 26 to 51 in polyploids comprising A-chromosomal sets and one to six Bs as well as two aneuploid individuals (Figs. 1 and 2; Table 1). Five ploidy levels were found including diploids with two distinct base chromosome numbers (x = 8 or 9; Fig. 2A, F), triploids (2n = 26; Fig. 2L), tetraploids (2n = 34, 35; Fig. 2M–T), pentaploids (2n = 42, 43; Fig. 2U–Y), and hexaploids (2n = 51; Fig. 2Z). Allotetraploid individuals were found most frequently (78 individuals, 59%), whereas allotri- and allohexaploids (one individual each, 1%) were found only sporadically (Figs. 1 and 2). B-chromosomes, easily distinguishable from the A-complement chromosomes by their smaller size (Fig. 2), were frequently found in both diploids (Fig. 2B–E, G–K) and polyploids (Fig. 2N–O, Q–T, V, X–Y) of the B. japonica complex (Table 1). Supernumerary chromosomal segments (SCSs), physically integrated into the standard chromosome complement, were occasionally found in polyploids (Fig. 2S, W), whereas aneuploidy was encountered in two individuals (Fig. 2V; Table 1).

Fig. 2
figure 2

Mitotic metaphase chromosomes of diploid and polyploid cytotypes of the Barnardia japonica complex. Black and white arrows indicate B-chromosomes and supernumerary chromosomal segments (SCSs), respectively. Scale bar = 10 μm

Representative karyotypes of all eight cytotypes found among analyzed plants: two diploids (i.e., AA and BB), one triploid (i.e., ABB), two tetraploids (i.e., ABBB and AABB), two pentaploids (i.e., AABBB and AAABB) and one hexaploid (i.e., AAABBB) are presented in Fig. 3. Each polyploid cytotype was easily distinguishable due to different karyotype morphology of the parental diploid cytotypes A and B (Fig. 3). Karyotypes of diploid AA and BB cytotypes were easily distinguishable due to the presence of larger chromosomes in AA (4.00–11.52 μm) compared to BB (2.35–7.45 μm) cytotype (Table S1). The haploid chromosome set of AA cytotype was composed of a large metacentric chromosome (A1), the longest subtelocentric chromosome (A2), five subtelocentric chromosomes (A3–A7) and a small metacentric chromosome (A8) (Fig. 3). Haploid chromosome set of BB cytotype possessed a subtelocentric chromosome with a satellite at the short arm (B1), four subtelocentric chromosomes (B2–B5) and four small meta- or submetacentric chromosomes (B6–B9) (Fig. 3). The total haploid chromosome length (HCL) of BB cytotype was lower than that of AA cytotype (42.09 vs. 53.22 μm; Table S1). Karyotype lengths of allopolyploids were not additive compared to their diploid progenitors and ranged from 61.85 μm in allotriploid (expected HKL: ca. 68 μm) to 119.20 μm in allohexaploid (expected HKL: ca. 142 μm; Table S1). Karyotype morphologies of allopolyploids (ABB, ABBB, AABB, AABBB, AAABB, and AAABBB) were additive compared to their diploid progenitors (Fig. 3). Additional genetic materials including both B-chromosomes (one to seven Bs) and supernumerary chromosomal segments (SCSs) have been found in three individuals (Figs. 2 and 4). B-chromosomes were present in both diploids and polyploids. All SCSs were found exclusively in the A subgenome of allopolyploids and were located distally, either within the long or short arms of chromosomes 1, 6, and/or 8 (Figs. 2 and 4).

Fig. 3
figure 3

Karyotypes of diploids and polyploids in the Barnardia japonica complex. Scale bar = 10 μm

Fig. 4
figure 4

Karyotypes of the selected individuals carrying B-chromosomes (boxes) and supernumerary chromosomal segments (black arrows) in the Barnardia japonica complex. Scale bar = 10 μm

The DNA amounts (1 C value) of all individuals ranged from 9.05 to 23.23 pg, indicating 2.5-fold difference (Fig. 5A). The coefficient of variation (CV) for the internal standard and sample peaks ranged from 0.92 to 4.49%. The 1 C values were relatively similar with some variation within all cytotypes (Fig. 5B; Table 1). Variation among individuals without B-chromosomes ranged from 9.05 to 22.95 pg and among bulbs with at least one B-chromosome from 9.13 to 23.23 pg (Table 1). Diploid AA cytotype had average genome size of 1 C = 12.23 pg (11.87–12.46 pg), diploid BB cytotype had average 1 C value of 9.17 pg (9.05–9.28 pg), tetraploid ABBB (2n = 35) had 1 C of 18.84 pg (18.79–18.90 pg), tetraploid AABB (2n = 34) had 1 C of 18.76 pg (17.05–20.26 pg), and pentaploid AABBB (2n = 43) had genome size of 22.55 pg (22.03–22.95 pg) (Table 1). The genome sizes of each cytotype with and without B chromosome were not significantly different (AA, P = 0.7891; BB, P = 0.597; ABBB, P = 0.5221; AABBB, P = 0.4514) regardless of the number of Bs: AA + 5Bs (12.28 pg), BB + 1B (9.30 pg), BB + 3Bs (9.13–9.52 pg), ABBB + 1B (19.02 pg), ABBB + 2Bs (18.67 pg), AABB + 1B (20.32 pg), AABB + 2Bs (20.03–20.19 pg), AABB + 6Bs (21.38 pg), and AABBB + 3Bs (22.42–23.23 pg), except for AABB and AABB with B chromosomes (P < 0.0001) (Fig. 5; Table 1). The 1 C values were cytotype-specific in diploids, but not in polyploids (Fig. 5).

Fig. 5
figure 5

Fluorescence histogram (A) and genome size variation (B) among analysed Barnardia japonica complex individuals (Table 1). Red boxes in (B) indicate the individuals carrying B-chromosomes. For collection numbers please refer to Table 1

Leaf and pollen micromorphological characters

The pollen morphological characters including viability, size, apertures and exine ornamentation have been reported herein for the first time for the different ploidy levels (i.e., 2x, 3x, 4x, 5x and 6x) of B. japonica (Figs. 6 and 7, Fig. S2). Pollen grains of all cytotypes were fertile, although sterile pollen grains were frequently found at odd-ploidy levels (ABB and AABBB). The pollen size ranged from 35.13 × 25.98 μm (long × short axis) in diploid BB to 50.33 × 31.99 μm in pentaploid AABBB indicating nearly 1.5-fold difference (Fig. 6). The pollen size and ploidy levels were positively correlated (r = 0.78, P < 0.0001). Regardless of pollen size variation among the all cytotypes, the exine ornamentations were all reticulate-perforate, and perforations were the major ornamentation in sulcus margin (Fig. 7D–I). The pollen size was significantly different between diploids and polyploids (P < 0.0001). The pollen apertures of all cytotypes were consistently monosulcate with granular membranes possessing tiny perforations (Fig. S2).

Fig. 6
figure 6

Pollen grain size variation among the diploid and polyploid cytotypes. For collection numbers please refer to Table 1

Fig. 7
figure 7

SEM (scanning electron microscope) micrographs of leaf and pollen micromorphological traits of diploids (AA and BB) and AABB allotetraploids of the Barnardia japonica complex. (AC) Stomatal complex of leaf surface. (DI) Pollen grains and exine ornamentation. (A, D, G) AA cytotype, (B, E, H) BB cytotype, (C, F, I) AABB cytotype

Sixteen individuals were selected for detailed analyses of leaf surface characters using SEM: four AA diploids, one BB diploid, one ABB triploid, three allotetraploids AABB, two pentaploids (one of each: AABBB and AAABB) and one AAABBB hexaploid, as well as four plants carrying B-chromosomes (Fig. 7A–C; Table 1). The leaf surface of both diploids and polyploids of the Barnardia japonica complex lacked trichomes (Fig. S1). AA diploids possessed corrugated cuticle (Fig. 7A, Fig. S1) whereas both BB diploid (Fig. 7B) and AABB allotetraploid (Fig. 7C) had slightly ribbed thickenings of leaf epidermis (Fig. S1). Slightly sunken anomocytic stomata were observed in all the cytotypes and no significant guard cell size variation between abaxial and adaxial surface of the leaves was found in all investigated bulbs of the B. japonica complex (Table S2).

Discussion

High levels of karyotype and genome size diversity in the Barnardia japonica complex

The current study reports for the first time comprehensive analysis of chromosome numbers and genome size variation of all hitherto known cytotypes in the perennial bulbous Korean B. japonica complex. The results of the extensive chromosomal survey were largely in agreement with previous reports from Chinese, Japanese, and Korean populations [24, 25, 71]. Genome sizes reported here, however, deviated from those reported by Shibata et al. [22]. The discrepancies (i.e., quite consistent and approximately 2 pg/1 C deviation) are likely to represent methodological problems, as the previous genome size values have been obtained via flow cytometry with less accurate staining of the DNA with DAPI that binds preferentially to AT-rich regions [72] as well as using Allium fistulosum as an internal standard [73]. It has been shown that DAPI staining may lead to larger errors in nuclear DNA content evaluation [62], and therefore propidium iodide is used as a standard dye for reliable estimation of genome size in plants using flow cytometry [74,75,76].

In contrast to the gross morphological uniformity [77], the B. japonica complex exhibits remarkable chromosome number and ploidy level variations with occasionally found mixed-ploidy populations, as also reported earlier [25]. Despite the considerable range of chromosome numbers and genome sizes observed in the B. japonica complex, this variation is well structured. AA diploids have consistently larger genome sizes than BB diploids (P < 0.0001) and there is a positive correlation between genome sizes and ploidy levels. The genome sizes measured in the current study have also provided evidence of genome downsizing in allopolyploids on all ploidy levels (3x, 4x, 5x, and 6x), as frequently reported in other plant groups [78,79,80]. The changes in genome size may be caused by the independent accumulation or reduction of repetitive DNA amounts (e.g., satellite DNAs and/or transposable elements), which contribute to both genome size increases and genome downsizing, which seems to be a general trend in polyploids [81,82,83]. The repetitive DNA composition of genomes in polyploids has often been shown to be fast-evolving compared to their diploid progenitors. This is often related to the process of genome diploidization [82]. Diploidization processes transform the polyploid genomes back into functional diploids through chromosomal rearrangements, genome downsizing, and gene loss or silencing [84,85,86]. There is not experimental evidence so far of the processes involved in the genome diploidization in the polyploids of the B. japonica complex. However, numerical and structural changes of chromosomes, including genome downsizing, are well studied in allopolyploid genomes as shown in Brassica [10], Melampodium [11, 12], Nicotiana [13], Prospero [14], or Tragopogon [15]. Further studies employing more sensitive techniques, including molecular phylogenetic analysis using both cpDNA and nrDNA sequences at the populational level, FISH (fluorescence in situ hybridization), and comparative analyses of the repetitive DNA fraction of the genome using next-generation sequencing (NGS) data, will allow for more in-depth analyses of the dynamics of the genome accompanying polyploidization. This will also allow us to identify major repeat types responsible for creating these levels of observed variation.

The occurrence of B-chromosomes reported in the current study is in agreement with earlier studies of the complex from China [24], Japan [71], and Korea [25]. The occurrence of supernumerary genetic material (either B-chromosomes and/or SCS) has frequently been reported not only for Barnardia, but also for many other species of the family Hyacinthaceae [41, 42, 87, 88]. Although B-chromosomes have often been reported to contribute to significant variation of genome size [45, 89], no positive correlation between the presence of supernumerary genetic materials and genome size variation has been inferred here for Barnardia japonica. This is likely due to rather large absolute genome sizes of Barnardia japonica [45]. Similarly, several other studies have also found no significant impact of B-chromosomes or aneuploidy on genome size variation among individuals with larger genome sizes [90,91,92,93].

The utility of micromorphological characters for taxonomical context

Species delimitation is rather complicated and problematic in some genera of the Hyacinthaceae family [14, 27, 60, 94, 95]. This is particularly true for the genera in which rather uniform morphology contrasts with a striking chromosomal variation, caused by dysploidy, polyploidy as well as presence of supernumerary DNA material [96]. Although micromorphological characters including leaf indumentum, reproductive organs, and pollen grains are often useful for species identification in taxonomically complex groups [56,57,58, 70, 97, 98], no diagnostic structured micromorphological variation has been found in the B. japonica complex with exception of some quantitative morphological characters (e.g., leaf diameter, guard cell size, and pollen grain size). Whereas the leaves of BB diploids tend to be narrower with fewer veins when compared to wider leaves with more veins of cytotype AA, these quantitative characters lose their significance once polyploids are included in the comparison. The variation of vegetative morphological characters has been suggested to be associated with genome size variation/polyploidy [99, 100] or with various environmental factors, most likely the availability of the nutrients and water [101], or combination of both. A correlation between ploidy levels and the size of guard cells have been extensively studied [102,103,104,105] and has been shown to be related to both genome size changes and to environmental stimuli at the same time, thus being very plastic and not allowing direct correlations of stomatal cell sizes and polyploidy [106,107,108,109].

Although we have not analysed populations from the entire range of distribution in the B. japonica complex, which extends to China and Japan [22, 24, 26, 28, 35, 59], our data suggest that the complex in Eastern Asia splits into several distinct evolutionary units on the basis of the available cytological and micromorphological evidence. Further studies involving molecular phylogenomic and cytogenetic analysis of populations representing whole distribution range of the complex, however, are still required to understand the genetic variation as well as the evolutionary history of the complex and test the existing taxonomic treatments.

Odd-ploidy level polyploids (3x, 5x, 7x, 9x, etc.) generally face difficulties in the production of functional gametes owing to irregularities in meiosis and are, therefore, expected to produce high-sterile pollen grains and seeds [110, 111]. On the other hand, triploid may also act as triploid bridge relevant for recurrent origin of higher (and even) ploidy levels [8]. Similarly to other various species complex groups [8], the odd-ploidy cytoytpes occurred in low frequency in the B. japonica complex, and pollen sterility in triploids was higher than that of pentaploids, as reported also in the Senecio carniolicus complex [112]. Although odd-ploidy individuals have been considered an instant postzygotic barrier due to the uneven number of chromosome sets caused by irregularities during meiosis [113, 114], the barrier is often leaky and gene flow may be mediated by at least partially fertile interploidy hybrids [8, 110,111,112, 115].

Pollen micromorphological analyses provided some evidence of correlation of ploidy levels and pollen sizes, with diploid taxa having significantly smaller pollen grains compared to all polyploid cytotypes, as reported also for other plant groups [116,117,118]. These differences were, however, not very pronounced. Although pollen exine ornamentations have been shown to be of great taxonomic significance in Hyacinthaceae [50, 95, 119, 120], the exine sculpturing in Barnardia japonica has shown no variation with reticulate-perforate sculpturing in all analyzed cytotypes. Thus, none of the analyzed leaf and pollen micromorphological characters were of any taxonomic significance and thus cannot be used for identification of the cytotypes in the B. japonica complex. This strongly contrasts with very structured and established karyotypical variation of the group. Clearly, chromosomal changes are one of the driving forces of the evolution of this species complex. Further population genetic and cytogenomic analyses of all ploidy levels from whole distribution area of the species should provide more information about the role of chromosomal changes in the diversification within this complex [121,122,123,124].

Conclusions

The present study constitutes the first step towards a better understanding of the evolutionary history of the Barnardia japonica complex. Although no structured diagnostic micromorphological variation has been found in the B. japonica complex, the reported structured variation of chromosome numbers and genome sizes, provides a baseline for further molecular cytogenetic analyses. Genome skimming NGS data will allow for characterization of complete repeat profiles of all cytotypes and will guide the selection of various repeat types that can be used as chromosomal markers. Further molecular phylogenetic analyses of the FISH-karyotyped individuals will allow testing the existing hypotheses on taxonomic treatments of this species complex and will eventually allow to resolve the taxonomic status of the cytotypes of the B. japonica complex.

Data Availability

All relevant data are within the paper and its Supporting Information files. The data presented in this study will be available on request from the corresponding authors.

References

  1. Guerra M. Chromosome numbers in plant cytotaxonomy: concepts and implications. Cytogenet Genome Res. 2008;120:339–50.

    CAS  PubMed  Google Scholar 

  2. Weiss-Schneeweiss H, Schneeweiss GM. Karyotype diversity and evolutionary trends in angiosperms. In: Leitch IJ, Greilhuber J, Doležel J, Wendel JF, editors. Plant Genome diversity volume 2, physical structure, Behavior and Evolution of Plant Genomes. Vienna: Springer Verlag; 2013. pp. 209–30.

    Google Scholar 

  3. Wendel JF. Genome evolution in polyploids. In: Doyle J, Gaut B, editors. Plant molecular evolution. Dordrecht: Springer; 2000. pp. 225–49.

    Google Scholar 

  4. Otto SP. The evolutionary consequences of polyploidy. Cell. 2007;131:452–62.

    CAS  PubMed  Google Scholar 

  5. Soltis DE, Albert VA, Leebens-Mack J, Bell CD, Paterson AH, Zheng C, et al. Polyploidy and angiosperm diversification. Am J Bot. 2009;96:336–48.

    PubMed  Google Scholar 

  6. Van de Peer Y, Maere S, Meyer A. The evolutionary significance of ancient genome duplications. Nat Rev Genet. 2009;10:725–32.

    PubMed  Google Scholar 

  7. Parisod C, Holderegger R, Brochmann C. Evolutionary consequences of autopolyploidy. New Phytol. 2010;186:5–17.

    CAS  PubMed  Google Scholar 

  8. Kolář F, Čertner M, Suda J, Schönswetter P, Husband BC. Mixed-ploidy species: progress and opportunities in polyploid research. Trends Plant Sci. 2017;22:1041–55.

    PubMed  Google Scholar 

  9. Mairal M, Šurinová M, Castro S, Münzbergová Z. Unmasking cryptic biodiversity in polyploids: origin and diversification of Aster amellus aggregate. Ann Bot. 2018;122:1047–59.

    PubMed  PubMed Central  Google Scholar 

  10. Książczyk T, Kovařík A, Eber F, Huteau V, Khaitova L, Tesarikova Z, et al. Immediate unidirectional epigenetic reprogramming of NORs occurs independently of rDNA rearrangements in synthetic and natural forms of a polyploidy species Brassica napus. Chromosoma. 2011;120:557–71.

    PubMed  Google Scholar 

  11. Weiss-Schneeweiss H, Blöch C, Turner B, Villaseñor JL, Stuessy TF, Schneeweiss GM. The promiscuous and the chaste: frequent allopolyploid speciation and its genomic consequences in american daisies (Melampodium sect. Melampodium; Asteraceae). Evolution. 2012;66:211–28.

    CAS  PubMed  Google Scholar 

  12. McCann J, Jang TS, Macas J, Schneeweiss GM, Matzke NJ, Novak P et al. Dating the species network: allopolyploidy and repetitive DNA evolution in American Daisies (Melampodidum sect. Melampodium, Asteraceae). Syst Biol. 2018;67:1010–24.

  13. Koukalova B, Moraes AP, Renny-Byfield S, Matyasek R, Leitch AR, Kovařík A. Fall and rise of satellite repeats in allopolyploids of Nicotiana over c. 5 million years. New Phytol. 2010;186:148–60.

    CAS  PubMed  Google Scholar 

  14. Jang TS, Parker JS, Emadzade K, Temsch EM, Leitch AR, Weiss-schneeweiss H. Multiple origins and nested cycles of hybridization result in high tetraploid diversity in the monocot Prospero. Front Plant Sci. 2018;9:433.

    PubMed  PubMed Central  Google Scholar 

  15. Malinská H, Jennifer A, Tate JA, Mavrodiev E, Matyášek R, Lim KY, et al. Ribosomal RNA genes evolution in Tragopogon: a story of New and Old World allotetraploids and the synthetic lines. Taxon. 2011;60:348–54.

    Google Scholar 

  16. Gaeta RT, Chris Pires J. Homoeologous recombination in allopolyploids: the polyploidy ratchet. New Phytol. 2010;186:18–28.

    CAS  PubMed  Google Scholar 

  17. Buggs RJA, Renny-Byfield S, Chester M, Jordon-Thaden IE, Facio Viccini LF, Chamala S, et al. Next-generation sequencing and genome evolution in allopolyploids. Am J Bot. 2012;99:372–82.

    PubMed  Google Scholar 

  18. Park I, Choi B, Weiss-Schneeweiss H, So S, Myeong HH, Jang TS. Comparative analyses of complete chloroplast genomes and karyotypes of allotetraploid Iris koreana and its putative diploid parental species (Iris series Chinenses, Iridaceae). Int J Mol Sci. 2022;23:10929.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Stebbins GL. Chromosomal evolution in higher plants. Edward Arnold Ltd; 1971.

  20. Angiosperm Phylogeny Group, Chase MW, Christenhusz MJ, Fay MF, Byng JW, Judd WS, et al. An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG IV. Bot J Linn Soc. 2016;181:1–20.

    Google Scholar 

  21. Chen JS, Gao Q, Zhou H, Huang YS, Ogisu M, Cao M. Distribution, karyomorphology, and morphology of Aspidistra subrotata (Asparagaceae) at different ploidy levels in limestone areas of Asia. Comp Cytogenet. 2017;11:1–13.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Shibata F, Hizume M, Ohashi H, Furukawa S. Multiple origin amphidiploids in the complex species of Scilla scilloides, Asparagaceae. Cytologia. 2017;82:413–22.

    CAS  Google Scholar 

  23. Nath S, Sarkar S, Patil SD, Saha PS, Lekhak MM, Ray S, et al. Cytogenetic diversity in Scilloideae (Asparagaceae): a comprehensive recollection and exploration of karyo-evolutionary trends. Bot Rev. 2022;1–43. https://doi.org/10.1007/s12229-022-09279-1.

  24. Ding K, Ge S, Hong D, Yu Z. Cytotype variation and cytogeography of Scilla sinensis (Louriro) Merrill (Hyacinthaceae) in China. Hereditas. 1998;129:151–60.

    Google Scholar 

  25. Choi HW, Lee WK, Choi EY, Park JH, Bang JW. Geographical distribution of cytotypes and genomic structures in natural populations of the Scilla scilloides complex in Korea. J Plant Biol. 2004;47:322–9.

    Google Scholar 

  26. Shibata F, Hizume M, Ohashi H, Furukawa S. An event preceding genome differentiation in the A genome populations of Scilla scilloides. Cytologia. 2017;82:307–16.

    Google Scholar 

  27. Hamouche Y, Amirouche N, Misset MT, Amisouche R. Cytotaxonomy of autumnal flowering species of Hyacinthaceae from Algeria. Plant Syst Evol. 2010;285:177–87.

    Google Scholar 

  28. Araki H. The distribution of diploids and polyploids of the Scilla scilloides complex in Korea. Genetica. 1985;66:3–10.

    Google Scholar 

  29. Bang JW, Choi HW. Cytogenetic studies of Scilla scilloides complex from Korea III. Karyotype of cytotype BB and B-chromosome composition. Korean J Bot. 1994;36:281–4.

    Google Scholar 

  30. Pfosser M, Speta F. Phylogenetics of Hyacinthaceae based on plastid and DNA sequences. Ann Mo Bot Gard. 1999;86:852–75.

    Google Scholar 

  31. Ali SS, Yu Y, Pfosser M, Wetschnig W. Inferences of biogeographical histories within subfamily Hyacinthoideae using S-DIVA and bayesian binary MCMC analysis implemented in RASP (reconstruct ancestral state in Phylogenies). Ann Bot. 2012;109:95–107.

    PubMed  Google Scholar 

  32. Araki H. Cytogenetics of Scilla scilloides complex IV. Eu- and aneuploid offspring from allotriploids in a natural population. Jpn J Genet. 1972;47:73–83.

    Google Scholar 

  33. Choi HW, Bang JW. Cytogenetic studies of Scilla scilloides complex from Korea I. distribution of genomes and composition and frequencies of B chromosome. Korean J Bot. 1990;33:237–42.

    Google Scholar 

  34. Yu Z, Araki H. The distribution of diploids and polyploids of the Scilla scilloides complex in the northeastern district of China. Bot Mag Tokyo. 1991;104:183–90.

    Google Scholar 

  35. Araki H, Hidaka S, Takahashi S. Cytogenetics of the Scilla scilloides complex VI. The structures of natural populations. Bot Mag Tokyo. 1976;89:83–91.

    Google Scholar 

  36. Choi HW, Kim JS, Lee SH, Bang JW. Physical mapping by FISH and GISH of rDNA loci and discrimination of genomes A and B in Scilla scilloides complex distributed in Korea. Korean J Bot. 2008;51:408–12.

    CAS  Google Scholar 

  37. Hizume M, Araki H. Discrimination of chromosomes belonging to the genomes A and B in polyploids of Scilla scilloides, Liliaceae by genomic in situ hybridization. La Kromosomo. 1996;83:2885–92.

  38. Jones N, Houben A. B chromosomes in plants: escapes from the A chromosome genome? Trends Plant Sci. 2003;8:417–23.

    CAS  PubMed  Google Scholar 

  39. Levin DA, Palestis BG, Jones RN, Trivers R. Phyletic hot spots for B chromosomes in angiosperms. Evolution. 2005;59:962–9.

    PubMed  Google Scholar 

  40. Houben A, Banaei-Moghaddam AM, Klemme S. Biology and evolution of B chromosomes. In: Doyle J, Gaut B, editors. Plant molecular evolution. Dordrecht: Springer; 2013. pp. 149–65.

    Google Scholar 

  41. Jang TS, Parker J, Weiss-Schneeweiss H. Structural polymorphisms and distinct genomic composition suggest recurrent origin and ongoing evolution of B chromosomes in the Prospero autumnale complex (Hyacinthaceae). New Phytol. 2016;210:669–79.

    CAS  PubMed  Google Scholar 

  42. Jang TS, Parker JS, Weiss-Schneeweiss H. Euchromatic supernumerary chromosomal segments – remnants of ongoing karyotype restructuring in the Prospero autumnale complex? Genes. 2018;9:468.

    PubMed  PubMed Central  Google Scholar 

  43. Houben A, Jones N, Martins C, Trifonov V. Evolution, composition and regulation of supernumerary B chromosomes. Genes. 2019;10:161.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Jones RN, Rees H. B chromosomes. New York: Academic; 1982.

    Google Scholar 

  45. Trivers R, Burt A, Palestis BG. B chromosomes and genome size in flowering plants. Genome. 2004;47:1–8.

    CAS  PubMed  Google Scholar 

  46. Jones RN, Viegas W, Houben A. A century of B chromosomes in plants: so what? Ann Bot. 2008;101:767–75.

    PubMed  Google Scholar 

  47. Chase MW, Reveal JL, Fay MF. A subfamilial classification for the expanded asparagalean families Amaryllidaceae, Asparagaceae and Xanthorrhoeaceae. Bot J Linn Soc. 2009;161:132–6.

    Google Scholar 

  48. Lynch AH, Rudall PJ, Cutler DF. Leaf anatomy and systematics of Hyacinthaceae. Kew Bull. 2006;61:145–59.

    Google Scholar 

  49. Martínez-Azorín M, Crespo MB, Juan A. Taxonomic revision of Ornithogalum subg. Ornithogalum (Hyacinthaceae) in the Iberian Peninsula and the Baleric Islands. Plant Syst Evol. 2010;289:181–211.

    Google Scholar 

  50. Heidarian M, Masoumi SM, Amiri S. Palynological study on selected species from Hyacinthaceae with focus on taxonomical implications in Iran. Palynology. 2022;46(1):1–16.

    Google Scholar 

  51. Moon HK, Hong SP, Smets E, Huysmans S. Phylogenetic significance of leaf micromorphology and anatomy in the tribe Mentheae (Nepetoideae: Lamiaceae). Bot J Linn Soc. 2009;160:211–31.

    Google Scholar 

  52. Kong MJ, Hong SP. Leaf micromorphology of the Persicaria sect. Cephalophilon (Polygonaceae) and its systematic re-evaluation. Phytotaxa. 2019;391:167–84.

    Google Scholar 

  53. Do Carmo A, Pacifico R, Gonçales-Silva R, Fidanza K, de Souza LA. Epidermal micromorphology and venation patterns of Microlicieae (Melastomataceae) leaves: looking for new characters for the taxonomy of this neotropical tribe. Flora. 2019;261:151494.

    Google Scholar 

  54. Choi B, Kim SY, Jang TS. Micromorphological and cytological comparisons between Youngia japonica and Youngia longiflora using light and scanning electron microscopy. Microsc Res Tech. 2020;83:1456–63.

    CAS  PubMed  Google Scholar 

  55. Leandro TD, Scatena VL, Clark LG. Comparative leaf blade anatomy and micromorphology in the systematics and phylogeny of Bambusoideae (Poaceae, Poales). Bot J Linn Soc. 2020;192:165–83.

    Google Scholar 

  56. Choi B, Ahn Y, Jang TS. Implications of foliar epidermal micromorphology using light and scanning electron microscopy: a useful tool in taxonomy of Korean irises. Microsc Res Tech. 2022;85:2549–57.

    PubMed  Google Scholar 

  57. Choi B, Ryu J, Jang TS. Can pollen exine ornamentation contribute to species delimitation in Korean Iris L. taxa (Iridaceae)? Palynology. 2022;46(4):1–9.

  58. Choi B, Park I, So S, Myeong HH, Ryu J, Ahn YE, et al. Comparative analysis of two Korean irises (Iris ruthenica and I. uniflora, Iridaceae) based on plastome sequencing and micromorphology. Sci Rep. 2022;12:9424.

  59. Araki H. Cytogenetics of Scilla scilloides complex. III. Homoeology between genomes A (x = 8) and B (x = 9). Jpn J Genet. 1971;46:265–75.

  60. Jang TS, Emadzade K, Parker J, Temsch EM, Leitch AR, Speta F, et al. Chromosomal diversification and karyotype evolution of diploids in the cytologically diverse genus Prospero (Hyacinthaceae). BMC Evol Biol. 2013;13:136.

    PubMed  PubMed Central  Google Scholar 

  61. Choi B, Yang S, Song JH, Jang TS. Karyotype analysis and genome size variation in Ajuga L. (Ajugoideae-Lamiaceae). Nord J Bot. 2019;37:e02337.

    Google Scholar 

  62. Doležel J, Sgorbati S, Lucretti S. Comparison of three DNA fluorochromes for flow-cytometric estimation of nuclear DNA content in plants. Physiol Plant. 1992;85:625–31.

    Google Scholar 

  63. Doležel J, Greilhuber J, Suda J. Estimation of nuclear DNA content in plants using flow cytometry. Nat Protoc. 2007;2:2233–44.

    PubMed  Google Scholar 

  64. Choi B, Weiss-Schneeweiss H, Temsch EM, So S, Myeong MH, Jang TS. Genome size and chromosome number evolution in Korean Iris L. species (Iridaceae Juss). Plants. 2020;9:1284.

  65. Choi B, Gang GH, Kim H, Byun H, Kwak M, So S, et al. Cytological study of Cypripedium japonicum Thunb. (Orchidaceae Juss.): an endangered species from Korea. Plants. 2021;10:1978.

    CAS  PubMed  PubMed Central  Google Scholar 

  66. Pellicer J, Leitch IJ. The application of flow cytometry for estimating genome size and ploidy level in plants. In: Besse P, editor. Molecular Plant Taxonomy. Methods in Molecular Biology. Volume 1115. New Jersey: Humana Press; 2014. pp. 279–307.

    Google Scholar 

  67. Khatun S, Flowers TJ. The estimation of pollen viability in rice. J Exp Bot. 1995;46:151–4.

    CAS  Google Scholar 

  68. Jang TS, Moon HK, Hong SP. Sex expression, population structure, and floral dimorphism in a gynodioecious herb, Agastache rugosa (Lamiaceae) in Korea. Flora. 2015;215:23–32.

    Google Scholar 

  69. Punt W, Hoen PP, Blackmore S, Nilsson S, le Thomas A. Glossary of pollen and spore terminology. Rev Paleobot Palynol. 2007;143:1–81.

    Google Scholar 

  70. Kim H, Choi B, Lee C, McAdam SAM, Paik JH, Jang TS. Micromorphological differentiation of Korean Disporum species using light and scanning electron microscopy. Microsc Res Tech. 2021;84:2614–24.

  71. Haga T, Noda S. Cytogenetics of the Scilla scilloides complex I. karyotype, genome, and population. Genetica. 1976;46:161–76.

    Google Scholar 

  72. Portugal J, Waring MJ. Assignment of DNA binding sites for 4’,6-diamidino-2-phenylindole and bisbenzimide (Hoechst 33258). A comparative footprinting study. Biochim Biophys Acta. 1988;949:158–68.

    CAS  PubMed  Google Scholar 

  73. Temsch EM, Koutecký P, Urfus T, Šmarda P, Doležel J. Reference standards for flow cytometric estimation of absolute nuclear DNA content in plants. Cytometry A. 2021;101:710–24.

    PubMed  PubMed Central  Google Scholar 

  74. Doležel J, Greilhuber J, Lucretti S, Meister A, Lysák MA, Nardi L, et al. Plant genome size estimation by flow cytometry: inter-laboratory comparison. Ann Bot. 1998;82:17–26.

    Google Scholar 

  75. Noirot M, Barre P, Louarn J, Duperray C, Hamon S. Consequences of stoichiometric error on nuclear DNA content evaluation in Coffea liberica var. dewevrei using DAPI and propidium iodide. Ann Bot. 2002;89:385–9.

  76. Doležel J, Bartos JAN. Plant DNA flow cytometry and estimation of nuclear genome size. Ann Bot. 2005;95:99–110.

    PubMed  PubMed Central  Google Scholar 

  77. Ding K, Ge S, Yu Z, Hong D. Allozyme variation in the diploid (A genome) populations of Scilla scilloides (Hyacinthaceae). Plant Syst Evol. 1999;218:23–31.

  78. Renny-Byfield S, Chester M, Kovařík A, Le Comber SC, Grandbastien MA, Deloger M, et al. Next generation sequencing reveals genome downsizing in allotetraploid Nicotiana tabacum, predominantly through the elimination of paternally derived repetitive DNAs. Mol Biol Evol. 2011;28:2843–54.

    CAS  PubMed  Google Scholar 

  79. Simonin KA, Rody AB. Genome downsizing, physiological novelty, and the global dominance of flowering plants. PLoS Biol. 2018;16:e2003706.

    PubMed  PubMed Central  Google Scholar 

  80. Wang X, Morton JA, Pellicer J, Leitch IJ, Leitch AR. Genome downsizing after polyploidy: mechanisms, rates and selection pressures. Plant J. 2021;107:1003–15.

    CAS  PubMed  Google Scholar 

  81. Leitch IJ, Bennett MD. Genome downsizing in polyploid plants. Biol J Linn Soc. 2004;82:651–63.

    Google Scholar 

  82. Weiss-Schneeweiss H, Emadzade K, Jang TS, Schneeweiss GM. Evolutionary consequences, constraints and potential of polyploidy in plants. Cytogenet Genome Res. 2013;140:137–50.

    CAS  PubMed  Google Scholar 

  83. Báez M, Vaio M, Dreissig S, Schubert V, Houben A, Pedrosa-Harand A. Together but different: the subgenomes of the bimodal Eleutherine karyotypes are differentially organized. Front Plant Sci. 2019;10:1170.

    PubMed  PubMed Central  Google Scholar 

  84. Dodsworth S, Chase MW, Leitch AR. Is post-polyploidization diploidization the key to the evolutionary success of angiosperms? Bot J Linn Soc. 2016;180:1–5.

    Google Scholar 

  85. Mandáková T, Gloss AD, Whiteman NK, Lysak MA. How diploidization turned a tetraploid into a pseudotriploid. Am J Bot. 2016;103:1187–96.

    PubMed  Google Scholar 

  86. Soltis DE, Visger CJ, Marchant DB, Soltis PS. Polyploidy: pitfalls and paths to a paradigm. Am J Bot. 2016;103:1146–66.

    PubMed  Google Scholar 

  87. Garrido-Ramos MA, Jamilena M, de la Herrán R, Ruiz Rejón C, Camacho JPM, Ruiz-Rejón M. Inheritance and fitness effects of a pericentric inversion and a supernumerary chromosome segment in Muscari comosum (Liliaceae). Heredity. 1998;80:724–31.

    Google Scholar 

  88. Weiss-Schneeweiss H, Riha K, Jang CG, Puizina J, Scherthan H, Schweizer D. Chromosome termini of the monocot plant Othocallis siberica are maintained by telomerase, which specifically synthesizes vertebrate-type telomere sequences. Plant J. 2004;37:484–93.

    CAS  PubMed  Google Scholar 

  89. Rosado TB, Clarindo WR, Carvalho CR. An integrated cytogenetic, flow and image cytometry procedure used to measure the DNA content of Zea mays A and B chromosomes. Plant Sci. 2009;176:154–8.

    CAS  Google Scholar 

  90. Poggio L, Rosato M, Chiavarino AM, Naranjo CA. Genome size and environmental correlation in maize (Zea mays ssp. mays, Poaceae). Ann Bot. 1998;82:107–15.

    Google Scholar 

  91. Leong-Škorničková J, Šída O, Jarolímová V, Sabu M, Fér T, Trávníček P, et al. Chromosome numbers and genome size variation in indian species of Curcuma (Zingiberaceae). Ann Bot. 2007;100:505–26.

    PubMed  PubMed Central  Google Scholar 

  92. Leitch IJ, Kahandawala I, Suda J, Hanson L, Ingrouille MJ, Chase MW, et al. Genome diversity in orchids: consequences and evolution. Ann Bot. 2009;104:469–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  93. Chumová Z, Mandáková T, Trávníček P. Are B-chromosomes responsible for the extraordinary genome size variation in selected Anthoxanthum annuals? Plant Syst Evol. 2016;302:731–8.

    Google Scholar 

  94. Moret J. Numerical taxonomy applied to a study of some ploidy levels within the Ornithogalum umbellatum complex (Hyacinthaceae) in France. Nord J Bot. 1992;12:183–95.

    Google Scholar 

  95. Yadav PB, Manning JC, Lekhak MM. A cytotaxonomic revision of Drimia Jacq. (Hyacinthaceae: Urgineoideae) in India. S Afr J Bot. 2019;123:62–86.

    Google Scholar 

  96. Jang TS. Chromosomal evolution in Prospero autumnale complex. PhD thesis, University of Vienna, Vienna, Austria; 2013.

  97. Leandro TD, Scatena VL, Clark LG. The contribution of foliar micromorphology and anatomy to the circumscription of species within the Chusquea ramosissima informal group (Poaceae, Bambusoideae, Bambuseae). Plant Syst Evol. 2017;303:745–56.

    Google Scholar 

  98. Ketjarun K, Traiperm P, Suddee S, Watthana S, Gale SW. Labellar anatomy of the Nervilia plicata complex (Orchidaceae: Epidendroideae) in tropical Asia. Kew Bull. 2019;74:1.

    Google Scholar 

  99. Zonneveld BJM, Duncan GD. Genome size for the species of Nerine Herb. (Amaryllidaceae) and its evident correlation with growth cycle, leaf width and other morphological characters. Plant Syst Evol. 2006;257:251–60.

    CAS  Google Scholar 

  100. Hanušová K, Ekrt L, Vít P, Kolář F, Urfus T. Continuous morphological variation correlated with genome size indicates frequent introgressive hybridization among Diphasiastrum species (Lycopodiaceae) in Central Europe. PLoS ONE. 2014;9:e99552.

    PubMed  PubMed Central  Google Scholar 

  101. Hodálová I, Mártonfiová L, Skokanová K, Majerová M, Somlyay L, Mereďa P Jr. The utility of genome size in plant identification: a case study on Sesleria (Poaceae) from Croatia and Slovenia. Plant Syst Evol. 2020;306:87.

    Google Scholar 

  102. Beaulieu JM, Leitch IJ, Patel S, Pendharkar A, Knight CA. Genome size is a strong predictor of cell size and stomatal density in angiosperms. New Phytol. 2008;179:975–86.

    PubMed  Google Scholar 

  103. Hodgson JG, Sharafi M, Jalili A, Díaz S, Montserrat-Martí G, Palmer C, et al. Stomatal vs. genome size in angiosperms: the somatic tail wagging the genomic dog? Ann Bot. 2010;105:573–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  104. Veselý P, Bureš P, Šmarda P, Pavlíček T. Genome size and DNA base composition of geophytes: the mirror of phenology and ecology? Ann Bot. 2012;109:65–75.

    PubMed  Google Scholar 

  105. Veselý P, Šmarda P, Bureš P, Stirton C, Muasya AM, Mucina L, et al. Environmental pressures on stomatal size may drive plant genome size evolution: evidence from a natural experiment with Cape geophytes. Ann Bot. 2020;126:323–30.

    PubMed  PubMed Central  Google Scholar 

  106. Franks PJ, Drake PL, Beerling DJ. Plasticity in maximum stomatal conductance constrained by negative correlation between stomatal size and density: an analysis using Eucalyptus globulus. Plant Cell Environ. 2009;32:1737–48.

    PubMed  Google Scholar 

  107. Doheny-Adams T, Hunt L, Franks PJ, Beerling DJ, Gray JE. Genetic manipulation of stomatal density influences stomatal size, plant growth and tolerance to restricted water supply across a growth carbon dioxide gradient. Philos Trans R Soc Lond B Biol Sci. 2012;367:547–55.

    CAS  PubMed  PubMed Central  Google Scholar 

  108. Lomax BH, Hilton J, Bateman RM, Upchurch GR, Lake JA, Leitch IJ, et al. Reconstructing relative genome size of vascular plants through geological time. New Phytol. 2014;201:636–44.

    CAS  PubMed  Google Scholar 

  109. Sun Y, Yan F, Cui X, Liu F. Plasticity in stomatal size and density of potato leaves under different irrigation and phosphorus regimes. J Plant Physiol. 2014;171:1248–55.

    CAS  PubMed  Google Scholar 

  110. Ramsey J, Schemske DW. Pathways, mechanisms, and rates of polyploid formation in flowering plants. Annu Rev Ecol Syst. 1998;29:467–501.

    Google Scholar 

  111. Morgan EJ, Certner M, Lucanova M, Deniz U, Kubikova K, Venon A, et al. Disentangling the components of triploid block and its fitness consequences in natural diploid-tetraploid contact zones of Arabidopsis arenosa. New Phytol. 2021;232:1449–62.

    CAS  PubMed  Google Scholar 

  112. Peskoller A, Silbernagl L, Hulber K, Sonnleitner M, Schonswetter P. Do pentaploid hybrids mediate gene flow between tetraploid Senecio disjunctus and hexaploid S. carniolicus s. str. (S. carniolicus aggregate, Asteraceae)? Alpine Bot. 2021;131:151 – 60.

  113. Felber F. Establishment of a tetraploid cytotype in a diploid population: effect of relative fitness of the cytotypes. J Evol Biol. 1991;4:195–207.

    Google Scholar 

  114. Köhler C, Scheid OM, Erilova A. The impact of the triploid block on the origin and evolution of polyploid plants. Trends Genet. 2010;26:142–8.

    PubMed  Google Scholar 

  115. Husband BC. Constraints on polyploid evolution: a test of the minority cytotype exclusion principle. Proc R Soc Lond B Biol Sci. 2000;267(1440):217–23.

    CAS  Google Scholar 

  116. Altmann T, Damm B, Frommer WB, Martin T, Morris PC, Schweizer D, et al. Easy determination of ploidy level in Arabidopsis thaliana plants by means of pollen size measurement. Plant Cell Rep. 1994;13:652–6.

    CAS  PubMed  Google Scholar 

  117. De Storme N, Zamariola L, Mau M, Sharbel TF, Geelen D. Volume-based pollen size analysis: an advanced method to assess somatic and gametophytic ploidy in flowering plants. Plant Reprod. 2013;26:65–81.

    CAS  PubMed  Google Scholar 

  118. Srisuwan S, Sihachakr D, Martín J, Vallès J, Ressayre A, Brown SC, et al. Change in nuclear DNA content and pollen size with polyploidisation in the sweet potato (Ipomoea batatas, Convolvulaceae) complex. Pl Biol (Stuttg). 2019;21:237–47.

    CAS  Google Scholar 

  119. Ghavami E, Jamzad Z, Tavasoli A. Evaluation of pollen morphology as a taxonomic character for generic delimitation in Scilla s.l. (Hyacinthaceae). Nord J Bot. 2009;27:510–5.

    Google Scholar 

  120. Citak BY, Uysal T, Bozkurt M, Demirelma H, Aksoy A, Ertuğrul K. The comparative palynomorphological studies on turkish Muscari genera and close relatives (Asparagaceae-Scilloideae) with their taxonomic implications. Microsc Res Tech. 2022;85:2692–707.

    PubMed  Google Scholar 

  121. Rivera-Lugo M, García-Mendoza A, Simpson J, Solano E, Gil-Vega K. Taxonomic implications of the morphological and genetic variation of cultivated and domesticated populations of the Agave angustifolia complex (Agavoideae, Asparagaceae) in Oaxaca, Mexico. Plant Syst Evol. 2018;304:969–79.

    Google Scholar 

  122. Caković D, Cresti L, Stešević D, Schönswetter P, Frajman B. High genetic and morphological diversification of the Euphorbia verrucosa alliance (Euphorbiaceae) in the Balkan and Iberian peninsulas. Taxon. 2021;70:286–307.

    Google Scholar 

  123. López-González N, Bobo-Pinilla J, Padilla-García N, Loureiro J, Castro S, Rojas-Andrés BM, et al. Genetic similarities versus morphological resemblance: unraveling a polyploid complex in a Mediterranean biodiversity hotspot. Mol Phylogenet Evol. 2021;155:107006.

    PubMed  Google Scholar 

  124. Suissa JS, Kinosian SP, Schafran PW, Bolin JF, Taylor WC, Zimmer EA. Homoploid hybrids, allopolyploids, and high ploidy levels characterize the evolutionary history of a western north american quillwort (Isoëtes) complex. Mol Phylogenet Evol. 2022;166:107332.

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Soonku So (Korea National Park Research Institute) and CNUK lab members for their help with field work. The authors thank Dr. Sangho Choi (Korea Research Institute of Bioscience and Biotechnology) to allow using the scanning electron microscope (SEM). The authors would like to thank Dr. Eva M. Temsch (University of Vienna, Austria) for the kind provision of seeds for reference standard for the genome size measurements.

Funding

Open access funding provided by University of Vienna. This work was financially supported by KR03/2020 (Austria) to H. W.-S. and the National Research Foundation of Korea (NRF) funded by the Korea government (NRF-2019K1A3A1A1811607011) to T.-S. Jang.

Open access funding provided by University of Vienna.

Author information

Authors and Affiliations

Authors

Contributions

H.W.-S. and T.-S.J. conceived and designed the study. B.C., H.W.-S., and T.-S.J. wrote the manuscript. H.K., B.C., and T.-S.J. performed the cytological and micromorphological work. B.C., H.K., H.W.-S., and T.-S.J. analyzed the data. H.K., B.C., C.L., J.-H.P., C.-G.J., and T.-S.J. collected data and samples. All authors contributed to the data interpretation and approved the final version of the manuscript.

Corresponding authors

Correspondence to Hanna Weiss-Schneeweiss or Tae-Soo Jang.

Ethics declarations

Ethics approval and consent to participate

All methods were performed in accordance with the relevant guidelines and regulations of institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

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.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, H., Choi, B., Lee, C. et al. Does the evolution of micromorphology accompany chromosomal changes on dysploid and polyploid levels in the Barnardia japonica complex (Hyacinthaceae)?. BMC Plant Biol 23, 485 (2023). https://doi.org/10.1186/s12870-023-04456-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04456-9

Keywords