Genome and chromosome size variation in the study group Poodae
2C values. Representatives of the supertribe Poodae, including the tribes Aveneae, Festuceae, Poeae and the intertribe hybrids, had 2C values, i.e. holoploid diplophasic or sporophytic genome sizes of the non-replicated nuclei, ranging between 1.49 pg in diploid Poa supina and 44.75 pg in 16–18x Anthoxanthum amarum, respectively. All of the 214 species and 399 accessions sampled in our study fell between these minimum and maximum values, which have been found in previous studies also using FCM + PI (Mao and Huff 2012; Chumová et al. 2015). Following the genome size categories of Leitch et al. (1998), 2.3% of our studied Poodae species had “very small” genomes of ≤ 2.8 pg/2C, 31.3% had “small genomes” (≤ 7 pg/2C), while the vast majority of 49% had “medium-sized” (> 7 and < 14 pg/2C), 16.4% “large” (≥ 14 and < 35 pg/2C) and 1.4% “very large genomes” of ≥ 35 pg/2C (Tables 2, 3; Figs. 1, 2A; Online Resource 1). Compared to the rest of the subfamily Pooideae, the 2C values of the Poodae were mostly lower than those of their sister supertribe Triticodae but higher than those of the phylogenetically ‘early-diverging’ lineages of the subfamily Pooideae (see Winterfeld et al. 2024).
1Cx values. The genome size of the monoploid non-replicated Poodae chromosome sets found in this study ranged from 1.13 pg/1Cx in one of the Holcus mollis accessions to 8.30 pg/1Cx in Echinaria capitata (Tables 2, 3; Figs. 1, 2B; Online Resource 1). However, these values were exceeded by a minimum value of 0.75 pg/1Cx found in Poa supina (Mao and Huff 2012) and the maximum value of 9.19 pg/1Cx found in Anthoxanthum gracile (Chumová et al. 2015). Compared to the supertribe Triticodae, the sister lineage of the Poodae, the values were therefore not significantly different from those of the tribe Bromeae but mostly smaller than those of the Triticeae (Winterfeld et al. 2024).
Mean chromosome DNA content (MC). The calculated mean chromosomes sizes of the studied Poodae ranged from 0.16 pg in the same Holcus mollis accession to 1.16 pg in Tricholemma jahandiezii (Tables 2, 3; Figs. 1, 2C; Online Resource 1), but the values were surpassed by the minimum MC of 0.11 pg previously recorded for Poa supina and the maximum value of 1.84 pg for Anthoxanthum gracile (Mao and Huff 2012; Chumová et al. 2015). The majority of the values were between 0.2 and 0.5 pg. As in the case of the 1Cx values, the MCs of the Poodae did not differ significantly from those of the tribe Bromeae in the sister supertribe Triticodae, but clearly from those of the Triticeae with mostly significantly larger chromosome sizes (Winterfeld et al. 2024).
Tribe Aveneae and its subtribes
The 2C values in the tribe Aveneae varied between 2.48 pg (MV) in diploid Torreya pallida and 33.69 pg in 14x Helictotrichon filifolium subsp. filifolium (Tables 2, 3; Figs. 1, 2A; Online Resource 1) with a maximum value of this tribe of 44.75 pg mentioned above for 16–18x Anthoxanthum amarum, implying an approximately 18-fold variation. The majority of holoploid genome sizes was between 5 pg and 10 pg.
The monoploid genome sizes were mostly between 1.22 pg/1Cx found for Amphibromus nervosus and 8.1 pg/1Cx in Tricholemma jahandiezii, the latter value being exceede only by 9.19 pg in Anthoxanthum gracile (see above). The majority of monoploid genome sizes was 2–3 pg/1Cx (Tables 2, 3; Figs. 1, 2B; Online Resource 1).
With respect to chromosome sizes, most MCs were 0.2–0.5 pg, with the outliers of exceptionally large chromosome sizes of > 1.1 pg again in T. jahandiezii and A. gracile (Tables 2, 3; Figs. 1, 2C; Online Resource 1).
Considering only the diploid Aveneae taxa, the maximum 1Cx value of species with x = 7 was 4.98 pg/1Cx in Avena hispanica (2n = 14) (Table 2; Online Resource 1), while the diploids with x = 5 had up to 9.19 pg in Anthoxanthum gracile (2n = 10) (Chumová et al. 2015).
The tribe Aveneae comprises about 54 genera and 1022 species and has a cosmopolitan distribution with a focus outside of the tropics and subtropics, namely in the temperate climatic zones of the Earth. Representatives of eight subtribes of the Aveneae were examined, arranged in alphabetical order (Tables 1–3; Figs. 2, 3; Online Resource 1).
Agrostidinae. The genera Agrostis s.l. (including Chaetopogon, Polypogon), Calamagrostis (including Ammophila, Deyeuxia) and Gastridium were studied for this predominantly northern hemisphere subtribe, comprising a total of 22 examined species and 31 accessions. The 2C values ranged from 3.83 pg for the diploid Agrostis mediterranea (syn. Polypogon maritimus) to 24.06 pg for the polyploid Calamagrostis scabrescens. For this species, only a tetraploid chromosome number reported (see CCDB), whereas judging by the 2C value, our specimen must have a higher ploidy of probably 10–12x. The 1Cx values of Agrostidinae ranged from 1.68 pg in Agrostis schraderiana to 2.83 pg (MV) in Gastridium ventricosum, and the MC from 0.24 pg to 0.40 pg in the same species. The highest 1Cx values and MCs of the whole subtribe were found in G. phleoides and G. ventricosum with 2.72–2.82 pg/1Cx and an MC of about 0.40 pg.
Among the Agrostis species studied, A. gigantea and A. subspicata (syn. Chaetopogon fasciculatus) had the highest 1Cx values of 2.35 pg and 2.38 pg, respectively (MVs). The two studied species of the former Polypogon, a genus included in Agrostis (Röser and Tkach 2024), the diploid A. mediterranea and the tetraploid A. alopecuroides largely agreed with the majority of Agrostis species with respect to their 1Cx values (1.87–1.92 pg) and MCs (both 0.27 pg). Agrostis species have often been placed in different sections (reviewed by Saarela et al. 2017; Peterson et al. 2020), mainly based on the presence/absence and length of the palea. However, the genome size data of representatives of both groups do not differ significantly, e.g. A. rupestris (short palea) with a 1Cx value of 1.77 pg and a MC of 0.25 pg compared to A. capillaris and A. castellana (long palea) with 1.71–1.75 pg/1Cx and MCs of 0.24–0.25 pg.
Most Calamagrostis species had 1Cx values of 1.85–2.27 pg/1Cx and MCs of 0.26–0.32 pg. Thus, the overall variation of 1Cx values and chromosome sizes (MC) was moderate and the representatives of the different lineages identified in Calamagrostis by molecular phylogenetics (Peterson et al. 2022), such as C. arundinacea (Deyeuxia group, clade B1), C. epigejos and C. pseudophragmites (Epigejos group, clade B2), C. scabrescens (Orientalis group, clade C1), C. purpurea (Purpurea group, clade C2), C. canescens and C. villosa (Calamagrostis group, clade C3), did not show clear differences, considering that chromosome numbers are not known or not reliably known for all taxa.
Calamagrostis rivalis, endemic to Saxony, Germany, has been reported to be octoploid (Heine 1970, 1972; Schiebold et al. 2009 as C. pseudopurpurea), but the accession we studied (11.45 pg/2C) probably had a lower ploidy level and therefore 1Cx and MC values that are probably consistent with those of the other Calamagrostis species. The same is true for the C. scabrescens accession studied (24.06 pg/2C). This species has been reported as tetraploid (CCDB), but probably has a higher ploidy.
Previous studies using FCM + PI in Agrostis capillaris, A. gigantea, A. rupestris, Calamagrostis arundinacea, C. canescens and C. epigeios in accessions with probably the same ploidy level as our samples (Šmarda et al. 2013, 2019) found genome sizes comparable in magnitude to this study, although these values were consistently lower (7–16%), except for A. rupestris, which differed by as much as 27%, a case that needs to be reexamined. The genome size of 5.42 pg/2C for C. arundinacea, also determined by FCM + PI (Pustahija et al. 2013), was about 20% lower than our estimate, although both accessions were presumably tetraploid. However, the values recorded for A. capillaris and C. canescens (Zonneveld 2019) differ from our estimates by only 0.1–0.35%, while A. gigantea was reported to have 13.2 pg/2C instead of 9.40 pg (MV) in our study, which is 1.4 times higher, indicating that it was obtained in an accession belonging to the 6x cytotype with 2n = 42, whereas in this study the 4x cytotype with 2n = 28 was used, as verified by chromosome counting for one of our three accessions studied.
However, the 2C values of 1.1 pg and 1.5 pg reported by Bai et al. (2012) for A. gigantea, which were also determined by FCM + PI appear to be incorrect and are likely based on misidentification of the plants studied.
Anthoxanthinae. This subtribe consists only of the genus Anthoxanthum s.l. if Hierochloe is included, and its species have either x = 7 (traditional Hierochloe) and x = 5 (Anthoxanthum s.s.) as chromosome base numbers. The 2C values of the studied species ranged from 7.83 pg/2C (MV) in diploids to 25.44 pg/2C in polyploids. However, the whole range in Anthoxanthum s.l. is even higher, ranging from 5.52 pg/2C in diploid A. alpinum to 44.75 pg/2C in 16x–18x A. amarum (Chumová et al. 2015), both with x = 5.
For the x = 7 taxa, the 1Cx values found were 4.73 pg in the diploid Central European A. australe (MV) and 4.88 pg in the presumably tetraploid A. nitens, a widespread Holarctic species. The 2C value of 25.44 pg for A. monticola from Alaska suggests that this accession was hexaploid (6x) rather than octoploid (8x) as usually found in this species (see CCDB). The 1Cx value of A. monticola would thus be 4.24 pg, which fits the other x = 7 species. The MC of all these x = 7 species was 0.61–0.70 pg.
Genome sizes of Anthoxanthum species from the southern hemisphere were estimated by FCM + PI in five polyploid species (4x–12x, 2n = 28–84) from New Zealand (as Hierochloe) (Murray et al. 2005). Their 2C values ranged from 18.10 pg to 29.97 pg, the (here calculated) 1Cx values from 2.32 pg to 3.14 pg and the MCs from 0.33 pg to 0.45 pg. Monoploid genome sizes (1Cx values) and chromosomes (MC) are thus significantly smaller than in the northern hemisphere x = 7 taxa. This is most likely real and is not an artifact based on methodological issues, since the DNA C-values reported by Murray et al. (2005) always agree very well with ours when the same species have been studied in both, e.g. in Deschampsia, Koeleria (Trisetum), Pentapogon, (Dichelachne) and Poa (see below), and also in the grasses of the tribe Stipeae (see Winterfeld et al. 2024).
The perennial tetraploid Anthoxanthum odoratum, belonging to the x = 5 taxa, had a 1Cx value of 3.27 pg (MV). The MC was 0.65 pg (MV) and therefore not significantly different from that of the x = 7 taxa. For A. odoratum, 1Cx values of 3.22 pg and 3.27 pg were previously found, which were also determined by FCM (Chumová et al. 2015; Zonneveld 2019). These values agree with our results, whereas the values of 2.83 pg and 2.89 pg found by Šmarda et al. (2013, 2019) are somewhat lower. The highly polyploid perennial A. amarum (16x–18x) with a 2C value of 44.75 pg (Chumová et al. 2015) had a 1Cx value of 2.48–2.80 pg and an MC of 0.50–56 pg, which is lower than in tetraploids and diploids.
The 1Cx values of diploids based on x = 5 in the perennials A. alpinum and A. maderense averaged 2.76 pg and 3.48 pg, respectively (Chumová et al. 2015). Anthoxanthum alpinum has therefore the smallest monoploid genome found in this genus. The diploid annual A. aristatum had 3.92 pg/1Cx (MV) (Table 2), which corresponds to the previously reported 3.83 pg (Chumová et al. 2015), while the diploid annual A. gracile had 9.19 pg/1Cx (Chumová et al. 2015), the highest 1Cx value found so far in the genus.
Aveninae. A total of 15 genera, 29 species and 57 accessions of the mainly northern hemisphere subtribe Aveninae were studied. The 2C values ranged from 3.66 pg in the diploid Gaudinia fragilis (MV) to 33.69 pg in the dodecaploid (12x) Helictotrichon filifolium subsp. filifolium when both diploids and polyploids are included. The 1Cx values of the diploids ranged from 1.83 pg (MV) in G. fragilis to 4.98 pg in Avena hispanica, both annuals. Their MCs ranged from 0.26 pg to 0.71 pg, respectively. In marked contrast, the 1Cx values of the polyploid species, considering those with x = 7, ranged from 2.14 pg in the verified octoploid Acrospelion distichophyllum (syn. Trisetum distichophyllum) to 8.10 pg in the tetraploid Tricholemma jahandiezii, a narrowly distributed endemic of the Moroccan Middle Atlas with an exceptionally large monoploid genome. Both are perennials.
Comparatively high 1Cx values (4.14–5.48 pg) were found in Arrhenatherum elatius (MV), Trisetopsis elongata, Sibirotrisetum sibiricum (syn. Trisetum sibiricum) and the aforementioned Avena hispanica and A. macrostachya, all perennials except for A. hispanica. Their MCs were 0.59–0.78 pg.
In the majority of the Aveninae studied (genera Graphephorum, Helictotrichon, Koeleria, Lagurus, Limnodea, Rostraria, Sphenopholis) the 1Cx values ranged from 2.15 pg/1Cx in Koeleria litvinowii to 3.52 pg/1Cx in Helictotrichon sempervirens.
Within the genera Helictotrichon and Koeleria, several species were examined. In Helictotrichon, six species and seven accessions covering the diploid to dodecaploid level, had 2C values ranging from 5.64 pg (MV) in H. sedenense subsp. sedenense to 33.69 pg in H. filifolium subsp. filifolium. The 1Cx values ranged from 2.82 pg/1Cx (MV) in H. sedenense to 3.52 pg/1Cx in H. sempervirens. The MCs were 0.40–0.50 pg. The chromosome numbers of H. aff. tianschanicum and H. tibeticum are not known, but judging from the 2C values, these species are probably tetra- and hexaploid, respectively.
A total of six species and ten accessions of the genus Koeleria were sampled, covering the ploidy levels of 2x–10x and including also K. spicata (syn. Trisetum spicatum). Chromosome numbers were counted for two accessions, taken from CCDB in the other cases, and inferred from the 2C value in the case of 3x K. vallesiana. The 2C values ranged from 5.34 pg/2C in diploid K. glauca to 22.03 pg/2C (MV) in decaploid K. pyramidata. The 1Cx values were 2.67 pg/1Cx in K. glauca and slightly lower, 2.15–2.58 pg/1Cx, in the polyploids, and the MCs were relatively uniform in all species (0.31–0.38 pg).
Our sample of the Aveninae included two species with a dysploid chromosome number, Trisetum flavescens (2n = 24) and Rostraria cristata (2n = 26) with 5.79 pg/2C and 7.43 pg/2C, respectively. Assuming that they were hypotetraploids, the 1Cx values would be about 1.45 pg and 1.86 pg, respectively, which is significantly lower than in all other Aveninae polyploids. The MCs were 0.24 pg and 0.29 pg.
The previously reported 2C values of 5.20–22.89 pg for 2x–10x taxa of Koeleria, also estimated using FCM + PI, agree with our results, as do the 1Cx values of (recalculated) 2.29–2.60 pg/1Cx (Pečinka et al. 2006). The 2C values of 5.01–5.42 pg and 8.72 pg (Zonneveld 2019) most likely refer to diploid and tetraploid Koeleria accessions, respectively. 4.06 pg/2C was recorded for diploid K. glauca, 4.43 pg/2C for a diploid K. macrantha accession, and 19.94 pg/2C for 10x K. pyramidata (Šmarda et al. 2019), which are values slightly lower than ours, as mentioned above. The 2C value of 9.92 pg estimated for tetraploid K. spicata (as Trisetum spicatum) from New Zealand (Murray et al. 2005) agrees with our result of 10.08 pg for an accession from the Tibetian Plateau.
The comparatively lew 2C values determined by FCM for other Aveninae taxa are 14.74 pg and 15.15 pg in Arrhenatherum elatius (Šmarda et al. 2013, 2019), 8.80 pg in Avena hispanica and 21.78 pg in A. macrostachya (both Yan et al. 2016), 3.44 pg in Gaudinia fragilis, 6.68 pg in Lagurus ovatus, 7.49 pg in R. cristata and 5.53 pg in Trisetum flavescens (Zonneveld 2019), which is in good agreement with our estimates, while the 4.71 pg/2C recorded for T. flavescens (Šmarda et al. 2013) is slightly lower.
Brizinae. The studied representatives of the widespread Eurasian subtribe Brizinae, excluding Macrobriza maxima (see below), had 2C values of 5.73 pg (MV) in the annual Mediterranean Briza minor (2n = 2x = 10) and 6.59 pg/2C and 13.55 pg/2C (MV), respectively, in the widespread Eurasian perennial B. media with x = 7. The former value refers to a diploid accession of B. media with 2n = 14 verified by chromosome counting from the Swiss Alps, the latter to tetraploids (2n = 4x = 28) of unknown origin. The 1Cx values were 2.87 pg (MV) in B. minor, 3.29 pg and 3.39 pg (MV) in the two cytotypes of B. media. The MCs were 0.57 pg (MV) in B. minor and 0.47–0.48 pg (MV) in B. media. Briza minor with x = 5 thus has an approximately 0.9-fold smaller monoploid genome than B. media with x = 7, but approximately 1.2-fold larger chromosomes. Previously reported 2C values were 6.68 pg, 11.92 pg and 12.38 pg for B. media (Šmarda et al. 2013, 2019; Zonneveld 2019), which obviously also refer to diploids and tetraploids, respectively. The value of 0.58 pg for B. minor (Siljak-Yakovlev et al. 2019) is most likely a calculation or typing error, as 5.8 pg would correspond to our result.
Calothecinae. Chascolytrum subaristatum with the verified chromosome number of 2n = 4x = 28 from this mainly South American subtribe had a 2C value of 8.94 pg, a 1Cx value of 2.24 pg and an MC of 0.32 pg.
Echinopogoninae. This Australasian and South American subtribe was represented by three Pentapogon species. Pentapogon crinitus and P. micranthus, both previously mostly assigned to their own genus Dichelachne, were decaploid with 2n = 10x = about 70, verified by chromosome counting, and had similar 2C values of 16.14 pg and 15.83 pg and 1Cx values of 1.61 pg and 1.58 pg, respectively, and an MC of 0.23 pg each. The octoploid Pentapogon quadrifidus var. quadrifidus (2n = 8x = about 56) had a 2C value of 14.91 pg, a 1Cx value of 1.86 pg and an MC of 0.27 pg. These values are consistent with previous genome size reports of 16.36 pg/2C in the decaploid P. crinitus (as D. crinita) and 16.85 pg/2C in the also decaploid P. micranthus (as D. micrantha) (Murray et al. 2005).
Phalaridinae. Within the monogeneric, almost cosmopolitan subtribe Phalaridinae, the three species studied had 2C values between 9.43 pg in the diploid annual P. canariensis, and 10.50 pg in the tetraploid perennial P. aquatica (MVs each). Phalaris canariensis is characterized by the dysploid monoploid chromosome number of x = 6, its 1Cx was 4.72 pg and the MC was 0.79 pg (MVs each). Phalaris aquatica and P. arundinacea, both with x = 7, had significantly lower 1Cx values of 2.63 pg (MV) and 2.46 pg and MCs of 0.38 pg (MV) and 0.35 pg, respectively. Our results are consistent with previous estimates obtained using FCM + PI, i.e. 9.61 pg/2C for P. canariensis (Zonneveld 2019), and 9.80–10.50 pg, 9.08 pg and 10.50 pg/2C for P. arundinacea (Bai et al. 2012; Šmarda et al. 2019; Zonneveld 2019).
Torreyochloinae. Both genera of this subtribe were studied, the disjunct Australasian–South American Amphibromus and the North American–East Asian Torreyochloa. The hexaploid Australian A. nervosus had a 2C value of 7.30 pg, a 1Cx value of 1.22 pg and an MC of 0.17 pg (MVs each), while the diploid American T. pallida had 2.48 pg/2C, 1.24 pg/1Cx and an MC of 0.18 pg (MVs each). The monoploid genome sizes and mean genome sizes of both taxa are thus quite similar.
Macrobriza. This genus probably originated as a hybrid taxon between a subtribe Brizinae (maternal line) and an ancestor of the subtribe Aveninae (Tkach et al. 2020). The only species is the annual Mediterranean M. maxima (syn. Briza maxima). The mean 2C value of diploid M. maxima (2n = 14) was 9.32 pg, in good agreement with a previous genome size estimate of 9.0 pg (Zonneveld 2019). The 1Cx value of M. maxima was 4.66 pg and the MC was 0.67 pg. Thus, the monoploid genome of M. maxima is approximately 1.4–1.6 times larger and the chromosomes approximately 1.2–1.4 times larger than in B. minor and B. media.
Tribe Festuceae and its subtribes
The holoploid 2C genome sizes in the studied Festuceae taxa ranged from 3.01 pg/2C in the diploid Lamarckia aurea to 20.11 pg/2C in the hexaploid Lolium giganteum (syn. Festuca gigantea). The apparently highest 2C value of 30.48 pg in Festuceae was found in 14x Festuca yvesii (Martínez-Sagarra et al. 2021). The monoploid 1Cx values varied between 1.51 pg/1Cx in Lamarckia aurea and 5.44 pg/1Cx in Festuca lachenalii (syn. Micropyrum tenellum), and the chromosome sizes (MCs) varied between 0.22 pg and 0.78 pg in the same species, respectively (Tables 2, 3; Online Resource 1). The variation of all of these quantitative genome parameters was therefore very similar to that of the tribe Aveneae, as shown in the corresponding plots (Figs. 1, 2).
The tribe Festuceae, comprising about 22 genera and 721 species (Table 1), is almost cosmopolitan, with the center of diversity in the temperate regions of the northern hemisphere, including the Mediterranean region. The tribe extends to the tropics and to the southern hemisphere along the Andes in South America and into the high mountains of Africa.
Representatives of all subtribes of the Festuceae, except of the Ammochloinae, have been studied and are listed below in alphabetical order (Tables 1–3; Figs. 2, 3; Online Resource 1).
Cynosurinae. This subtribe, which has been merged with Parapholiinae into a single subtribe based on molecular phylogenetic analyses (Rasti et al. 2023), comprises 10–11 genera, seven of which have been studied, with a total of eleven species and 41 accessions. The 2C values of this Eurasian group, which evolved mainly in drier regions from the Mediterranean to the Middle East, ranged from 5.74 pg in the diploid Desmazeria sicula (2n = 14) to 17.50 pg (MVs) in the polyploid Parapholis incurva (2n = 38).
The 2C values of all diploids (2n = 14) were only moderately variable (i.e., Catapodium rigidum, Ciliochloa effusa (syn. Cynosurus effusus), Cutandia maritima, Cynosurus cristatus, Desmazeria sicula, Falona echinata (syn. Cynosurus echinatus) and Parapholis filiformis. The minimum values were again found in Desmazeria sicula (2.87 pg/1Cx, MC of 0.41 pg), the maximum values in Falona echinata (7.77 pg/2C, 3.88 pg/1Cx, MC of 0.56 pg (MVs each).
The tetraploid Catapodium marinum (2n = 28) had a 2C value of 12.88 pg, which was about twice that of the diploids. It had a similar 1Cx value of 3.22 pg and an MC of 0.46 pg (MVs each).
Within the genus Parapholis (including Hainardia), the two euploid (orthoploid) species examined were two accessions of the diploid P. filiformis (2n = 14), and one of the tetraploid P. strigosa (2n = 28) with 2C values of 6.59 pg (MV) and 14.73 pg, 1Cx values of 3.30 pg (MV) and 3.68 pg, and MCs of 0.47 pg (MV) and 0.53 pg, respectively. In the dysploid and presumably hypotetraploid P. cylindrica (syn. H. cylindrica) with 2n = 4x = 26 and the presumably hypohexaploid P. incurva (2n = 6x = 38) the 2C values were 14.44 pg and 17.50 pg, and the MCs were 0.56 pg and 0.46 pg, respectively (MVs each). Assuming that these two dysploid species have essentially tetraploid and hexaploid chromosome complements, respectively, the approximate 1Cx values would be 3.61 pg for P. cylindrica and 2.92 pg for P. incurva.
The only taxa in Cynosurinae that have been previously studied with FCM so far appear to bey Cynosurus cristatus with 2C values of 5.38 pg, 5.43 pg, 5.94 pg and 6.10 pg (Šmarda et al. 2008, 2013, 2019; Zonneveld 2019) and P. strigosa with 13.0 pg (Zonneveld 2019). Some of these values are broadly consistent with our results, while the relatively low value of 5.38 pg represents a 1.3-fold difference from our estimate.
Dactylidinae. This small subtribe is widespread in the Old World and includes the consistently perennial genus Dactylis, which is widespread in Eurasia and the Mediterranean, and Lamarckia with the single species L. aurea, a characteristic annual of the Saharo-Sindian and Mediterranean regions. The 2C values of the Dactylidinae ranged from 3.01 pg/2C (MV) in zje diploid L. aurea to 8.52 pg/2C in the tetraploid cytotype of D. glomerata subsp. glomerata.
Dactylis was represented by accessions of diploid D. glomerata subsp. himalayensis and D. polygama and the aforementioned tetraploid D. glomerata subsp. glomerata. The diploids had holoploid genome sizes of 4.31–4.64 pg/2C (MVs), the latter slightly less than twice as much (8.52 pg/2C). The 1Cx values and MCs of both cytotypes were almost identical, i.e. 2.24 pg/1Cx versus 2.13 pg/1Cx, while the MCs were 0.32 pg versus 0.30 pg, respectively. Several previous genome size studies in Dactylis using FCM + PI, and partly recognizing narrowly defined seperate species or using taxonomic synonyms of D. glomerata, are in broad agreement with our results and recorded 2C values of 4.04–4.53pg/2C in diploids (Šmarda et al. 2008, 2019) and 7.81–9.04 pg/2C in tetraploids (Greilhuber and Baranyi 1999; Zonneveld et al. 2005; Šmarda et al. 2008, 2019; Pustahija et al. 2013; Vallès et al. 2014; Zonneveld 2019). However, slightly different and lower values have been reported for the 4x cytotype of D. glomerata without specifying the fluorescent dye (Marie and Brown 1993) and for the 2x, 4x and 6x cytotypes using DAPI (Horjales et al. 1995). In the latter case, the low estimates were most likely caused by the preferential binding of DAPI to AT base pairs, as opposed to the intercalating and not base pair-specific propidium iodide (Sumner 1990).
The second genus of the Dactylidinae, the monospecific Lamarckia with the chromosome number of 2n = 2x = 14 verified in L. aurea, had a 1Cx value of 1.51 pg (see above) and a MC of 0.22 pg (MVs each), thus a significantly smaller genome and chromosome size than D. glomerata.
Loliinae. This large subtribe, which mainly contributes to the almost worldwide distribution of the tribe Festuceae, was sampled by 17 exemplary species of the genus Festuca (including Micropyrum, Psilurus, Vulpia) and six of Lolium, comprising a total of 32 accessions.
The 2C values ranged from 3.32 pg in the diploid annual F. octoflora (syn. V. octoflora) (2n = 14) to 20.11 pg in the hexaploid perennial L. giganteum (syn. F. gigantea) (2n = 42). The 1Cx values varied between 1.66 pg in F. octoflora and 5.44 pg (MV) in F. lachenalii (syn. M. tenellum), also a diploid annual.
Otherwise, the 1Cx values showed a more or less continuous range of variation between 2.07 pg (MV) in the tetraploid F. incurva (syn. P. incurvus) (2n = 28) and 5.15 pg in the diploid L. subulatum (2n = 14), of which two accessions each were examined. These species are both annuals like F. octoflora and F. lachenalii. In contrast, the perennials have intermediate 1Cx values with also continuous variation from 2.23 pg in F. valesiaca to 4.34 pg in F. altissima and 4.85 pg in F. drymeja, all of which are diploids with 2n = 14.
Accordingly, the chromosomes (MC) were comparatively small in the annual F. octoflora (0.24 pg) and F. incurva (MV of 0.30 pg), followed by the perennial F. alpina, F. trachyphylla, F. vaginata and F. valesiaca, which had MCs of 0.32–0.35 pg. The perennials F. altissima and F. drymeja had comparatively large chromosomes (MCs of 0.62–0.69 pg), surpassed only by those of the annuals L. subulatum (MV 0.74 pg) and F. lachenalii (MV 0.78 pg).
Chromosome numbers were verified in this study for eleven Loliinae accessions including what appears to be the first chromosome count for L. subulatum (2n = 14). One of the two accessions of F. heterophylla studied is apparently hexaploid and not tetraploid as is common in this species. The new ploidy level of this accession is suggested by its 2C value, which is about 1.5 times higher than that of the other (18.13 pg/2C versus 11.38 pg/2C).
The 2C values of Festuca species extensively studied by Šmarda et al. (2008), also using FCM + PI, are mostly in good agreement with our results, e.g. 4.25–4.59 pg for F. alpina (4.76 pg in this study), 8.04 pg for F. altissima (8.68 pg in this study), 9.78 pg for F. drymeja (9.69 pg in this study), 16.39 pg for F. heteromalla (16.95 pg in this study), 11.32 pg for tetraploid F. heterophylla (11.38 pg in this study), 4.90 pg for F. vaginata (4.87 pg in this study), 16.98–17.22 pg for Lolium arundinaceum (as F. arundinacea) (16.63 pg in this study), 20.75 pg for L. giganteum (as F. gigantea) (20.11 pg in this study) and 5.51 pg for the diploid cytotype of L. perenne (5.67 pg in this study). However, the 2C values of 5.72 pg and 5.04 pg recorded for L. temulentum (Šmarda et al. 2008, 2019), possibly from the same accession, do not agree with our estimates of 7.98 pg/2C and 8.25 pg/2C for two different accessions, which also confirm a previous estimate of 8.27 pg/2C (Zonneveld 2019).
Furthermore, the 2C values of 9.15 pg for F. heterophylla (Šmarda et al. 2019), 4.02 pg for F. valesiaca (4.45 pg in this study) (Šmarda et al. 2019), 18.3 pg and 17.3 pg for L. giganteum (as F. gigantea) (Zonneveld 2019; Šmarda et al. 2019) (20.11 pg in this study), 15.59 pg, 15.94 pg, 17.45 pg and 17.3 pg for L. arundinaceum (as F. arundinacea) (Arumuganathan et al. 1999; Loureiro et al. 2007; Kopecký et al. 2010; Zonneveld 2019), 5.36 pg, 4.93 pg, 5.26 pg, 5.56 pg and 5.51 pg for diploid L. perenne (Kopecký et al. 2010; Šmarda et al. 2019; Zonneveld 2019; Frei et al. 2021 using a double-haploid plant; Moreno-Aguilar et al. 2022) and 6.4 pg for L. persicum (Moreno-Aguilar et al. 2022) (7.51 pg in this study) are broadly in agreement with our estimates.
Genome sequencing studies have reported lengths of approximately 2.26 Gbp (≈ 4.6 pg/2C), 2.467 Gbp (≈ 5.04 pg/2C) and 2.55 Gbp (≈ 5.21 pg/2C) for the Lolium perenne genome (Byrne et al. 2015; Frei et al. 2021; Nagy et al. 2022), with the more recent data closer to the value of 5.67 pg/2C in this study (see above) than the older ones, and moreover, Frei et al. (2021) also considered the genome size of 5.44 pg/2C determined by them with FCM to be more reliable than their estimate by sequencing (see above).
The above values for L. arundinaceum all refer to the hexaploid cytotype (2n = 6x = 42), which was also examined in this study. However, L. arundinaceum includes additional 4x, 8x and 10x cytotypes (Bulińska-Radomska and Lester 1986). The octo- and decaploids are often referred to as taxonomically distinct varieties or subspecies. Their 2C values were also determined by FCM + PI as 16.22 pg in the octoploid subsp. atlantigenum and 19.70 pg in the decaploid var. letourneuxianum, while the hexaploid subsp. corsicum had 16.62 pg (Ezquerro-López et al. 2017 sub Festuca). The latter agrees well with the values for the two hexaploid accessions we studied (16.54–16.72 pg/2C), although they belonged to the type subspecies/variety.
Considering the whole range of variation in holoploid genome sizes within the subtribe Loliinae, the highest known 2C values occur in dodeca- and tetradecaploid accessions of F. yvesii subspecies from the Iberian Peninsula, namely 26.69 pg and 30.40 pg (MVs calculated here), respectively (Martínez-Sagarra et. al. 2021), while the lowest value seems to be that of the aforementioned diploid F. octoflora, resulting in a total range of 3.25–29.7 Gbp, making the total variation in the Loliinae somewhat larger than previously thought (see Moreno-Aguilar et al. 2022).
Regarding the classification of Festuca and its seperate genera, it can be noted that some of the different infraspecific taxa (sections or subsections) seem to have relatively characteristic and different monoploid genome sizes: The species of F. subsect. Festuca had 1Cx values of 2.23–2.44 pg (2x F. alpina, 6x F. trachyphylla, 2x F. vaginata, 2x F. valesiaca), those of sect. Phaeochloa 4.34–4.85 pg (2x F. altissima, 2x F. drymeja) and those of sect. Aulaxyper 2.83–3.02 pg (6x F. heteromalla, 4x and 6x F. heterophylla, 6x F. rubra).
Considering the taxa often ascribed to Vulpia, the low 2C value of 3.32 pg of the diploid F. octoflora is interesting from a biogeographical and phylogenetic point of view, as this species is one of only four native American species (Soreng et al. 2000 onwards, Stace 2022) of this segregated genus. Vulpia has a center of diversity mainly in the Mediterranean, a region, where 17 of the 21 species occur (Euro + Med PlantBase 2006 onwards; Stace 2022). Compared to F. octoflora, the only American Vulpia species whose genome size has been studied so far, the diploid Old World Vulpias have 1.5–1.9 times larger genomes of about 5.11–6.27 pg/2C (see also Šmarda et al. 2008, 2019; Zonneveld 2019). Festuca octoflora (1.16 pg/1Cx) is taxonomically placed in the same section, i.e. the sect. Vulpia (Stace 2022), as 2x V. bromoides (5.12–5.86 pg/2C and 2.66–2.93 pg/1Cx), 4x V. ciliata (8.28–8.99 pg/2C and 2.06–2.25 pg/1Cx) and 6x V. myuros (12.1–13.86 pg/2C and 2.21–2.31 pg/1Cx). However, molecular phylogenetic data did not support this relationship, but showed that F. octoflora (North America and Cono Sur) and the other American Vulpias studied to date are closer to American fescues than to Old World Vulpia, revealing that traditional Vulpia is polyphyletic and most likely polytopic in origin (Inda et al. 2008; Díaz-Pérez et al. 2014; Minaya et al. 2017). This phylogenetic background is consistent with the large genome size difference between F. octoflora and the Old World Vulpias, so that at least the phylogenetically close South American 2x F. australis (syn. V. australis) may also have a small genome, possibly also the North American 6x F. microstachys (syn. V. microstachys) and other related American taxa (see clade IFL6 of Díaz-Pérez et al. 2014). This suggests that the origin of the American F. octoflora (and its relatives?) was accompanied by a genome size reduction similar to that found in the Old World Vulpia (Šmarda et al. 2008), but even larger.
In the Old World Vulpia, it has been suggested that the origin of the hexaploids is based on allopolyploidy (Stace 2005). In particular, in the case of 6x V. myuros, an origin from a diploid species such as V. muralis and a tetraploid member of the V. ciliata-Psilurus incurvus clade has been suggested. Although the genome size of V. muralis is still unknown, that of the presumably closely related 2x V. bromoides is known to be 5.48 pg/2C, that of 4x V. ciliata 8.70 pg/2C (MV), that of 6x V. myuros 13.60 pg/2C (MV) (Šmarda et al. 2008; Zonneveld 2019) and that of 4x F. incurva (syn. Psilurus incurvus) 8.27 pg/2C. The additivity of the found genome sizes of the diploids and tetraploids (Table 5 in Online Resource 2) would be compatible with the proposed origin of the hexaploid V. myuros, with the caveat that the genome sizes of other potential parental species are not known.
Tribe Poeae and its subtribes
The 2C values ranged from 2.41 pg in the diploid Poa persica to 36.73 pg in the 12x cytotype of Arctagrostis latifolia (Tables 1–3; Figs. 2, 3; Online Resource 1). The apparently smallest 2C value recorded for Poeae to date was 1.49 pg/2C in P. supina (Mao and Huff 2012).
The 1Cx values ranged from 0.75 pg/1Cx in Poa supina, 1.21 pg/1Cx in Poa persica to 6.07/1Cx in Ventenata macra, all of which are diploids with 2n = 2x = 14. Their MCs were accordingly 0.11–0.87 pg.
The largest chromosomes known so far for the Poeae occurred in the dysploid Colpodium biebersteinianum (2n = 2x = 4) with an MC of 0.90 pg (MV), corresponding to the recorded genome sizes of 3.5–3.7 pg/2C (as Zingeria biebersteiniana), also estimated by using FCM + PI (Houben et al. 2003; Kotseruba et al. 2003, 2010).
The variation in all of these quantitative genome parameters therefore was therefore very similar to that of the tribe Aveneae, as shown in the corresponding plots (Figs. 1, 2).
The tribe Poeae comprises about 41 genera and 878 species (Table 1) and is distributed in temperate to often rather cool climates of both the northern and southern hemispheres, reaching the Arctic and Antarctic zones. Ecologically, this group appears to have a wide amplitude. It often grows in wetlands, but also occurs in forests, at high altitudes, and in relatively arid habitats.
Alopecurinae. The almost cosmopolitan Alopecurinae were represented in this study by four species of the genus Alopecurus, which comprises altogether about 44 species worldwide in temperate climates, while the small Asian, mainly Siberian to Arctic genus Limnas was not sampled. The 2C values ranged from 7.16 pg and 7.61 pg in the diploid annuals A. myosuroides and A. aequalis (both 2n = 14), respectively, to 12.73 pg (MV) in a tetraploid accession of the perennial A. pratensis (2n = 28). Alopecurus myosuroides had a 1Cx value of 3.58 pg and an MC of 0.51 pg, which were similar to those of A. pratensis with 3.18 pg/1Cx and an MC of 0.45 pg and those of A. aequalis with 3.81 pg/1Cx and an MC of 0.54 pg.
Previous FCM + PI data for Alopecurus were 7.2 pg/2C and 7.06 pg/2C for A. aequalis, 7.74 pg/2C and 7.52 pg/2C for A. myosuroides, 13.25 pg/2C and 12.6 pg/2C for A. pratensis (Wentworth et al. 2004; Zonneveld 2019), which are in agreement with our results, while those of 5.88 pg in A. aequalis and 6.58 pg in A. myosuroides (Šmarda et al. 2019) and 11.19 pg in A. pratensis (Šmarda et al. 2013) are again lower. Taxonomically, all of our studied species represent at the same time different sections of Alopecurus:
The annual A. cucullatus (syn. Cornucopiae cucullatum) had a higher 1Cx value of 4.32 pg and a higher MC of 0.62 pg (MVs) than the other studied Alopecurus species.
Although our sample of Alopecurus species is limited to four species, three different sections of this genus (Tzvelev 1976; Doğan 1999; Cabi et al. 2017; Tzvelev and Probatova 2019; Gnutikov et al. 2024) are represented: (1) sect. Pseudophalaris with diploid A. myosuroides (3.58 pg/1Cx), (2) sect. Alopecurium with diploid A. aequalis (3.81 pg/1Cx) and (3) sect. Alopecurus with tetraploid A. pratensis (3.18 pg/1Cx) and diploid A. cucullatus (4.32 pg/1Cx), a species that presumably also belongs to this section based on its placement in the trees of molecular phylogenetic analyses (Boudko 2014). The tetraploid A. pratensis has a smaller calculated monoploid genome size than all diploids of Alopecurus examined, especially with respect to the presumably closely related A. cucullatus. This may indeed reflect “genome downsizing” following polyploidization, as previously suggested for the tetraploid A. geniculatus (Wentworth et al. 2004), but molecular phylogenetic data support a possible hybrid origin of the sect. Alopecurus and the presence of different maternal lines in A. pratensis (Cabi et al. 2017; Gnutikov et al. 2024). This implies that allopolyploidy involving the donor of a small genome, rather than genome downsizing, may indeed be the case in A. pratensis (Table 5 in Online Supplement 2), which is supported by meiotic studies suggesting segmental allopolyploidy for this species (Sieber and Murray 1979; Wentworth et al. 2004).
Avenulinae. Avenula pubescens (2n = 2x = 14), the only species of this monospecific subtribe, had a comparatively large holoploid genome (10.67 pg/2C), a 1Cx value of 5.34 pg and an MC of 0.76 pg. A quite similar value of 10.1 pg/2C was found previously (Zonneveld 2019 as Helictotrichon pubescens), while 8.38 pg and 9.40 pg (Šmarda et al. 2013, 2019) again were lower.
Beckmanniinae. This small subtribe with only four genera was studied using the genera Beckmannia, with two species in the Holarctic, and the monospecific Pholiurus from western Eurasia, both characterized by a chromosome number of 2n = 2x = 14. The two Beckmannia species examined, with two accessions each of perennial B. eruciformis and the perennial to annual B. syzigachne, had fairly uniform 2C values of 6.21–6.59 pg, 1Cx values of 3.11–3.30 and MCs of 0.44–0.47 pg (MVs each). For B. eruciformis, 5.37 pg/2C has been previously reported (Šmarda et al. 2019), a value lower than ours. Pholiurus pannonicus had a larger holoploid genome than Beckmannia, with 8.14 pg/2C, a 1Cx value of 4.07 pg and an MC of 0.58 pg (MVs).
Cinninae. The subtribe Cinninae, which comprises the widely distributed Holarctic Cinna and four other genera with disjunct distributions (Gillespie et al. 2022), was represented in this study only by C. arundinacea with a verified 2n = 4x = 28, the 2C value of 12.40 pg, the 1Cx value of 3.10 pg and an MC of 0.44 pg. Three previously examined C. arundinacea accessions were found to have 10.6–10.8 pg/2C (Bai et al. 2012).
For the endemic New Zealand genus Simplicia, a southern hemisphere outlier of the Cinninae, 2C values of 11.07 pg and 10.24 pg were found for two of its three species, S. buchananii and S. laxa, respectively (Murray et al. 2005). Both are tetraploid, therefore their 1Cx values are 2.77 pg and 2.56 pg and the MCs are 0.40 pg and 0.37 pg, all in remarkable similarity to the values of the northern hemisphere Cinna.
Coleanthinae. This subtribe has an almost worldwide distribution and represents a lineage with about ten genera, of which five were studied. The 2C values ranged from 2.71 pg in the annual diploid Coleanthus subtilis to 8.80 pg (MVs) in the hexaploid Puccinellia distans and 9.19 pg in the, judging by the 2C value, also hexaploid P. leiolepis. An octoploid accession of P. maritima had a 2C value of 10.70 pg (Zonneveld 2019), the highest value recorded so far for the subtribe Coleanthinae. However, available chromosome counts indicate that up to 11-fold polyploidy occurs in the Coleanthinae (CCDB), suggesting that the maximum 2C value in this subtribe may actually be somewhat higher.
The subtribe Coleanthinae is characterized by the occurrence of three different chromosome base numbers, x = 2, 5 and 7. The presumably phylogenetically original chromosome base number of x = 7 is found in the genera Puccinellia with five sampled species and Sclerochloa with only one species (monospecific genus). The short-lived small Coleanthus subtilis had the lowest 1Cx value of 1.36 pg and an MC of 0.19 pg, followed by diploid to hexaploid, perennial Puccinellia species with 1Cx values of 1.37–1.53 pg and an MC of 0.20–0.22 pg. The six examined accessions of the diploid annual S. dura had the highest 1Cx values of 1.63 pg and an MC of 0.23 pg (MVs each).
Catabrosa aquatica, characterized by a monoploid chromosome set of x = 5, was represented by one tetraploid and one hexaploid accession, as judged by the 2C values. The 1Cx values were 1.38 and 1.28 pg, the MCs were 0.28 pg and 0.26 pg, respectively.
The two Colpodium (syn. Zingeria) species studied, the diploid C. versicolor (2n = 4) and the tetraploid C. trichopodum (2n = 8) represent the x = 2 lineage. Their 2C values were 2.85 pg and 4.98 pg. The 1Cx value of C. versicolor was 1.43 pg and its MC was 0.71 pg, while these values were not calculated for 4x C. trichopodum, an allopolyploid containing differently sized monoploid genomes (Kotseruba et al. 2003), as detailed below (see chapter “Dysploidy…”). These data are broadly consistent with a previous genome size estimate also using FCM + PI of 2.40 pg/2C in C. versicolor (Kotseruba et al. 2005). However, an estimate of 5.30 pg/2C recorded for C. trichopodum (as Z. trichopoda) (Kotseruba et al. 2003) actually refers to C. pisidica (see Kotseruba et al. 2010), a close relative of C. trichopodum that was previously not considered as separate species. In C. biebersteinianum (2n = 4), genome sizes of 3.5 pg/2C and (recalculated) 3.7 pg/2C, also determined by FCM + PI, were recorded (Kotseruba et al. 2003; Houben et al. 2003), while the hexaploid C. kochii (as Z. kochii) (2n = 12) with 6.96 pg/2C had the largest holoploid genome found so far in the x = 2 lineage (Kotseruba et al. 2010).
The previously reported DNA C-values using FCM + PI for Sclerochloa dura (3.12 pg/2C (Šmarda et al. 2008), Catabrosa aquatica (6.29 pg/2C), Puccinellia distans (8.70 pg/2C) (Zonneveld 2019) as well as P. limosa (2.59 pg/2C for the diploids of Central Germany, taxonomically treated under P. distans agg.) (Kúr et al. 2023) mostly agree well with our data, while the values recorded for Coleanthus subtilis (2.30 pg/2C), Catabrosa aquatica (5.12 pg/2C) and Sclerochloa dura (2.61 pg/2C) (Šmarda et al. 2019) are again up to 20% lower.
Hookerochloinae. Hookerochloinae comprises five genera, four of which are scattered throughout the southern hemisphere (Gillespie et al. 2022), while the genus Arctagrostis is widespread in the Holarctic and comprises only two species, one of which was sampled. Different specimens examined from a collection of A. latifolia from Alaska had 2C values of 24.51 pg and 36.73 pg. Their chromosome numbers of 2n = 56 and 2n = 84, representing the octo- and dodecaploid levels, respectively, were determined by counting. The 1Cx value of 3.06 pg and the MC of 0.44 pg were identical in both.
Miliinae. The monogeneric subtribe Miliinae has about six species in the Holarctic region, of which the perennial forest species Milium effusum, the most widespread species, was sampled. This tetraploid species (2n = 28), as verified by chromosome counting, had a 2C value of 9.23 pg, a 1Cx of 2.31 pg and an MC of 0.33 pg. The 2C value is comparable to one of the previous FCM + PI estimates, namely 8.30 pg/2C (Šmarda et al. 2019), while 5.17 pg/2C could refer to a diploid accession of M. effusum (Zonneveld 2019) or another species, as this accession is also listed as ‘M. cf. effusum’ (Zonneveld 2019: Electr. Suppl. Table 5). However, the four other estimates of 8.86–9.22 pg/2C for M. effusum (Zonneveld 2019: Electr. Suppl. Table 5) fit our result much better and should therefore actually belong to the tetraploid M. effusum.
For M. vernale, an annual diploid with the diverging dysploid chromosome number of 2n = 8, a 2C value of 6.28 pg was previously obtained (Zonneveld 2019), which largely agrees with the value of 5.73 pg previously estimated by Feulgen densitometry (Bennett and Thomas 1991). The 1Cx value of M. vernale would therefore be 3.14 pg and the MC would be 0.79 pg (Table 3).
Phleinae. The Holarctic and Andean subtribe Phleinae comprises only the genus Phleum, when its segregate genera are included. Of the approximately 15 species, nine were examined in this study, including diploid (2n = 14) to hexaploid (2n = 42) species and cytotypes, with 2C values ranging from 2.77 pg in the diploid P. cf. alpinum to 8.74 pg (MV) in the hexaploid P. pratense. The other diploids studied (P. arenarium, P. bertolonii, P. phleoides, P. rhaeticum, P. subulatum) had 3.13–4.01 pg/2C, the tetraploids P. alpinum and P. paniculatum had 6.08–6.55 pg/2C. Previous FCM + PI genome size data mostly agree well with these values, e,g. 8.1 pg/2C in presumably also 6x P. pratense (Bai et al. 2012), 2.44 pg/2C in 2x P. alpinum, 3.49 pg/2C in 4x P. phleoides and 7.63 pg–7.99 pg/2C in 6x P. pratense (Šmarda et al. 2013, 2019), 3.13 pg/2C for 2x P. arenarium and 8.47–8.78 pg/2C in 6x P. pratense (Zonneveld 2019). Several diploid accessions of P. alpinum (as P. commutatum) and P. rhaeticum had mean 2C values of 2.68 pg and 2.61 pg, respectively, while tetraploid P. alpinum accessions (as P. commutatum) from America and Europe uniformly had about 6.14–6.20 pg (Kula et al. 2006), which is also consistent with our results. The 1Cx values of all Phleum accessions examined in our study varied between 1.39 pg and 2.01 pg, the MCs between 0.20 and 0.29 pg.
The subtribe Poinae consists only of the genus Poa, into which several segregate genera have been included. Poa is the largest of all grass genera with about 570 species (Soreng et al. 2022). It is distributed almost worldwide, even reaching the Arctic region and Antarctica.
The total of 24 Poa species and 35 accessions analyzed had 2C values ranging from 2.41 pg in a probably diploid accession of P. persica (syn. Eremepoa persica) to 15.29 pg in 12x–14x P. cita. Two other diploids examined, P. chaixii and P. trivialis, had 3.15 and 3.36 pg/2C, respectively.
The tetraploids examined were the annual P. annua with 4.40 pg/2C (MV) and the perennials P. badensis with 6.01 pg/2C (MV) and P. billardierei (syn. Austrofestuca littoralis) with 7.74 pg/2C. Chromosome counts were not performed for the other species and accessions, but they all were probably polyploid. Also P. palustris with only 4.49 pg/2C (MV) was probably not diploid but tetraploid (CCDB).
The 1Cx values of the diploids ranged from 1.21 to 1.68 pg, those of the polyploids from 1.07 to 1.94, the latter in the Australasian P. billardierei. The MCs showed similar variation, ranging from 0.17 to 0.24 pg in diploids, and from 0.15 to 0.28 pg in polyploids, the latter value again in P. billardierei, indicating that this species has the largest monoploid genome and the largest chromosomes found in Poa so far.
For most of the species we studied, the 1Cx values and MCs cannot be given because the chromosome numbers of these accessions are unknown. All these species are known to have more than one cytotype (ploidy level), some show in additional aneuploidy and facultative apomixis (e.g. P. bulbosa and P. pratensis), and it is not clear which cytotypes our accessions used for genome size estimates belong to.
Our genome size data in Poa are largely consistent with the results of several previous studies also using FCM + PI, such as 4.21 pg/2C and 4.19 pg/2C (MVs) recorded for P. annua (Zonneveld 2019; Siljak-Yakovlev et al. 2020), 14.71 pg/2C (MV) for P. cita and 7.42 pg/2C for P. billardierei (as Austrofestuca littoralis) (both Murray et al. 2005), 4.2–5.8 pg/2C for P. palustris (Bai et al. 2012), and 3.32 pg/2C and 3.48 pg/2C for P. trivialis (Šmarda et al. 2008; Zonneveld 2019). For P. annua, 3.88 pg/2C (Mao and Huff 2012) and 3.87 pg/2C (Šmarda et al. 2019) have also been found, which is closer to the genome size of 1,778 Mbp (≈ 3.64 pg/2C) suggested by sequencing (Robbins et al. 2023; Benson et al. 2023) than our FCM estimate. The same discrepancy between FCM- and sequencing-based genome size estimates as in P. annua is evident in the case of P. trivialis, whose sequenced genome length of 1,350 Mbp (≈ 2.76 pg/2C) (Brunharo et al. 2024) is about 17–21% lower than the available FCM genome size estimates of 3.32–3.48 pg/2C (see above).
However, it is difficult to say which are the more accurate estimates, since sequencing approaches, although depending on the method used, typically underestimate the true genome size due to the mostly incomplete representation of the amount of repetitive DNA as has often been noted (e.g. Bennett et al. 2003; Doležel et al. 2018; Kapustová et al. 2019; Blommaert 2020; Becher et al. 2022; Winterfeld et al. 2023, 2024; Tkach et al. 2024), due to the difficulty of assembling repeats with monomers longer than the sequencing read length and the representation of two or more copies, which may be of different lengths, in only one sequence of the assembly.
Some other previously reported C-values of Poa species were also about 6–18% lower than ours, such as 2.74 pg/2C and 2.97 pg/2C instead of 3.36 pg/2C for P. trivialis (Wieners et al. 2006: p. 1537; Šmarda et al. 2013), 4.38 pg/2C instead of 5.13 pg/2C for P. angustifolia, 2. 83 pg/2C instead of 3.15 pg/2C for P. chaixii, 5.28 pg/2C instead of 6.25 pg/2C for P. nemoralis and 3.73 pg/2C instead of 4.49 pg/2C for P. palustris (Šmarda et al. 2019), while these samples, either verified by chromosome counting or most likely, have the same ploidy level as the accessions we used.
In other cases, different ploidy levels were investigated, resulting in 2.48 pg/2C for diploid (Šmarda et al. 2019) and 6.01 pg/2C for tetraploid P. badensis, and 4.08 pg/2C for tetraploid P. laxa (Šmarda et al. 2019), suggesting that our accessions of P. laxa was octoploid, judging by its genome size of 8.12 pg/2C.
The three studied accessions of P. pratensis with a genome size of 4.93–7.61 pg/2C had probably different chromosome numbers. The numerous cytotypes of this facultatively apomictic species (Bonos and Huff 2013), comprising 4x–22x ploidy levels, cause considerable variation in the size of the holoploid genomes (Eaton et al. 2004; Wieners et al. 2006; Dennhardt et al. 2016; Zonneveld 2019; Ghanbari et al. 2023; Phillips et al. 2023). Using FCM + PI, an octoploid accession of P. pratensis subsp. angustifolia was found to have 3,525.29 Mbp (≈ 7.21 pg/2C), while another accession of unknown ploidy of subsp. angustifolia had 3,248.04 Mbp (≈ 6.64 pg/2C) and of subsp. pratense had 4,030.61 Mbp (≈ 8.24 pg/2C) and 4,856.38 Mbp (≈ 9.93 pg/2C), respectively (Phillips et al. 2023).
The highest value of 18.66 pg/2C recorded so far in the genus Poa also belongs to P. pratensis (Raggi et al. 2015, using DAPI as fluorescent dye), while for most other accessions of P. pratensis have 2C values below 11 pg/2C have been recorded.
The lowest 2C values recorded in this genus refer to the perennial P. supina and the annual P. infirma, which are the diploid parental species of the aforementioned weedy annual to perennial, cosmopolitan allotetraploid P. annua (Nannfeld 1937 using the synonym P. exilis for P. infirma; Tutin 1952, 1957; Soreng et al. 2010; Mao and Huff 2012; Chen et al. 2016; Nosov et al. 2019; Benson et al. 2023).
The genome size of P. infirma was 2.52 pg/2C and 2.84 pg/2C (MVs) according to FCM data (Mao and Huff 2012; Zonneveld 2019) and was 1,125.5 Mbp (≈ 2.30 pg/2C) according to genome sequencing (Benson et al. 2023; Robbins et al. 2023). Poa supina had an even smaller genome size estimated as 1.48 pg/2C and 1.49 pg/2C by FCM (Mao and Huff 2012; Šmarda et al. 2019) and 636 Mbp (≈ 1.30 pg/2C) by DNA sequencing (Benson et al. 2023), respectively.
Unclear are some other genome size estimates for Poa species (Joshi et al. 2016), which were also performed using FCM + PI. For all eight species considered diploid (P. asiae-minoris, P. badensis, P. bucharica, P. chaixii, P. hybrida, P. pumila, P. sibirica, P. trivialis with two accessions), the identical 2C value of 2.5 pg was found, which is hardly credible.
Ventenatinae. This subtribe comprises six genera and is distributed in western Eurasia, the Mediterranean and the Middle East. Representatives of three genera were examined, all them diploid (2n = 14). The perennial Bellardiochloa violacea had the smallest holoploid genome size of 7.77 pg/2C, a 1Cx value of 3.89 pg and an MC of 0.56 pg.
Two species and three accessions of the annual genus Apera had 8.68–9.60 pg/2C, 1Cx values of 4.34–4.80 pg and MCs of 0.62–0.69 pg (A. interrupta and A. spica-venti). Previously recorded genome sizes were 8.66 pg/2C (Šmarda et al. 2019) and 10.3 pg/2C for A. spica-venti and 8.73 pg/2C for A. interrupta (both Zonneveld 2019), in agreement with our results.
The also annual Ventenata, in which two species and four accessions were studied, showed some variation, as V. dubia had 10.65 pg/2C, a 1Cx of 5.33 pg and an MC of 0.76 pg (MVs), while V. macra had a larger holoploid genome of 12.14 pg/2C, a 1Cx of 6.07 pg and an MC of 0.87 pg.
Arctopoa. This hybrid taxon of the tribe Poeae, has ancestors from the subtribes Poinae and Cinninae each (Gillespie et al. 2010, 2022; Tkach et al. 2020). The studied hexaploid Arctopoa eminens (2n = 42) had a genome size of 11.75 pg/2C and an MC of 0.28 pg. The 1Cx value of 1.96 pg and the MC of 0.28 pg are average values calculated for this allopolyploid, but the combination of very different sized monoploid genomes is also possible in this species.
Intertribe hybrid groups
The subtribes, which are taxonomically unplaced due to presumable hybrid origin or not yet assigned to a tribe, had 2C values ranging from 2.86 pg (MV) in the diploid Corynephorus canescens to 38.89 pg in 18x Helictochloa pratensis (Tables 2, 3; Figs. 1, 2A; Online Resource 1.
The largest holoploid genome size among the diploids with x = 7 was found in Mibora minima with 10.05 pg/2C (MV). In the diploid Echinaria capitata, characterized by the dysploid number of x = 9, the 2C value was even higher and amounted to 16.59 pg/1Cx (MV) (Table 2, 3; Fig. 1, 2B).
Monoploid genome sizes ranged from 1.13 pg/1Cx in a tetraploid of Holcus mollis accession to 5.03 pg/1Cx in Mibora minima, both with x = 7, exceeded by 8.30 pg/1Cx (MVs) in Echinaria capitata with x = 9. Similarly, the MCs varied between the minimum of 0.16 pg in H. mollis and the maximum of 0.72 pg in M. minima and 0.92 in E. capitata, respectively (Fig. 1, 2C).
The size distribution of 2C values, 1Cx values, and MCs of the intertribe hybrid groups as a whole closely resembles that of the tribe Aveneae (Fig. 1), which is consistent with the involvement of Aveneae as one of the parental taxa in their origin, as shown in Table 1, which updates the parental tribe designations for two subtribes compared to Tkach et al. (2020: p. 255).
The intertribe hybrid groups comprise 18 genera and approximately 159 species, representing seven subtribes, six of which were sampled in this study (Table 1). Most of their genera are distributed in the temperate regions of Eurasia and the Mediterranean, only a few are Holarctic or nearly so (Helictochloa, Scolochloa, Vahlodea), while Dryopoa is SE Australian, and only Deschampsia is nearly cosmopolitan. The subtribes of the intertribe hybrid groups are listed below in alphabetical order (Tables 1–3; Figs. 2, 3; Online Resource 1).
Airinae. This subtribe with four genera and nearly worldwide distribution had 2C values ranging from 2.86 pg in the diploid perennial Corynephorus canescens to 11.95 pg in the tetraploid perennial Avenella flexuosa (MVs each). The consistently annual Aira species studied had 2C values of 5.72 pg (MV) in A. elegans and 6.92 pg in A. praecox (MV), both of which were diploid. Tetraploid A. caryophyllea had 12.45 pg/2C (MV). The 1Cx values and MC of the Airinae species varied accordingly, with the lowest values occurring in C. canescens (1Cx value of 1.43 pg, MC of 0.20 pg) and the highest in the diploid A. praecox (1Cx value of 3.46 pg, MC of 0.49 pg). Due to its allopolyploid (amphidiploid) status (Albers 1973, 1978, 1980a, b), the 1Cx and MC were not calculated for A. caryophyllea (Table 5 in Online Resource 2).
Genome size data from most previous FCM studies were largely in agreement with our estimates and amounted to values of 6.12 pg/2C and 6.78 pg/2C for Aira praecox, 2.24 pg/2C and 3.07 pg/2C for C. canescens, 10.87 pg/2C and 12.20 pg/2C for Avenella flexuosa (as Deschampsia flexuosa) (Šmarda et al. 2019; Zonneveld 2019). Aira caryophyllea has been reported with 12.70 pg/2C (Zonneveld 2019), which agrees with our estimate, as well as the previous Feulgen-densitometric estimates of 5.87 pg for A. praecox and 12.05 pg for A. caryophyllea (Albers 1980a).
Using the non-intercalating DAPI as a fluorescent dye for FCM, the ratio between the genome sizes of A. praecox and A. caryophyllea was 0.56–0.59 pg, calculated from the data of Gregor et al. (2023), which is consistent with the 0.56 pg (MV) found in our estimates.
Aristaveninae. The genus Deschampsia, with about 62 species, is the only representative of the subtribe, which is distributed almost worldwide, including Antarctica. Deschampsia is widespread in the Holarctic, but also has a center of species diversity in South America. The 2C values varied from 7.58 pg in D. danthonioides to 18.61 pg in D. littoralis (MVs). All but one of the species studied chromosomally so far have 2n = 26 or 52 (CCDB) or slightly different chromosome numbers due to aneuploidy or B chromosomes (Albers 1972), which at first glance implies a chromosome base number of x = 13.
The only exception to x = 13 is found in D. setacea, a rare species from the comparatively humid regions of western Europe, which is diploid and has 2n = 14. It was split off at some point as the genus Aristavena. All other species with x = 7 that were occasionally listed under Deschampsia (e.g., in the CCDB), such as D. atropurpurea, D. flexuosa and D. minor actually belong to other genera, namely Avenella, Holcus and Vahlodea. However, based on molecular phylogenetic data (Tkach and Röser, unpublished data), D. setacea forms a maximally supported monophyletic group with all other Deschampsia species, making it plausible to accept its inclusion within Deschampsia. The 2C value of D. setacea was 5.25 pg (MV), the 1Cx value 2.63 pg and the MC 0.38 pg.
The Deschampsia species with 2n = 26 are either hypotetraploids, i.e. tetraploids with 2n = 28, in which two chromosomes have been lost in some way, e.g. by centric fusion of two acrocentric chromosomes, or they are allopolyploids derived from a diploid ancestor with x = 7 such as D. setacea and an unknown diploid with x = 6 (Albers 1980a, Garcia Suarez et al. 1997; Gonzales et al. 2021). Deschampsia argentea from Madeira (three accessions studied), the D. cespitosa (five accessions), D. media (one accession) and D. wibeliana (one accession specifically from the Elbe estuarine zone) from central Europe had 8.63–9.96 pg/2C, 1Cx values of about 2.16–2.49 pg when considered tetraploid and an MCs of 0.33–0.38 pg. Deschampsia koelerioides from China had a larger genome of 11.24 pg/2C, a 1Cx of 2.81 pg and an MC of 0.43 pg (MVs each).
Deschampsia littoralis from northwestern Switzerland and D. rhenana from the Lake Constance are both octoploids with 2n = 52 (CCDB). They had 18.61 and 18.30 pg/2C, 2.33 and 2.29 pg/1Cx and MCs of 0.36 and 0.35 pg, respectively (MVs for D. littoralis). The 1Cx and MC values were similar to those of the 2n = 26 species mentioned above.
The annual North American D. danthonioides has a smaller holoploid genome, 7.58 pg/2C, than the other Deschampsia species with 2n = 26. Its 1Cx value was 1.90 pg and the MC was 0.29 pg (MV each).
Previous estimates of 2C genome size using FCM in Deschampsia were 4.94 pg for D. setacea with 2n = 14 (Zonneveld 2019). Considering the taxa with 2n = 26, the 2C values were 9.51 pg for D. argentea (Greimler et al. 2022), 7.51–10.88 pg for D. cespitosa (Murray et al. 2005; Šmarda et al. 2019; Zonneveld 2019; Greimler et al. 2022 including many accessions) and 9.73–10.47 pg for D. koelerioides (Greimler et al. 2022). In the taxa with 2n = 52, the 2C values were 17.95–18.08 pg for D. littoralis and 16.96–17.72 pg for D. rhenana (Greimler et al. 2022). All these estimates are in good agreement with our data. Some variation in genome size in D. cespitosa could also be due to B chromosomes, which occur in variable numbers and sizes in this species (Albers 1972). In addition, occasional accessions of D. cespitosa with 2n = 52 have also been found, with holoploid genome sizes of 15.89–17.95 pg/2C (Greimler et al. 2022).
Helictochloinae. The genus Helictochloa, all species of which are perennial and characterized by a chromosome base number of x = 7), forms the subtribe Helictochloinae together with the dysploid Mediterranean annual Molineriella (x = 4; not studied). Helictochloa is widespread in the northern hemisphere and is characterized by a wide range of ploidy levels from 2x to 22x (e.g., Gervais 1973 as Avenochloa; Röser 1996, 1998 as Helictotrichon). The 2C values of the six studied species ranged from 7.31 pg in the diploid Helictochloa versicolor (2n = 14) to 38.89 pg in the highly polyploid H. pratensis (2n = 18x = 126). The 1Cx values varied between 3.66–4.29 pg in diploids and 2.16–3.07 in polyploids, which had correspondingly smaller chromosome sizes (an MC of 0.52–0.61 vs. 0.31–0.44 pg). For 18x H. pratensis accessions, previous estimates using FCM were 30.29–33.22 pg/2C and 35.8 pg/2C (Šmarda et al. 2019 and Zonneveld 2019 as Helictotrichon pratense), which are in the same order of magnitude as our data.
Holcinae. Both genera of this small subtribe were studied, the mainly European genus Holcus and the mainly amphi-Beringian and arctic Vahlodea. The 2C values were 3.36–7.5 pg in the two Holcus species studied and 6.32 pg (only one measurement made) in V. atropurpurea. Three accessions of the diploid H. lanatus examined had 3.36 pg/2C, 1.68 pg/1Cx, and MC of 0.24 pg (MV each). These values are consistent with 3.41 pg/2C (Zonneveld 2019), while 2.97–2.99 pg/2C are slightly lower values (Šmarda et al. 2013, 2019). A verified tetraploid H. mollis accession had 4.53 pg/2C, a 1Cx value of 1.13 pg and an MC of 0.16 pg. Another accession had 7.50 pg/2C, indicating higher polyploidy. Holcus mollis is known to be 3x–7x (CCDB), which explains why values of 5.42 pg/2C and 5.99 pg/2C were also previously found for this species (Šmarda et al. 2019; Zonneveld 2019). Vahlodea atropurpurea has been reported as diploid, but there seems to be only one chromosome available for this species (CCDB), while our accession is most likely polyploid, considering the high 2C value.
Scolochloinae. This subtribe is characterized by a remarkably disjunct distribution of its only two genera, the northern hemisphere genus Scolochloa with two species and the southeastern Australian monospecific genus Dryopoa. The tetraploid S. festucacea (2n = 28) had 9.75 pg/2C, 2.44 pg/1Cx and an MC of 0.35 pg (MVs each). The hexaploid S. marchica, a rare local endemic of eastern Germany (2n = 42) had 14.80 pg/2C, 2.47 pg/1Cx and an MC of also 0.35 pg (MVs). The genome size data suggest that S. marchica is autopolyploid and arose from S. festucacea by fusion of reduced and unreduced gametes, without the involvement of another parental taxon (Röser et al., unpubl. data).
The Australian Dryopoa dives has a much larger holoploid genome of 30.76 pg/2C. The sampled accession was verified to be decaploid (2n = 70), resulting in a comparatively high 1Cx value of 3.08 pg and an MC of 0.44 pg, consistent with the large size of the chromosomes under the microscope (not shown).
Seslerieae. This European to Mediterranean subtribe was represented in this study by all of its six genera. Five of them have a monoploid chromosome set of x = 7, namely the annual, species-poor genus Mibora and the perennial, also species-poor genera Oreochloa, Psilathera and Sesleriella as well as the perennial but species-rich genus Sesleria with about 36 species. The 2C values ranged from 6.81 pg in the diploid P. ovata to 19.87 pg in the hexaploid Sesleria comosa.
x = 7. Mibora minima, a diploid short-lived, ephemeral and winter annual species (2n = 14) had a 2C value of 10.05 pg, a 1Cx of 5.03 pg and an MC of 0.72 pg (MVs), confirming a previous estimate of 9.92 pg/2C (Zonneveld 2019). The perennial diploids (2n = 2x = 14) Oreochloa, Psilathera and Sesleriella had 9.02 pg/2C (MV), 6.81 pg/2C and 8.10 pg/2C, respectively, with 1Cx values of 4.51 pg (MV), 3.41 pg and 4.05 pg, and MCs of 0.64 pg (MV), 0.49 and 0.58 pg, respectively. The genome size of P. ovata (as Sesleria ovata) was previously reported to be 5.95 pg/2C, also determined by FCM + PI (Lazarević et al. 2015).
Interestingly, M. minima is a diminutive, short-lived annual of the winter-mild regions of Western Europe that germinates already in the fall and survives the winter as a plantlet. Its large genome size may be related to such a life form and the ecological characteristics as discussed for the ‘neotenic’ centrolepids of the family Restionaceae and other examples mentioned by Winterfeld et al. (2023).
The consistently polyploid genus Sesleria was represented by seven species and nine accessions. The four tetraploid species examined had 8.62–10.14 pg/2C, the three octoploid species 16.37–19.78 pg/2C. The values of 9.26 pg/2C, 9.78 pg/2C (MV) and 9.55 pg/2C for tetraploid S. albicans (Lysák and Doležel 1998 and Lysák et al. 2000, also using FCM + PI; Zonneveld 2019) are largely in agreement with our estimate of 10.14 pg/2C (MV). For the tetraploids S. alba and S. argentea we obtained 9.04 pg/2C (MV) and 8.62 pg/2C, which is consistent with previous estimates of 8.58 pg/2C and 8.97 pg/2C (MV), respectively (Lazarević et al. 2015). Similarly, our estimate of 18.49 pg/2C in the octoploid S. sadleriana is broadly consistent with that of 18.00 pg/2C (MV) (Lysák and Doležel 1998).
The 1Cx values of all Sesleria species examined ranged from 2.05 pg to 2.54 pg (MVs), the MC from 0.29 pg to 0.36 pg (MVs).
x = 9. The remaining genus of the subtribe Sesleriinae, the annual monospecific Echinaria with 2n = 18 has the divergent monoploid chromosome number of x = 9. Echinaria capitata had a comparatively high 2C value of 16.59 pg, a 1Cx of 8.30 pg and an MC of 0.92 pg (MVs each).
Are there any trends in the evolution of Poodae genome sizes?
Compared to the tribes Aveneae and Festuceae, the Poeae are not fundamentally different with respect to the genome size parameters (Fig. 1), but have a comparatively high proportion of taxa with low 2C values < 5 pg/2C, small monoploid genomes < 2.0 pg and show the prevalence of small chromosomes with MCs < 0.3 pg (see Figs. 1, 2, 3). This is mainly caused by the monogeneric subtribes Phleinae (Phleum) and Poinae (Poa), and most of the Coleanthinae, while one of the late-diverging lineages of the tribe Poeae has mostly large 2C values, but also large monoploid genomes (1Cx values) and large chromosomes (MCs). This lineage consists of the subtribes Alopecurinae, Beckmanniinae, Cinninae, Hookerochloinae and Ventenatinae (Fig. 3), forming the ‘ABCV clade’ (Tkach et al. 2020), which was subsequently renamed the ‘Alopecurinae superclade’ (Gillespie et al. 2022). However, this genome expansion is also seen in the monogeneric subtribe Avenulinae (Avenula), a relative of both, Phleinae/Poeae and the ABCV clade (see also Tkach et al. 2020: Fig. 7).
Regarding the tribe Festuceae, comparatively large 2C values and monoploid genome sizes and MCs are found in the subtribe Cynosurinae/Parapholiinae compared to the Dactylidinae and the Loliinae, a subtribe that apparently diverged earlier in the phylogeny. This could also imply a ‘secondary’ genome and chromosome enlargement, but the variation within Loliinae is strong and there is no clear trend (Fig. 3).
The variation of genome and chromosome sizes in the tribe Aveneae is also very striking and there seems to be no clear direction of ‘evolutionary progress’. However, the enormous variability of these parameters within the monogeneric subtribe Anthoxanthinae (Anthoxanthum) (Tables 2, 3; Figs. 2, 3) is particularly remarkable (Chumová et al. 2015, 2017, 2021). This is also true for the occurrence of the relatively small genome and chromosome sizes in the subtribe Torreyochloinae (Tables 2, 3; Figs. 2, 3), which show a surprisingly uniform monoploid genome size (1Cx value) and MC in the studied species of the diploid Torreya pallida from Canada and the hexaploid Amphibromus nervosus from Australia, which are almost maximally geographically separated.
The variability of genomic parameters appears to be associated less with long-term phylogenetic and evolutionary events and possibly more with evolutionary events within smaller groups of related taxa or individual genera. Features such as polyploidization, which per se is associated with a sudden change in genome size, but also chromosomal mutations associated with dysploid changes in chromosome number, or the evolution of different life forms (annual vs. perennial) have been discussed previously as possible factors for changes in genome size and associated traits and biological properties of organisms (Greilhuber and Leich 2013; Šmarda et al. 2019; Liddell et al. 2021; Zhan et al. 2021; Heslop-Harrison et al. 2023). The available data on numerous representatives of a relatively narrowly defined group of plants, such as the investigated subtribes grouped in the supertribe Poodae, provide good opportunities to address these questions.
Polyploidy (whole-genome duplication, WGD) and genome size
Polyploidy, the multiplication of chromosome sets by whole genome duplication (WGD) as it is commonly referred to today, represents the most significant and drastic change in holoploid genome size (2C values) because it takes effect immediately. In contrast, changes caused by a gradual increase or decrease in the proportion of certain, typically repetitive elements in the DNA are much slower. Since all angiosperms have undergone one or more cycles of genome duplication in their evolutionary past, holoploid genome sizes should in principle have become larger and larger (Leitch and Bennett 2004; Bennett and Leitch 2005; Wendel 2015; Wendel et al. 2018; X. Wang et al. 2021). However, since many angiosperms living today have very small genomes, there must inevitably be mechanisms that counteract this genome expansion through polyploidy, which is very well documented by numerous examples. Within the investigated genome sizes of grasses, there are polyploids whose holoploid genome represents the sum of the constituent individual monoploid genomes (‘additivity’) as well as those in which a reduction of the holoploid genome (‘downsizing’) has already occurred (see Table 5 in Online Resource 2 for all calculations). Both processes can occur in both autopolyploids and allopolyploids, providing new insight into the long-standing debate as to whether alloploidy promotes downsizing, whereas autoploidy is associated with genome size additivity (Ozkan et al. 2001; 2006; Soltis et al. 2003; Leitch and Bennett 2004; Bennett and Leitch 2005; Johnston et al. 2005; Garnatje et al. 2006; Suda et al. 2007; Vaio et al. 2007; Eilam et al. 2008, 2009, 2010; Leitch and Leitch 2008; Pellicer et al. 2010; Husband et al., 2013; Zenil-Ferguson et al. 2016; Becher et al. 2021; X. Wang et al. 2021; Feldman and Levy 2023b).
Additivity of monoploid genome sizes in polyploids. Polyploids whose holoploid genome sizes reflect the sum of the genome sizes of their diploid or lower polyploid relatives document a lack of or ineffective, genome size reduction after polyploidization. They can be autopolyploids or allopolyploids.
Examples found of the most likely autopolyploid taxa with little or no downsizing include (see Table 6 in Online Resource 2):
(A) triploid (3x) Koeleria vallesiana,
(B) tetraploid (4x) Anthoxanthum nitens (syn. Hierochloe nitens), Briza media (as previously noted also by Murray 1975, 1976), Catapodium marinum, Dactylis glomerata subsp. glomerata, Festuca trachyphylla, F. fasciculata (syn. Vulpia fasciculata), Gastridium ventricosum, G. phleoides, Koeleria macrantha, K. spicata (syn. Trisetum spicatum), Poa badensis and P. laxa,
(C) hexaploid (6x) Amphibromus nervosus, Festuca heteromalla, F. heterophylla, F. rubra and Scolochloa marchica,
(D) octoploid (8x) Sesleria sadleriana,
(E) hypo-tetraploid (4x-2) Deschampsia argentea, D. cespitosa, D. koelerioides, D. media and D. wibeliana and
(F) hypo-octoploid (8x-4) Deschampsia littoralis, D. rhenana in comparison with (E).
Examples of allopolyploids with little or no downsizing, characterized by an almost invariant additivity of the different parental genomes include (see Table 7 in Online Resource 2):
(G) tetraploid (4x) Aira caryophyllea, Anthoxanthum odoratum, Poa annua and
(H) hexaploid (6x) Festuca muyros (syn. Vulpia myuros).
Auto- or allopolyploids (unresolved) with little or no downsizing are (see Table 8 in Online Resource 2):
tetraploid (4x) Phleum alpinum,
(J) hexaploid (6x) Puccinellia distans and P. leiolepis,
(K) octoploid (8x) Sesleria comosa.
Polyploids with significant genome downsizing and decrease in the mean chromosome DNA content. Genome downsizing (2C and 1Cx values) and decrease of chromosome size (MC) (MC) > 10% were found in some most likely autopolyploids (Table 9 in Online Resource 2):
(L) tetraploid (4x) Koeleria litvinowii and K. pyramidata,
(M) hypo-octoploid (8x-4) Deschampsia littoralis and D. rhenana in comparison with D. setacea (2x),
(N) decaploid (10x) Helictochloa agropyroides, Koeleria pyramidata, Pentapogon crinitus and P. micranthus.
Examples of allopolyploids with significant genome downsizing and decrease of MC include (see Table 10 in Online Resource 2) are:
(O) tetraploid (4x) Alopecurus pratensis, Colpodium pisidicum/C. trichopodum (syn. Zingeria pisidica/Z. trichopoda), Holcus mollis, Phleum paniculatum,
(P) hexaploid (6x) Phleum pratense,
(Q) hypohexaploid (6x-4) Parapholis incurva,
(R) octodecaploid (18x) Helictochloa praeusta, H. pratensis subsp. pratensis.
Taxa with unresolved status of their polyploidy, whether auto- or allopolyploid, that have a significant downsizing, include (see Table 11 in Online Resource 2):
(S) tetraploid (4x) Anthoxanthum odoratum, Koeleria litvinowii,
(T) hexaploid (6x) Anthoxanthum monticola (syn. Hierochloe alpina),
(U) octoploid (8x) Sesleria tenuifolia.
Polyploids with uncertain diploid ancestors. In some cases, the genome size data help to exclude certain diploid taxa as ancestors or genome donors of some diploids (Tables 5, 12 in Online Resource 2):
The genome size of 1.66 pg/1Cx and the MC of 0.24 pg of Festuca octoflora (syn. Vulpia octoflora) suggest that this New World diploid species was not involved in the origin of the sampled Old World 4x F. fasciculata (syn. V. fasciculata), which had 3.26/1Cx and an MC of 0.47 pg, i.e. approximately twice the chromosome size of F. octoflora. Comparable size data, however, were found in the Mediterranean diploid F. alopecuros (syn. V. alopecuros), i.e. 3.14 pg/2C and an MC of 0.45 pg, thus belonging to the putative diploid ancestors of the Old World polyploids.
The highly polyploid western Mediterranean Helictotrichon filifolium subsp. filifolium (12x) had 2C, 1Cx and MC values that almost exactly matched the parameters of the diploid Helictotrichon species studied, H. decorum and H. sedenense subsp. sedenense, and could therefore considered a good example of allopolyploidy without significant genome downsizing. However, the actual genome donors have been shown to be different species, namely H. parlatorei and H. sarracenorum or related taxa (Wölk et al. 2015; Winterfeld et al. 2016), for which genome size data are not yet available.
A similar case is the allohexaploid H. sempervirens from the Western Alps, which even showed an increase in 2C and 1Cx values and MC compared to the above diploid species, but again, the relevant diploid species that acted as actual genome donors, H. parlatorei and H. setaceum or related taxa (Wölk et al. 2015; Winterfeld et al. 2016), have not been studied.
In the genus Parapholis, the sampled diploid P. filiformis apparently did not represent a genome donor for the tetraploid P. strigosa and the hypotetraploid (4x-2) P. cylindrica due to its genome size parameters with too low a 1Cx value and too low a MC.
Dysploidy and change in genome and chromosome size
Dysploidy is a less dramatic change in chromosome number than polyploidy and usually refers to one or a few chromosomes of a monoploid set of chromosomes causing a gradual change in chromosome number, typically as a decrease, termed descending or reductional dysploidy (Rieger et al. 1991; Schubert and Lysák 2011; Mandáková and Lysák 2018; Lysák 2022). It is typically associated with structural rearrangements leading to fusion of two chromosomes, usually as nested fusions, less often as end-to-end fusions, whereas gain is caused by chromosome fissions and is apparently a rarer event than descending dysploidy (Stebbins 1950; Lysák et al. 2009; Escudero et al. 2014; Lysák 2014; Carta et al. 2020; Winterfeld et al. 2020a; Mayrose and Lysák 2021; Senderowicz et al. 2021; Chase et al. 2023; Xavier et al. 2024).
Dysploidy is less common than polyploidy in the studied Poodae taxa. Of the 98 grass taxa with known chromosome number, 39 (approximately 40%) are polyploid and 59 (approximately 60%) are diploid. Dysploidy occurs in only 20 cases (approximately 20%), sometimes combined with polyploidy (Anthoxanthum, Catabrosa, Colpodium, Deschampsia, Parapholis) (Table 2; Online Resource 1). In addition, many dysploid grass taxa are annuals, supporting the hypothesis “that a low chromosome number has a selective advantage in a cross-fertilized annual species, since it increases the amount of linkage and therefore the degree of constancy of a population” (Stebbins 1950: p. 458).
Increase of 1Cx and MC. Several dysploid taxa of the grasses studied show an increase in their monoploid genome size (1Cx value) and an increase in the mean chromosome DNA content (MC) compared to their euploid relatives (for the calculations see Table 13 in Online Resource 2). The most striking example is Phalaris canariensis (x = 6), whose 1Cx and MC are increased by 363% and 214%, respectively, compared to its euploid congeners (x = 7) (Table 4; Table 13 in Online Resource 2).
Within the genus Colpodium, the diploids C. biebersteinianum and C. versicolor have 2n = 2x = 4 chromosomes. This is one of the few cases of such low chromosome number in angiosperms. Their chromosomes are 310–430% larger than those of related euploid species with x = 7 from the same subtribe Coleanthinae. The size of the monoploid genomes (1Cx values) shows less striking differences. Colpodium biebersteinianum has an increase by 23%, while C. versicolor has a decrease by 10%. This shows that chromosome fusions were the main mechanism of chromosome enlargement in these x = 2 taxa.
Echinaria capitata (2n = 2x = 18) is distinguished by a 95% increase in the 1Cx value and a 35% increase in MC compared to the other diploid but euploid species within the subtribe Sesleriinae, which have a chromosome number of x = 7. This is a notable enlargement, although less pronounced than that observed in Phalaris canariensis. It appears to be associated with a change in chromosome number from x = 7 to x = 9, which may involve chromosome fissions, but this has not been investigated. It is the only example of ascending dysploidy found among the Poodae species studied.
Decrease of 1Cx and increase of MC. A moderate decrease in monoploid genome size (1Cx value) and a moderate increase in MC of about 13–23% each were found in Briza minor (x = 5) compared to its congener B. media (x = 7) and Catabrosa aquatica (x = 5) compared to the its related taxa with x = 7 in the subtribe Coleanthinae.
Decrease of 1Cx and MC and ambiguous cases. The hypotetraploid Rostraria cristata (2n = 4x-2 = 26) shows a decrease of the 1Cx value and the MC by 31–36% compared to its diploid euploid congener R. hispida (2n = 2x = 14). However, it is not certain that this is actually caused by dysploidy, as this species is allopolyploid (Tkach, unpublished data) and may possibly also contain a smaller chromosome set than R. hispida in its tetraploid chromosome complement, causing its lowered (mean) 1Cx value and MC.
Trisetum flavescens (2n = 4x = 26) stands out for its comparatively small 1Cx and MC values, which are the smallest within the comparatively extensively sampled subtribe Aveninae (Table 2). Both values are 43–59% lower than the mean values of the diploids and tetraploids of Aveninae used for comparison (Table 4). The values would be even lower, if Tricholemma jahandiezii, with its exceptionally high and rather atypical 1Cx value and MC (chromosome size), had not been excluded from the calculation of the mean values of the tetraploid Aveninae (see Appendix to Table 13 in Online Resource 2). It would be interesting to include other cytotypes of the cytogenetically rather heterogeneous Trisetum flavescens, which also comprises hexaploids (2n = 36) (Winterfeld 2006; Winterfeld and Röser 2007b) and euploid tetraploid populations with 2n = 28 based on x = 7 (CCDB). This would also provide a more reliable means of comparing whether genome size reduction in T. flavescens actually occurred in association with dysploidy than the average genome size values from other Aveninae genera used here, or whether it occurred independently of dysploidy.
The genus Anthoxanthum, including euploid taxa with x = 7 (former Hierochloe) such as A. nitens and dysploid taxa with x = 5 such as the A. aristatum, both of which are diploids (Tables 2, 3). This appears to show a decrease in 1Cx from 4.73 pg to 3.92 pg, a 17% reduction in genome size, and an increase in chromosome size from 0.68 pg to 0.78 pg assossiated with this dysploid chromosome change. The values recorded for other diploid x = 5 taxa (A. alpinum, A. maderense) seem to show the same or even stronger decrease of the monoploid genome size (1Cx value) and in- or decrease of MCs (Chumová et al. 2015). However, the situation in Anthoxanthum is rather ambiguous, as one of the x = 5 diploids, A. gracile, has a very large genome of 9.19 pg/1Cx and an MC of 1.84 pg (Chumová et al. 2015), implying a dramatic increase. Due to the contradictory data on genome up- and down-scaling, it must be assumed that very different processes of genome evolution with very different outcomes overlap in Anthoxanthum, a conclusion supported by cytogenetic and molecular analyses (Chumová et al. 2017, 2021) (see also the following chapter).
In the genus Deschampsia, there is no clear difference between D. setacea, the only euploid and diploid species (2.63 pg/1Cx and MC of 0.38 pg), and the dysploid and at the same time polyploid taxa (2n = 4x-2 or 8x-4 and 2.16–2.81 pg/1Cx, MC of 0.33–0.43 pg). All of them are perennials, whereas D. danthonioides, the only annual species of this genus studied, has a much lower 1Cx and MC (1.90 pg/1Cx; MC 0.29 pg). As D. danthonioides also has 2n = 4x-2, its small genome size values represent rather a change related to the life form (see below).
In summary, the enlargement of genomes (1Cx values) and chromosomes (MC) associated with a dysploid change in chromosome number, as found in about half of the cases of dysploidy among the grasses studied, seems to be the most common pattern observed in angiosperms in general (Lysák et al. 2006; Cheng et al. 2013; Vaio et al. 2013; Winterfeld et al. 2018, 2020b). Such an increase in genome size has been shown to be caused by massive transposon amplification due to loss of repression and elimination of transposable elements as a result of chromosome rearrangements (Theuri et al. 2005; Lou et al. 2012; Pellicer et al. 2014; Yang et al. 2014; Rockinger et al. 2016; Ferraz et al. 2023).
Genome size reduction in the course of dysploid chromosome number change, as apparently observed in some of the studied taxa (Rostraria, Trisetum), is a pattern that has rarely been observed in other angiosperms (Chase et al. 2023).
Nevertheless, the current results seem to indicate that there is no simple causal relationship between in chromosome number changes and genome size variation in dysploidy.
Different life forms and genome size
Plants that are able to complete their life cycle in one year or less, i.e. annuals (therophytes), are typically adapted to a highly seasonal climate where there is usually a period of drought, which requires the production of seeds that can survive this hostile season. Therefore, the Mediterranean region and the Middle East, with its winter rainfall regime, is one of the most important centers for the development of annual plants, including grasses. In summary, annuals should be favored when adult mortality is high and seed persistence and seedling survival are relatively high (Hjertaas et al. 2023; Poppenwimer et al. 2023).
Ploidy and genome sizes of perennials and annuals. As a characteristic of the studied annual grasses it can be noted that they are mostly diploid and only more rarely polyploid with ploidy levels of mostly 4x and only very rarely 6x, the latter occurring among the studied taxa only in the presumably hypohexaploid Parapholis incurva (2n = 6x-4 = 38). Hexaploid annuals also are otherwise quite rare in the Poodae, but are known, for example, from a few wild annual species of Avena (oat) and of Aegilops from the Triticodae (wheat and relatives), the sister supertribe of the Poodae, while hexaploids are quite frequent among the cultivated annuals (crop species) of the tribes Aveneae and Triticeae (Rajhathy and Thomas 1974; Baum 1977; Yan et al. 2016; Feldman and Levy 2023a).
The perennial species, on the other hand, do not seem to have such restrictive constraints regarding the ploidy level, being found up to 18x in the species studied (Helictochloa pratensis) (Table 2).
2C values. The likely explanation for this striking difference between annuals and perennials seems to be that annuals have an upper limit on the 2C value, i.e. the holoploid genome size, which slows down the cell cycle by requiring the replication of a larger amount of DNA during the interphase. Although larger genome size correlates with larger cell size (nucleus/plasma ratio), which could in principle be advantageous for faster growth through fewer cell divisions, the cell cycle retardation appears to be more negative for annuals, while it may be less significant for perennials (Bennett 1972, 1987; Grant 1987; Bennett and Leitch 1995, 1997; Greilhuber and Leitch 2013).
2C values of > 20 pg were found only in perennials among the species studied (Fig. 4), namely in 21 (≈ 9.0%) of all 234 species (perennials and annuals) and 31.3% of all 167 perennial species studied (excluding facultative perennials). For the lower 2C values, which account for the majority, there is little difference between perennials and annuals, and the lowest values of < 5 pg/2C are found in both.
1Cx values. The vast majority of monoploid chromosome sets had 1–5 pg/1Cx, again with apparently little difference between perennials and annuals (Fig. 4). The comparatively few large values of > 5 pg/1Cx also occur in both, as does the largest of 8–10 pg/1Cx in the tetraploid perennial Tricholemma jahandiezii, and the diploid annuals Echinaria capitata (Table 2) and Anthoxanthum gracile (Chumová et al. 2015).
All three species have special characteristics: Tricholemma jahandiezii is a taxonomically isolated endemic of the Middle Atlas in Morocco with a comparatively small range in the highlands; Anthoxanthum gracile is a phylogenetically early-diverging species of its genus (Chumová et al. 2021) with a huge genome compared to the other species. For both taxa, an ancient relict status and a small past or present population size could be assumed, which could be responsible for the spread of genome-enlarging repetitive elements in the DNA (genetic bottleneck). Although Echinaria capitata is instead a very widespread Mediterranean species, it is characterized by a dysploid chromosome set of x = 9 as opposed to x = 7, which prevails in its close relatives of the subtribe Sesleriinae. As in some other cases in the species studied, the dysploidy and the associated chromosome rearrangements probably contributed significantly to genome enlargement in the monospecific genus Echinaria (see above Phalaris).
Mean DNA content of the chromosomes. The MCs showed a similar distribution of comparable sizes of perennials and annuals as noted for the 1Cx values. Most values were between 0.1 and 0.8 pg. Only five species had chromosomes larger than 0.8 pg. These were, in addition to the three species just mentioned, the dysploid Colpodium biebersteinianum (2n = 2x = 4) and the euploid Ventenata macra with 2n = 2x = 14, which belongs to a phylogenetic lineage that also includes some other species with comparatively large genomes (Poeae subtribes Ventenatinae and Beckmanniinae) (Tables 2, 3; Online Resource 1).
Genera with perennial and annual species. Only a few of the genera studied contain both perennial and annual species.
Briza and Macrobriza. The annual B. minor (1Cx of 2.87 pg and MC of 0.57) has a smaller monoploid genome but slightly larger chromosomes than B. media (1Cx of 3.34 pg and MC of 0.47), the only perennial of this genus studied. This opposite change in both parameters is apparently related to the dysploidy in B. minor (x = 5) compared to the euploid B. media (x = 7) and the associated chromosome restructuring.
Macrobriza maxima, a diploid euploid annual (2n = 2x = 14), differs markedly from both Briza species due to its large genome of 9.32 pg/1Cx and large chromosomes (MC of 0.67 pg). However, due to its evolutionary hybrid origin (Tkach et al. 2020), it must be taxonomically excluded from Briza and therefore cannot serve as an example of genomic changes in the course of the emergence of the therophyte life form within Briza.
Deschampsia. The only annual species of Deschampsia examined, D. danthonioides (2n = 4x-2 = 26), has a 1Cx genome size and MC about 23% lower than its also dysploid but perennial congeners with (2n = 4x-2 = 26 or 8x-4 = 52) (Table 2; Online Resource 1).
Anthoxanthum. The opposite case occurs among the diploid Anthoxanthum species, where the perennial A. alpinum has the lowest 1Cx value of 2.70 pg/1Cx and the smallest chromosomes (MC of 0.55 pg), whereas the annual A. gracile has the highest 1Cx value of 9.19 pg/1Cx and the largest chromosomes (MC of 1.8 pg) (Chumová et al. 2015). This may be related to the different processes of genome evolution mentioned above that seem to occur in this genus.
Festuca and Lolium. Among the genera Festuca and Lolium as taxonomically used in study, there is no clear distinction between annuals and perennials based on genome and chromosome sizes. The smallest 1Cx values and chromosome sizes in Festuca actually occur in the annuals, namely the diploid F. octoflora (syn. Vulpia octoflora) and the tetraploid F. incurva (syn. Psilurus incurvus) (1.66 and 2.07 pg/1Cx; MC of 0.24 and 0.30 pg), but also the highest one, in the diploid F. lachenalii (syn. Micropyrum tenellum) (5.44 pg/1Cx; MC of 0.78 pg), while the numerous perennials studied are intermediate (Table 2; Online Resource 1).
The annuals of Lolium, which are all diploid (L. persicum, L. subulatum, L. temulentum), also have larger monoploid genomes and larger chromosomes (3.76–5.15 pg/1Cx; MC 0.54–0.74 pg) than the studied perennials (L. arundinaceum, L. giganteum, L. perenne), which comprise diploid and hexaploid species (2.77–3.35 pg/1Cx; MC 0.40–0.48 pg).
These results for the genera with both annual and perennial species support the conclusion drawn above from the overall distribution of monoploid genome sizes (1Cx values) across all species studied, that the two life forms do not differ per se in this genomic trait and the associated mean chromosome sizes (MCs). Thus, the transition from a perennial to an annual life form is not necessarily associated with either genome downsizing or genome inflation, as was also observed in the genus Hordeum from the related grass tribe Triticeae (Jakob et al. 2004).
Two species in our sample were facultative annuals/perennials (winter annuals), Beckmannia syzigachne (6.59 pg/2C) and Sphenopholis obtusata (5.58 pg/2C) (Finot et al. 2004; Clayton et al. 2006 onwards). Their genome sizes were largely consistent with that of their strictly perennial congeneric taxa B. eruciformis (6.21 pg/2C) and S. intermedia (5.35 pg/2C), respectively. All species were diploid. Thus, a significantly higher genome size in obligate compared to facultative perennial monocotyledons, as suggested by Bennett (1972), is not evident from these grass examples.