Next Article in Journal
In Vitro Regeneration of Miscanthus x giganteus through Indirect Organogenesis: Effect of Explant Type and Growth Regulators
Next Article in Special Issue
Morphological and Genome-Wide Evidence of Homoploid Hybridisation in Urospermum (Asteraceae)
Previous Article in Journal
Toxicity of meta-Tyrosine
Previous Article in Special Issue
Cytological Study of Cypripedium japonicum Thunb. (Orchidaceae Juss.): An Endangered Species from Korea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 12
10.3390/plants10122798
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome Size, Cytotype Diversity and Reproductive Mode Variation of Cotoneaster integerrimus (Rosaceae) from the Balkans

by
Faruk Bogunić
1,*,
Sonja Siljak-Yakovlev
2,
Irma Mahmutović-Dizdarević
3,
Alma Hajrudinović-Bogunić
1,
Mickaël Bourge
4,
Spencer C. Brown
4 and
Edina Muratović
3
1
Faculty of Forestry, University of Sarajevo, Zagrebačka 20, 71000 Sarajevo, Bosnia and Herzegovina
2
Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, CEDEX, 91405 Orsay, France
3
Laboratory for Research and Protection of Endemic Resources, Department of Biology, Faculty of Sciences, University of Sarajevo, Zmaja od Bosne 33-35, 71000 Sarajevo, Bosnia and Herzegovina
4
Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, CEDEX, 91198 Gif-sur-Yvette, France
*
Author to whom correspondence should be addressed.
Plants 2021, 10(12), 2798; https://0-doi-org.brum.beds.ac.uk/10.3390/plants10122798
Submission received: 28 November 2021 / Revised: 14 December 2021 / Accepted: 15 December 2021 / Published: 17 December 2021
(This article belongs to the Special Issue Advances in Genome Size Evolution of Plants)
Browse Figures
Versions Notes

Abstract

:
Cotoneaster integerrimus represents a multiploid and facultative apomictic system of widely distributed mountain populations. We used flow cytometry to determine genome size, ploidy level, and reproduction mode variation of the Balkan populations, supplemented by analysis of nuclear microsatellites in order to address: (i) geographic distribution and variation of cytotypes among the populations; (ii) variation of reproduction mode and the frequency of sexuality; (iii) pathways of endosperm formation among the sampled polyploids and their endosperm balance requirements; (iv) genotypic diversity and geographic distribution of clonal lineages of polyploids. The prevalence of apomictic tetraploid cytotype followed by sexual diploids and extremely rare triploids was demonstrated. This prevalence of tetraploids affected the populations’ structure composed from clonal genotypes with varying proportions. The co-occurrence of diploids and tetraploids generated higher cytotype, reproductive mode, and genotypic diversity, but mixed-ploidy sites were extremely rare. The endosperm imbalance facilitates the development and the occurrence of intermediate triploids in mixed-ploidy populations, but also different tetraploid lineages elsewhere with unbalanced endosperm. All these results showed that the South European populations of C. integerrimus have higher levels of cytotype and reproductive diversity compared to the Central European ones. Therefore, the South European populations can be considered as a potential reservoir of regional and global diversity for this species.

1. Introduction

Genome size, as a fundamental biological characteristic of evolutionary significance, has implications on overall biodiversity [1]. The causes of genome size variation in plants are varied, but the most striking is polyploidization. The evolutionary significance of polyploidy has tremendous and far-reaching implications for plant diversity [2]. It is assumed that almost 50% of plant groups have been affected by at least a single polyploidization event throughout their evolutionary history [3]. Polyploidy recomposes genome structure, alters gene expression, induces phenotypic and physiological changes, and provide adaptive potential to the polyploid plant [4,5,6]. Allopolyploidy is considered to be more frequent than autopolyploidy [7], but both processes profoundly shape plant diversity, either independently or in concert ([8], and references therein).
One of the consequences of being a polyploid is breakdown of self-incompatibility, allowing a shift towards asexual reproductive mode [9,10,11]. Apomixis (asexual seed formation) represents an effective life strategy by which apomicts preserve and maintain hybrid and heterozygous genetic lineages, cytotypes with unbalanced chromosomes, enabling their long-term persistence and dispersal [12]. Apomixis is highly correlated with polyploidy and hybridity. However, it is still unclear as to which of these factors precedes and what exact genetic and developmental mechanisms regulate the emergence of apomixis [13].
Gametophytic apomixis, the most common type of asexual reproduction, represents the formation and development of unreduced megagametophyte (parthenogenesis) coupled with the fertilization of unreduced central cell (pseudogamy) [14,15]. As a consequence of bypassing meiosis and lack of genome recombination, apomicts have a genetic structure identical to the maternal plant, thus creating clonal populations. Apomictic polyploids show a greater colonization ability by occupying more extreme ecological niches and persisting in larger distribution areas (geographic parthenogenesis, [14]) than their diploid sexual relatives [16,17,18]. Although genetically uniform, apomicts acquire genetic variation over time via accumulation of spontaneous mutations and via residual sexuality fostering diversity within clonal populations [19,20,21]. One of the mechanisms enhancing genetic diversity of apomictic populations comes from crosses between asexual and related sexual lineages, resulting in offspring that have a predominant apomictic mode of reproduction [17,22,23]. Given the long-life cycle of woody species, apomictic reproduction and independence of mates is extremely important in the early stages of the establishment of new populations and colonization of new environments.
Genus Cotoneaster Medik. (Rosaceae, Spiroaeoideae, Pyreae, [24]) includes shrubs and some small trees that are distributed in Europe, North Africa, and temperate areas of Asia (except Japan) [25]. The species number varies greatly among authors and ranges from 50–70 [26], to ≈90 [27], to ≈400 [25]. Evolutionary history and diversification of Cotoneaster have been affected by homoploid and heteroploid hybridization, as well as autopolyploidy and apomixis [25,28,29,30,31,32,33,34,35]. Both classical chromosome countings and cytometric analyses confirmed a series of ploidy levels in the genus ranging from diploids to hexaploids, with a prevalence of tetraploids.
Apomixis has been proven to be the most frequent reproductive mode in Cotoneaster that was primarily inferred from embryological data [36,37,38], breeding experiments [39], and later by molecular analyses [31,33]. Recently, apomixis has been confirmed using flow cytometry and different reproductive pathways of seed formation in tetraploids of C. integerrimus Medik. [29,30].
Key information on embryo sac development of several tri- and tetraploid Cotoneaster species (C. rosea Edgew., C. nitens Rehd. & Wils., C. bullata Bois, C. obscura Rehd. & Wils., C. acutifolia var. villosula Rehd. & Wils., C. racemiflora var. soongorica Schneid.) was provided in 1962 [36], and for C. melanocarpus Fisch. ex Blytt much later [37,38]. These species showed a strong tendency towards apomictic reproduction, namely, apospory and rarely diplospory, and the early development of their embryo sacs had general similarities with the other apomictic rosaceous genera. The most frequent development of unreduced embryo sac implied degeneration of the primary megaspore mother cell in the earliest developmental stages of the ovule and its replacement with unreduced embryo sac originating from nucellar cells (apospory). In several cases, the unreduced embryo sac developed from archesporial cells following the degeneration of the primary megaspore mother cell (diplospory). The secondary mother megaspore cell also occurred and was developed from unreduced embryo sac or via meiosis [36]. The reduced embryo sacs were rarely observed in the studied species [36,37,38]. The authors did not provide data on the polar nuclei. Conclusion of those findings is that most polyploid species of Cotoneaster are apomictic [25,32]. Sexual reproduction is less represented in Cotoneaster and follows a common pattern of development of Polygonum type of embryo sac: seeds of diploid mothers have diploid embryo and triploid endosperm; seeds of tetraploid mothers have tetraploid embryo and hexaploid endosperm [29,30].
In recent studies, only a tetraploid cytotype in C. integerrimus populations from Central Europe was found [30,40]. The same authors [30] proved that facultative apomixis was the main reproductive mode, followed by autonomous apomixis and haploid parthenogenesis, while the proportion of sexuality was 10% in the studied sample. The authors showed that different reproductive pathways involved the interaction of reduced and unreduced gametes documented in both sexuals and asexuals.
In this study, we analyzed the variation of genome size, ploidy level, and the diversity of reproductive modes of the C. integerrimus populations from the Balkans with additional samples from several other European regions. In addition, molecular analyses based on nuclear microsatellite markers were performed on the studied cytotypes. Our questions aimed to assess: (i) the geographic distribution and variation of cytotypes within populations in the sampled area; (ii) the most prevalent mode of reproduction in polyploid cytotypes and the frequency of sexuality; (iii) the pathways of endosperm formation among the sampled polyploids and their endosperm balance requirements; (iv) genotypic diversity and geographic distribution of clonal lineages of polyploids. Finally, we discuss the importance of South European populations of C. integerrimus within the context of cytotype diversity and variation in reproduction mode.

2. Results

2.1. Genome Size and Ploidy Level Variation

The flow cytometry (FCM) of 208 Cotoneaster integerrimus individuals from Bosnia and Herzegovina resulted in three distinct groups whose holoploid genome sizes (2C) corresponded to three different ploidy levels, namely, di-, tri- and tetraploid cytotypes (Table 1, Figure 1A). The two singular values corresponded to pentaploid (2C = 2.99 pg, 1Cx = 0.598 pg) and hexaploid cytotype (2C = 3.35 pg, 1Cx = 0.558 pg); thus, both values were excluded from further analyses. The 2C genome size mean values were 1.28 pg in diploids, 1.856 pg in triploids, and 2.481 pg in tetraploids. The significantly (F2, 207 = 10.98, p > 0.001) highest monoploid (1Cx) genome size value was recorded in diploid, and the lowest in triploid cytotype (Table 1).
Geographic distribution of cytotypes showed a clear prevalence of tetraploids (89% in the total sample) in each population except the Umoljani (Um) in Bosnia and Herzegovina (Table 2, Figure 1B). The sample for Figure 1B encompassed the ploidies obtained from flow cytometry as well as the estimated ploidy levels from microsatellite data for 65 individuals (explained in Table 2 and the Material and Methods section). Diploids (7.3% in total sample) were recorded in only three populations with a single individual in Borova glava (Bg) and Devečani (De), and 18 individuals in Umoljani (Um) (Table 2, Figure 1B). Triploids (3.7% in total sample) were observed in seven populations and were represented with just one individual in six populations (Vrbe-Vr, Rujište-Ru, Sovićka vrata-So, Vošac-Vo, Premužićeva staza-Pr, Rtanj-Rt) and four individuals in the population Umoljani (Table 2).

2.2. Flow Cytometric Seed Screening

Only unambiguous results of flow cytometric seed screening (FCSS) are presented. Approximately one quarter of analyzed seeds had minor additional peaks corresponding to doubled nuclear DNA levels mirroring endoreplicated embryo and/or endosperm nuclei. In total, the reproductive mode was characterized for 591 seeds from three C. integerrimus cytotypes (Table 2).
In general, inferred pathways of sexual and asexual seed formation showed great diversity and included both reduced and unreduced gametes within each cytotype (Figure 2, Table 3). The tetraploid cytotype was the most abundant, and it consequently resulted in the highest number of reproductive pathways, predominantly apomictic. Diploids involved only sexual reproduction (Table 3). The triploid cytotype was represented with only three seeds, one of sexual and two of apomictic origin (Table 3). In any case, sexual and asexual seed formation was clearly distinguished.

2.2.1. Seed Origin in Diploids

Diploids yielded exclusively sexually originated seeds having the same or increased ploidy level compared to the mother plant, which depended on sperm cell ploidies (Figure 2(S1,S2), Table 3). The profile 2x emb.:3x end. was as expected the most represented (103 seeds) and included both reduced female and male gametes. The profile 3x emb.:4x end. (eight seeds) was a result of interploid crosses between diploids and tetraploids. Those seeds originated from joint of reduced gametes (2x) from tetraploid males and a reduced egg cell (1x) of diploid mother.

2.2.2. Seed Origin in Triploids

All three analyzed triploid seeds had different profiles (Figure 2(S3,A1,A2), Table 3). The only sexual profile 3x emb.:4x end. from one seed included 2x male gamete and reduced (1x) egg cell (Table 3). The other two seeds were of apomictic origin. One included unreduced egg cell fecundated with reduced (1x) sperm cell creating 3x emb.:7x end. Profile, and the other included unreduced (3x) sperm cell creating 3x emb.:9x end. profile (Figure 2(A1,A2), Table 3).

2.2.3. Seed Origin in Tetraploids

Three pathways of sexual seed formation were observed in tetraploids (Figure 2(S4–S6), Table 3). The most frequent profile 4x emb.:6x end. (92 seeds) was formed by fertilization of a meiotically reduced egg cell (2x) with reduced 2x sperm cell (Figure 2(S4)). Fecundation of unreduced female gamete (4x) with reduced 2x male gamete resulted in the formation of Biii hybrids (6x emb.:10x end.) and was detected in 22 seeds (Figure 2(S5)). Two cases of gamete exchange between cytotypes were registered. The fecundation of a reduced egg cell (2x) of tetraploid female with reduced gamete (1x) of diploid male resulted in a 3x emb.:5x end. profile (Figure 2(S6), Table 3).
Pseudogamous seed formation had various pathways. The two most frequent apomictic profiles were 4x emb.:10x end. (148 seeds, Figure 2(A3)) and 4x emb.:12x end. (137 seeds, Figure 2(A4)). In the first case, unreduced polar nuclei were fertilized with reduced 2x male gametes (or, less likely, two 1x gametes). The second case included fecundation of unreduced polar nuclei with two 2x male gametes (Table 3). The following profiles were observed with lower frequencies (Table 3): 4x emb.:11x end. (22 seeds), 4x emb.:9x end. (11 seeds), 4x emb.:7x end. (9 seeds). The 4x emb.:11x end. profile included the fecundation of polar nuclei with 3x pollen or simultaneously with 2x and 1x sperm cells (Figure 2(A5), Table 3). The 4x emb.:9x end. (Figure 2(A6)) profile resulted from the fecundation of polar nuclei with 1x sperm cell. Exact pathway of origin for the 4x emb.:7x end. profile (Figure 2(A8)) remained unresolved, as well as several profiles with oddly high ploidies of endosperm (14x, 15x, 16x, 20x, data not shown).
Moreover, profiles that corresponded to autonomous apomixis (4x emb.:8x end.) have been documented in 32 seeds (Figure 2(A7)). Finally, two seeds originated via haploid parthenogenesis (2x emb.:6x end.) were observed. These cases involved reduced (2x) sperm cells fecundating solely central nuclei (Figure 2(A9), Table 3).

2.2.4. Endosperm Balance

All sexually derived seeds that involved the exchange of reduced gametes within diploid and tetraploid cytotypes had a balanced endosperm with a maternal to paternal genome ratio of 2m:1p (Table 3). Unbalanced endosperms with a ratio of 1m:1p and 4m:1p were observed in sexually produced seeds originating from interploid crosses between diploids and tetraploids (Table 3). On the other hand, within the sample of apomictically originated seeds (354 seeds), the apomictic profile 4x emb.:12x end. was the only one with the balanced endosperm (137 seeds, 38.7%). The remaining seeds of apomictic origin (217 seeds, 61.3%) had unbalanced endosperms (Table 3). Interestingly, the same tetraploid mother plants produced sexual and apomictic seeds having both balanced and unbalanced endosperm (Table S1).

2.3. Genotypic Variation in Cotoneaster integerrimus Populations

Five microsatellite loci were successfully amplified in the analyzed C. integerrimus populations and yielded a total of 99 alleles from 254 individuals (Table S2).
A total of 92 multilocus genotypes (MLGs) were obtained from five loci (Ng, Table 4). Each of the studied diploids (N = 13) had a unique genotype (Table S2). The clonal genotypes (N = 32) were found within a polyploid pool of 241 individuals (Table S2). Each population contained at least one clonal genotype shared by a different number of individuals within a population, indicating clonal reproduction (Ni = 194, Table 4). The proportion of the detected clones (Ng/N) ranged from 0.1 to 0.91 (Table 4), reflecting the different level of clonal structure of populations. In each population, at least two individuals shared the same genotype (Table 4). Moreover, certain populations were completely clonal (Pu–Mt. Trebević, Si-Sinjajevina, He–Austria; Ng = 1, Table 4). In the majority of populations, the effective number of genotypes was higher than one, indicating the presence of unique genotypes in populations (Table 4). Different values of genotypic diversity of populations (Table 4) showed their diverse structure of clonal and unique genotypes, suggesting different ratios of sexuality and asexuality in populations.
While the majority of clonal genotypes (N = 29) were geographically restricted within a single or two populations, the two clonal genotypes occurred at multiple sites (Table S2). The more widespread clonal genotype occurred in populations Bo–Borova glava (N = 20), Go–Mt. Igman (N = 9), Ru–Rujište (N = 9), So–Sovićka vrata (N = 2), and Um–Umoljani (N = 1), which was equivalent to 17% of all polyploids. The second clonal genotype was less widespread and was shared by plants from three populations: So–Sovićka vrata (N = 4), Um–Umoljani (N = 1), and Vr–Vrbe (N = 3). One hundred and thirty-seven (56.8%) individuals belonged to 25 spatially limited clonal genotypes, while other individuals possessed unique genotypes.
Principal coordinate analysis (PCoA) based on Jaccard distances was applied on C. integerrimus polyploid MLGs (Figure 3). PCoA revealed the pattern of grouping of MLGs dispersed along the first (PCo1) and the second (PCo2) coordinates (Figure 3). A weak regional geographic pattern of separation among MLGs groups was evident. The MLGs from western Balkans formed three groups along the PCo1: the first group consisted of MLGs from Bosnia and Herzegovina (Bg, Bo, Ru, Go, So, Um, Vr) and Croatia (Vo); the second consisted of several Bosnian MLGs (Go and Pu); and the third one was the most heterogeneous and consisted of MLGs originating from Bosnia and Herzegovina (Ca, De, Pr, Um, Vr), Croatia (Pr), and Montenegro (Si). MLGs from Italy (Mo), Austria (He), and Germany (Ro) were intermingled with Balkan MLGs. The southern Balkan MLGs from Greece (Ol and Pa) and North Macedonia (Su) clustered into a separate group along the second coordinate (Figure 3). Serbian MLGs (Mu and Rt) formed a group neighboring the southern Balkan cluster. Certain MLGs from De and Um populations were interspersed within the southern Balkan cluster.

3. Discussion

3.1. Genome Size Variation and Geographic Distribution of Cotoneaster integerrimus Cytotypes

Cotoneaster integerrimus showed stability of the monoploid genome size of three cytotypes, despite asymmetric sample size and significant differences between them. Monoploid values were quite similar to those obtained by Macková et al. [30] and Kšiňan et al. [40]. A slight decrease of monoploid genome size of polyploids compared to diploids was evident. Genome downsizing is a common phenomenon during polyploidization and inter-cytotype hybridization [42,43,44]. The obtained genome sizes and ploidies of seed embryos were consistent with leaf data. However, certain inconsistencies were present in precise determination of endosperm ploidy, which was calculated using the monoploid genome size of embryo. This issue has been often reported in studies dealing with flow cytometric seed screen [22,30]. The source of endosperm variation lies in possibility of the central cell’s fertilization by unbalanced and aneuploid gametes of variable genome size and in different number of nuclei in the central cells [22,23,45,46]. However, the observed variation in endosperm ploidy did not hinder the discrimination between the sexual and asexual origin of seeds but, in some cases, rather made it difficult to accurately infer the number and ploidy of sperm cells involved in the fertilization of central cell’s nuclei.
Cotoneaster integerrimus tetraploids are the dominant cytotype in the studied Balkan populations (Figure 1B), as also documented in Central Europe [30,40]. Populations from Central Europe were homogeneously tetraploid. Sexual diploids are detected only in the Western Alps [30]. In comparison, our study showed the greater presence of sexual diploids in the Balkans but stressed their rare occurrence as well. They were registered at the three populations in Bosnia and Herzegovina, of which two contained only one diploid individual. In contrast, it was the prevailing cytotype at the most numerous Umoljani population, outnumbering the polyploids (Figure 1B, Table 2). The observed co-occurrence of tetraploids and triploids in seven populations (Figure 1B) implicate the exchange of reduced gametes between tetraploid and diploid cytotypes. However, the sympatry of diploids, triploids, and tetraploids was documented only in the population Umoljani. We assume that diploids in other populations simply were not covered by sampling due to their overall rarity. Hence, a more extensive sampling strategy should be applied to test this assumption. If we consider the present results jointly with Macková et al. [30] and Kšiňan et al. [40], as a comprehensive dataset on C. integerrimus ploidy variation, it is apparent that diploids are a minority cytotype in the overall cytotype structure of the species. Diploids of C. integerrimus as a rare cytotype never form monoploid populations [47]. The coexistence of diploid and tetraploid cytotypes is also extremely rare in C. integerrimus. Diploid cytotype at the Umoljani site contributes markedly to ploidy richness, reproductive mode variation, and genotypic diversity. Namely, diploids affect cytotype dynamics through reproductive interactions with sexual and asexual tetraploids, generating intermediate triploid cytotype (Table 3). Furthermore, interploid crossings contribute to greater diversity of reproductive pathways, which was confirmed at the Umoljani site (Table S1). Finally, gene flow between di- and tetraploid cytotypes contributes to the higher genotypic diversity (Table 4).
The question remains as to why C. integerrimus diploids are so rare. Geographic distribution of C. integerrimus diploids in our sample confirmed a restricted and smaller geographic range relative to asexual tetraploids, which is a common pattern in different agamic groups [47,48,49,50]. These areas inhabited with diploids are considered as relict habitats in glacial refugia during the Pleistocene climate changes [10,51]. This assumes that the site Umoljani served as refugium for the persistence of diploids during glacial cycles. Such a scenario seems plausible for Umoljani as it is situated just above the huge Rakitnica gorge within the high mountains in Bosnia and Herzegovina, which has been confirmed as a refugium for many plant species [52]. Moreover, the Umoljani site has been identified as one of the hotspots of Sorbus cytotype diversity in the Balkans [23]. The current coexistence of the C. integerrimus diploids and tetraploids reflects the primary contact zone of cytotypes that involved in situ ancient formation of tetraploids via autopolyploidy or postglacial cytotype overlap of the two cytotypes (secondary contact zone). Nuclear microsatellites showed genetic divergence between diploids and tetraploids (Figure S1) and rather support the secondary contact zone. However, due to the small sample size and occurrence of clonal genotypes (six out of eight individuals) at the Umoljani site, it remains uncertain as to whether the current coexistence resulted from the primary or secondary contact. Such sites, where different sexual and asexual cytotypes coexist and interact, encouraging the formation of novel biodiversity, represent areas of particular conservation concern [53].

3.2. Asexual Seed Formation in Monoploid and Mixed-Ploidy Populations

Flow cytometric seed screening unequivocally confirmed pseudogamous apomixis as the most prevalent mode of reproduction in the C. integerrimus complex along the sampled area. While diploids are exclusively sexual, tetraploids combine asexual and sexual reproduction. Our results for tetraploids are fully consistent with those of Macková et al. [30], as well as with other related genera: Amelanchier [54], Crataegus [55], Potentilla [45], Rubus [56], and Sorbus [22,23]. All C. integerrimus tetraploids in monoploid populations showed similar patterns of seed formation involving both reduced and unreduced gametes (Table 3). Seed formation in tetraploids mainly involved central cell fertilization with one or two 2x sperm cells (two dominant profiles: 4x emb.:10x end. and 4x emb.:12x end., Table 3). The latter profile suggested dispermy of the central cell during endosperm formation [54,55,56], rather than fertilization with unreduced sperm cells. Namely, we did not notice the occurrence of male unreduced gametes in any of the supposed reproductive pathway of tetraploids (Table 3). The 10x endosperm was prevalent in Crataegus and Rubus but quite rare in Amelanchier and Sorbus, having 12x endosperm. Both 10x and 12x endosperm had almost equal frequency in C. integerrimus tetraploids, as confirmed by Macková et al. [30]. In addition, both ploidy levels of endosperm were present in seeds produced by the same plants, even yielding sexually originated seeds with balanced endosperm at the same time (Table S1). The capability of a single plant to simultaneously produce pseudogamous seeds with balanced/unbalanced endosperm, as well as seeds of sexual origin with a balanced endosperm, represents an advanced reproductive strategy. The potential to use the full benefits from both reproduction modes, accompanied with self-compatibility and self-pollination (10, 12), makes the Cotoneaster integerrimus tetraploid apomicts ecologically superior for rapid adaptation, colonization, and long-term persistence of populations in dynamic environments. The prevalence of tetraploid apomictic populations over diploid sexual ones goes in favor of geographic parthenogenesis along the sampled area (14).
The frequency of endosperm with odd ploidy levels was low (8.8%) in tetraploids of C. integerrimus, although this is common in related genera [22,23,30,53,55,56]. Such endosperm in our case would require the participation of triploid and haploid pollen that was not observed in the field, except for the population Umoljani. Lepší et al. [22] showed that the exact embryo and endosperm genome sizes primarily depend on the combination of actual gamete genome sizes involved in the fertilization and may markedly deviate from expected mean values, making profile interpretation difficult. We assume that these profiles with an unbalanced endosperm (4x emb.:7x end., 4x emb.:9x end. and 4x emb.:11x end.) represent pseudogamy that most likely included aneuploid sperm cells with variable genome size [55,56] or products of an irregular pseudogamous process [45]. However, our data showed that tetraploid cytotype apparently tolerate endosperm imbalance, but the relationship between endosperm imbalance tolerance and seed viability and germination remains uninvestigated.
In addition to pseudogamous development, we also observed a profile corresponding to autonomous apomixis (4x emb.:8x end.) (Table 3). Unlike Macková et al. [30], in our dataset, we detected some cases of G2 peaks of tetraploid embryos (Figure 2(S4)), and therefore we cannot precisely distinguish endosperm with endoreduplicated embryos. Finally, several cases of haploid parthenogenesis have been found in tetraploid seeds as a distinct rare pathway of asexual seed formation yielding unviable seeds and anormal development of the plant.

3.3. Sexuality in Cotoneaster integerrimus Populations

The fertilization via reduced (2x) gametes, resulting in seeds with a balanced endosperm (2m:1p), represents the main (19.4%, N = 92, Table 3) sexual reproductive pathway of tetraploid C. integerrimus, which is noticeably higher than in Central European populations (6.3%, [30]). As mentioned earlier, different seed formation pathways (apomictic and sexual) involving different combinations of gamete number can operate in one plant. On the other hand, the rate of sexuality within a single plant/population greatly varies, as evidenced by the proportion of clonal genotypes and complete clonality in some populations (Table 4).
The frequency (N = 22, 4.2%, Table 3) of unreduced female gametes (4x) in crosses with reduced male gametes (2x) resulted in Biii hybrid seeds, with a similar proportion found in Central European populations. Our screening revealed only one fruitless hexaploid individual in the field, and we assume that this cytotype is extremely rare. The viability of Biii hybrid seeds and seedlings is most likely low due to endosperm imbalance [22].
Sexuality in C. integerrimus also included bidirectional intercytotype crosses, resulting in 12 seeds with triploid embryos that accounted for ≈9% of all analyzed seeds from the Umoljani population. We expected that the number of such crosses would be much higher in a larger sample. Singular occurrence of triploid individuals in other populations (So, Ru, and Vr) also requires more extensive sampling to depict cytotype structure in populations. Our results provide evidence that gene flow occurs between diploid and tetraploid cytotypes in their sympatry, but the small number of triploids that have a weak or no fruit yield suggest a lower fitness than either parent (only three fertile pyrenes from two individuals). The occurrence of triploid cytotype in mixed-ploidy populations indicates the existing weak reproductive barriers between diploid and tetraploid cytotype, which nevertheless allow the formation of triploids but with a low frequency.
The key reproductive isolation mechanism in the formation of triploids is considered to be the triploid block, which operates in different polyploid systems [22,57,58,59,60]. The triploid block is caused by malfunction of the endosperm as a result of an imbalance between paternally and maternally imprinted genes during the endosperm development, and often leads to seed abortion [61]. The observed endosperm imbalance in all triploid seeds originated from di-, tri-, and tetraploid mother plants in the present dataset (Table 3), suggesting that a small portion of seeds are insensitive to the triploid block. These seeds, originated via intercytotype crosses in C. integerrimus, represent novel polyploid forms of genomic variation with potential for divergence. The frequency of repeated crosses between diploids and tetraploids accompanied by the relaxed endosperm balance would directly affect the frequency of triploids in mixed-ploidy populations. Although triploids are considered unstable and less fertile relative to diploids and tetraploids in Rosaceae [55,62], they play an important role as intermediates in the formation of tetraploids via ‘triploid bridge’ [23,59,63]. Our scenario involves bidirectional crosses between a sexual diploid cytotype and a predominantly apomictic tetraploid cytotype that produces a mostly apomictic triploid progeny of C. integerrimus. The reproductive success and persistence of the triploid cytotypes will primarily depend on the capabilities provided by their mating system [10]. Certainly, a future sampling design will encompass a larger geographic distribution of C. integerrimus populations and a larger number of individual plants/progenies to answer the question on the occurrence, dynamics, and distribution of triploids and their potential role in evolution of this complex.

3.4. Patterns of Genotypic Diversity in Cotoneaster integerrimus

The use of nuclear microsatellite markers confirmed cytometric results and shed a new light on the breeding system of C. integerrimus. The patterns of genotypic diversity within and among populations primarily reflected on the breeding system of C. integerrimus that favors pseudogamous reproduction. Favored asexuality is evident in the low level of genotypic diversity and the high proportion of clonal genotypes in most populations. While the studied diploids have unique MLGs, a trait of exclusive out-crossers [63,64], polyploids share the same genotypes in each population. Heteroploid populations showed higher levels of genotypic diversity. These patterns may reflect ancient hybridization between cytotypes that has taken place over time or current processes in which younger apomictic lineages tend to prevail via intercytotypic crosses [18,65]. The proportion of sexuality of tetraploids significantly differs among the populations, but the factor that governs the sexual reproduction in C. integerrimus remains unclear. On the other hand, the sexual reproduction of facultative apomicts may depend on changes in environmental conditions [56]. The interaction between the breeding system and environmental conditions allows for the facultative apomict to increase the rate of sexuality at a given time, reshuffling the overall genetic variability of populations in dynamic environments [66]. In any case, sexuality is an important and propulsive agent of genetic variability in monoploid and mixed-ploidy C. integerrimus populations. Additional factors contributing to genetic variability in different apomictic systems are mutations and residual recombination of the female gametophyte [20,63,67,68], but these aspects are beyond the scope of this study.
Although the most clonal MLGs were spatially restricted within the populations, few were dispersed among the populations. Moreover, the most widespread clonal MLG constituted the complete structure of the Bg population, and significantly in the Go and Ru sites. Geographical distances between the sites/populations with the clonal MLGs ranged between 40 and 70 km. The wide geographic occurrence of a certain clonal lineage confirmed its high abilities to colonize and spread in new geographic areas. This pattern might fit the idea of generalist clonal genotypes that are geographically widespread and occur in different habitats [18]. The widely distributed clonal MLGs may represent old C. integerrimus lineages that had the chance for multiple long-distance dispersal events, and spatially restricted ones may represent relatively new clonal genotypes that have not had time to disperse across environments [65]. In that context, the lack of clear geographic structure among the polyploid MLGs, as revealed by PCoA, might lie in their multiple origins via recurrent mating between widely distributed lineages and local genotypes, as well as multiple colonization of sites by different clones [18,68]. In such a scenario, recurrent admixture of divergent genotypes and their interaction will blur the geographic structuring of populations.

4. Materials and Methods

4.1. Plant Material

We sampled a total of 18 natural populations of C. integerrimus originating from the Balkans (covering Bosnia and Herzegovina and the neighboring regions of Serbia, Montenegro, Croatia, North Macedonia, and Greece (Figure 1A, Table 2)). Our sample was complemented with three populations from Italy, Austria, and Germany. More detailed sampling was performed in 10 Bosnia and Herzegovina populations encompassing 17–30 individuals per site, depending on the population size (Table 2). The distance between adult plants at all sites was at least 20 m to avoid collecting ramets. Altogether, 254 adult individuals from 21 populations were used for genome size, ploidy level determination, and molecular analysis. To determine the reproductive mode variation, we studied 591 seeds from the 63 mother individuals. The vouchers were deposited in the Herbarium of the National Museum of Bosnia and Herzegovina (SARA herbarium, voucher numbers 51809-51818).

4.2. Genome Size and Ploidy Level Determination

Absolute genome size and ploidy level determination was performed by flow cytometry (FCM) for 208 individuals from Bosnia and Herzegovina. Genome size determination followed the protocol previously used for Sorbus spp. [23]. Briefly, parts (around 1 cm2) of C. integerrimus fresh leaves and leaf material of internal standard (Solanum lycopersicum cv. Montfavet 63-5, 2C = 1.99 pg, [69]) were chopped together using a razor blade in cold Gif Nuclear Buffer [70]. The suspension was filtered through a 50 µm nylon mesh (CellTrics from Partec-Sysmex, Goertliz, Germany), and RNAse (Roche CustomBiotech, Mannheim, Germany) was added to 2.5 U mL−1. The nuclei were stained with propidium iodide (Sigma-Aldrich, Saint-Louis, United States of America) to a final concentration 50 µg mL−1 and incubated on ice for at least five minutes. The fluorescence of 5000 nuclei was recorded for each sample using a Partec CyFlow SL3 (Partec, Münster, Germany) 532 nm laser cytometer. Fluorescence histograms were analyzed using FloMax ver. 2.8 (Partec, Münster, Germany). The absolute 2C DNA values of Cotoneaster individuals were obtained by calculation of linear relationship between the fluorescence signals of unknown sample and known internal standard. Individual DNA ploidy levels [71] were inferred following earlier chromosome counts on Cotoneaster species and compared with obtained 2C DNA values. Each individual was assigned into particular cytotype: di-, tri-, and tetraploid. The total 2C DNA value/ploidy level ratio was used to obtain monoploid genome size (1Cx, [72]) of each cytotypes. The mean values of monoploid genome size of cytotype groups were tested using one-way ANOVA followed by Tukey’s HSD test. Prior to ANOVA, homogeneity of groups’ variances using Levene’s test and data distribution using Kolmogorov–Smirnov test were checked. Analyses were conducted using PAST 3.17 [73].
For individuals sampled outside Bosnia and Herzegovina (N = 65), the ploidy level was inferred from microsatellite data. The maximum number of alleles per locus was used to determine the ploidy level for those samples. This procedure was successfully applicable in certain groups [72,74,75,76], but this criterion was reliable for tetraploids only.

4.3. Flow Cytometric Seed Screening

Flow cytometric seed screening was successfully conducted on 591 seeds following the procedure described by Hajrudinović et al. [23]. Seeds (pyrenes) were collected from attached fruits on 14 diploid, 3 triploid, and 63 tetraploid mother individuals, previously cytotyped (Table S1). Only well-formed seeds, dried for 48 h at room temperature and kept in paper bags at 4 °C prior to analysis, were used. Each seed was analyzed separately. Only three healthy seeds were recovered from triploid plants, which indicates a high level of sterility of this cytotype. Entire seeds were co-chopped with part of fresh leaf (around 0.5 cm2) of internal standard Oryza sativa ssp. japonica cv. Nipponbare (2C = 0.9 pg, [77]) in cold Gif Nuclear Buffer following the steps in procedure described for genome size and ploidy level determination. Both linear and log histograms were recorded after the fluorescence was divided between two photomultipliers with a 50/50 mirror. Endosperm ploidy was calculated using the inferred monoploid genome size of embryo of the same seed. Estimated DNA ploidies of embryo and endosperm were compared to distinguish between sexual and apomictic origin of each analyzed seed and to deduce fertilization pathways according to [23,55,78].

4.4. Analysis of Nuclear Microsatellites

A modified CTAB-procedure [79] was used to extract total genomic DNA from around 20 mg of silica-dried leaf material. Amplification of five nuclear microsatellite markers was performed for 254 individuals (Table 2) following Robertson et al. [64]. The applied primers are primarily designed for Malus × domestica (CH01H10, CHO1FO2, and CHO2D11, [80]) and S. torminalis (MSS5 and MSS16, [81]) but were successfully applied in different rosaceous genera. An ABI PRISM 3500 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA) was used for electrophoretic separation of PCR products. Alleles were sized relative to the internal size standard TAMRA 500 (Applied Biosystems, Warrington, United Kingdom). Electropherograms were analyzed using GeneMapper v. 5 (Applied Biosystems).
To study genetic diversity in polyploid C. integerrimus populations, we determined the multilocus genetic genotypes (MLGs or simply genotypes) for each individual on the basis of microsatellite allele data for each of five loci in the software program GenoType v. 1.2 [82]. Assignment of individuals into a particular clone was performed using the Meirmans and Van Tienderen [82] algorithm according to the calculation of a genetic distance matrix and a clonal threshold (set to two for the study, after testing different threshold values upon recommendations by the authors of the programs) under the stepwise mutation model option. In order to analyze clonal diversity of C. integerrimus, we included following measures: the total number of genotypes (Ng) and effective genotypes (Eff), the number of unique genotypes (Un), and Simpson’s diversity index (D) calculated in the software GenoDive v. 1.2 [82]. The proportion of distinguishable genotypes (Ng⁄ N) [83] was calculated as well. Relationships among multilocus genotypes were visualized by principal coordinate analysis (PCoA) based on Jaccard distances using PAST 3.17 [73].
The maximum number of alleles per loci corresponded to ploidy levels obtained using cytometry. For 65 individuals with unknown genome size, these data were used to infer the ploidy level (Table 2 and Table S2).

5. Conclusions

From genome size, ploidy level, mode of reproduction, and genotypic variation, we demonstrated that South European C. integerrimus populations from the Balkans represent a dynamic polyploid aggregate with complex breeding. The prevalence of apomictic tetraploids determines the overall cytotype structure of populations, confirming the pattern of geographic parthenogenesis along the sampled area. Genotypic variability of populations mainly results from the interaction of predominant facultative apomixis and residual sexual reproduction in tetraploids on a large geographical scale. Sexual diploids, due to their rarity, also contribute to the cytotype, reproductive diversity, and genotypic diversity of polyploid populations but at the local scale. The endosperm imbalance facilitates the development and the occurrence of intermediate triploids in mixed-ploidy populations, but also different tetraploid lineages elsewhere with unbalanced endosperm. South European populations of C. integerrimus showed higher levels of cytotype and reproductive diversity compared to the Central European ones, representing a potential reservoir of regional and global diversity for the species. Future efforts should be focused on more extensive sampling strategy of C. integerrimus populations to reveal the underlying processes shaping the diversity of this multiploid species.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/plants10122798/s1, Table S1: Individual FCSS profiles. Table S2: Genome composition of the studied Cotoneaster integerrimus individuals at five microsatellite loci. Figure S1: Principal coordinate analysis of multilocus genotypes from the Umoljani site.

Author Contributions

Conceptualization, F.B. and E.M.; methodology, F.B., S.S.-Y., E.M. and A.H.-B.; validation, M.B. and S.S.-Y.; formal analysis, F.B. and E.M.; investigation, I.M.-D. and A.H.-B.; resources, F.B. and S.S.-Y.; writing—original draft preparation, F.B.; writing—review and editing, E.M., A.H.-B., S.S.-Y., I.M.-D., M.B., S.C.B. and F.B.; visualization, F.B. and A.H.-B.; supervision, S.S.-Y.; funding acquisition, F.B. All authors have read and agreed to the published version of the manuscript.

Funding

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are greatly thankful to Elma Silajdžić for her help with molecular laboratory work. Božo Frajman supported us with samples from Austria, Germany, and Italy; Nevena Kuzmanović with samples from Serbia and Greece; Danka Petrović with samples from Montenegro. We thank Dubravka Šoljan for her help during collection of samples and Mario Gomez Pacheco for help with cytometry. We thank three anonymous referees for valuable comments on the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pellicer, J.; Hidalgo, O.; Dodsworth, S.; Leitch, I.J. Genome size diversity and its impact on the evolution of land plants. Genes 2018, 9, 88. [Google Scholar] [CrossRef] [Green Version]
  2. Van de Peer, Y.; Mizrachi, E.; Marchal, K. The evolutionary significance of polyploidy. Nat. Rev. Genet. 2017, 18, 411–424. [Google Scholar] [CrossRef]
  3. Soltis, D.E.; Visger, C.J.; Soltis, P.S. The polyploidy revolution then and now: Stebbins revisited. Am. J. Bot. 2014, 101, 1057–1078. [Google Scholar] [CrossRef] [Green Version]
  4. Ďurkovič, J.; Kardošová, M.; Čaňová, I.; Lagaňa, R.; Priwitzer, T.; Chorvát, D.; Cicák, A.; Pichler, V. Leaf traits in parental and hybrid species of Sorbus (Rosaceae). Am. J. Bot. 2012, 99, 1489–1500. [Google Scholar] [CrossRef]
  5. Soltis, D.E.; Albert, V.A.; Leebens-Mack, J.; Bell, C.D.; Paterson, A.H.; Zheng, C.; Sankoff, D.; Depamphilis, C.W.; Wall, P.K.; Soltis, P.S. Polyploidy and angiosperm diversification. Am. J. Bot. 2009, 96, 336–348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hegarty, M.J.; Hiscock, S.J. Genomic clues to the evolutionary success of polyploid plants. Curr. Biol. 2008, 18, 435–444. [Google Scholar] [CrossRef] [Green Version]
  7. Parisod, C.; Holderegger, R.; Brochmann, C. Evolutionary consequences of autopolyploidy. New Phytol. 2010, 186, 5–17. [Google Scholar] [CrossRef] [PubMed]
  8. Barker, M.S.; Arrigo, N.; Baniaga, A.E.; Li, Z.; Levin, D.A. On the relative abundance of autopolyploids and allopolyploids. New Phytol. 2016, 210, 391–398. [Google Scholar] [CrossRef] [PubMed]
  9. Dobeš, C.; Scheffknecht, S.; Fenko, Y.; Prohaska, D.; Sykora, C.; Hülber, K. Asymmetric reproductive interference: The consequences of cross-pollination on reproductive success in sexual–apomictic populations of Potentilla puberula (Rosaceae). Ecol. Evol. 2019, 8, 365–381. [Google Scholar] [CrossRef] [Green Version]
  10. Hörandl, E. The evolution of self-fertility in apomictic plants. Sex. Plant Reprod. 2010, 23, 78–86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Barrett, S.C.H. The evolution, maintenance, and loss of self-incompatibility systems. In Plant Reproductive Ecology, 1st ed.; Lovett, J.D., Lovett, L.D., Eds.; Oxford University Press: New York, NY, USA, 1988; pp. 98–124. [Google Scholar]
  12. Hörandl, E. The complex causality of geographical parthenogenesis. New Phytol. 2006, 171, 525–538. [Google Scholar] [CrossRef] [PubMed]
  13. Hojsgaard, D.; Hörandl, E. The rise of apomixis in natural plant populations. Front. Plant Sci. 2019, 10, 358. [Google Scholar] [CrossRef]
  14. Hörandl, E.; Cosendai, A.-C.; Temsch, E.M. Understanding the geographic distributions of apomictic plants: A case for a pluralistic approach. Plant Ecol. Divers. 2008, 1, 309–320. [Google Scholar] [CrossRef]
  15. Ozias-Akins, P.; van Dijk, P.J. Mendelian genetics of apomixis in plants. Ann. Rev. Genet. 2007, 41, 509–537. [Google Scholar] [CrossRef] [PubMed]
  16. Mráz, P.; Mrázová, V. Greater reproductive assurance of asexual plant compared to sexual relative in a low density sympatric population–experimental evidence for pollen limitation. J. Evol. Biol. 2021, 34, 1503–1509. [Google Scholar] [CrossRef] [PubMed]
  17. Alonso-Marcos, H.; Nardi, F.D.; Scheffknecht, S.; Tribsch, A.; Hülber, K.; Dobeš, C. Difference in reproductive mode rather than ploidy explains niche differentiation in sympatric sexual and apomictic populations of Potentilla puberula. Ecol. Evol. 2019, 9, 3588–3598. [Google Scholar] [CrossRef] [Green Version]
  18. Coughlan, J.M.; Han, S.; Stefanović, S.; Dickinson, T.A. Widespread generalist clones are associated with range and niche expansion in allopolyploids of Pacific Northwest Hawthorns (Crataegus L.). Mol. Ecol. 2017, 26, 5484–5499. [Google Scholar] [CrossRef]
  19. Hojsgaard, D.; Hörandl, E. A little bit of sex matters for genome evolution in asexual plants. Front. Plant Sci. 2015, 6, 82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Majeský, L.; Vašut, R.J.; Kitner, M.; Trávníček, B. The pattern of genetic variability in apomictic clones of Taraxacum officinale indicates the alternation of asexual and sexual histories of apomicts. PLoS ONE 2012, 7, e41868. [Google Scholar] [CrossRef] [Green Version]
  21. Paun, O.; Stuessy, T.F.; Hörandl, E. The role of hybridization, polyploidization and glaciation in the origin and evolution of Ranunculus cassubicus complex. New Phytol. 2006, 171, 223–236. [Google Scholar] [CrossRef]
  22. Lepší, M.; Koutecký, P.; Nosková, J.; Lepší, P.; Urfus, T.; Rich, T.C.G. Versatility of reproductive modes and ploidy level interactions in Sorbus s.l. (Malinae, Rosaceae). Bot. J. Linn. Soc. 2019, 191, 502–522. [Google Scholar] [CrossRef]
  23. Hajrudinović, A.; Siljak-Yakovlev, S.; Brown, S.C.; Pustahija, F.; Bourge, M.; Ballian, D.; Bogunić, F. When sexual meets apomict–genome size, ploidy level and variation of reproductive mode in Sorbus aria s.l. and S. austriaca (Rosaceae) in Bosnia and Herzegovina. Ann. Bot. 2015, 116, 301–312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Potter, D.; Eriksson, T.; Evans, R.C.; Oh, S.H.; Smedmark, J.E.E.; Morgan, D.R.; Kerr, M.; Robertson, K.R.; Arsenault, M.P.; Dickinson, T.A.; et al. Phylogeny and classification of Rosaceae. Plant Syst. Evol. 2007, 266, 5–43. [Google Scholar] [CrossRef]
  25. Fryer, J.; Hylmö, B. Cotoneasters: A Comprehensive Guide to Shrubs for Flowers, Fruit, and Foliage, 1st ed.; Timber Press: London, UK, 2009; p. 334. [Google Scholar]
  26. Dickoré, W.B.; Kasperek, G. Species of Cotoneaster (Rosaceae, Maloideae) indigenous to, naturalising or commonly cultivated in Central Europe. Willdenowia 2010, 40, 13–45. [Google Scholar] [CrossRef]
  27. Lu, L.D.; Brach, A.R. Cotoneaster. In Flora of China, 1st ed.; Wu, Z.Y., Raven, P.H., Hong, D.Y., Eds.; Science Press: Beijing, China; Missouri Botanical Garden Press: St. Louis, MO, USA, 2003; Volume 9, pp. 85–108. [Google Scholar]
  28. Meng, K.K.; Chen, S.-F.; Xu, K.-W.; Zhou, R.-C.; Li, M.-W.; Dhamala, M.K.; Liao, W.-B.; Fan, Q. Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae). Mol. Phylogenetics Evol. 2021, 158, 107083. [Google Scholar] [CrossRef]
  29. Bourge, M.; Mahmutović, I.; Hajrudinović-Bogunić, A.; Siljak-Yakovlev, S.; Silajdžić, E.; Brown, S.C.; Muratović, E.; Bogunić, F. Cytotypes and reproductive modes of Cotoneaster integerrimus (Rosaceae) from Bosnia and Herzegovina: Preliminary results. In Proceedings of the Cytometrie: 19e Congrès annuel de l’association Française de Cytométrie, Antibes, France, 18–20 November 2015; p. 285. [Google Scholar]
  30. Macková, L.; Nosková, J.; Ďurišová, Ľ.; Urfus, T. Insights into the cytotype and reproductive puzzle of Cotoneaster integerrimus in the Western Carpathians. Plant Syst. Evol. 2020, 306, 58. [Google Scholar] [CrossRef]
  31. Li, M.; Chen, S.; Zhou, R.; Fan, Q.; Li, F.; Liao, W. Molecular evidence for natural hybridization between Cotoneaster dielsianus and C. glaucophyllus. Front. Plant Sci. 2017, 8, 704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Rothleutner, J.R.; Friddle, M.W.; Contreras, R.N. Ploidy levels, relative genome sizes, and base pair composition in Cotoneaster. J. Amer. Soc. Hort. Sci. 2016, 141, 457–466. [Google Scholar] [CrossRef] [Green Version]
  33. Nybom, H.; Bartish, I.V. DNA markers and morphometry reveal multiclonal and poorly defined taxa in an apomictic Cotoneaster species complex. Taxon 2007, 56, 119–128. [Google Scholar] [CrossRef]
  34. Kroon, G.H. Polyploidy in Cotoneaster II. Acta Bot. Neerl. 1975, 24, 417–420. [Google Scholar] [CrossRef]
  35. Sax, H.J. Polyploidy and apomixis in Cotoneaster. J. Arnold Arbor. 1954, 35, 334–365. [Google Scholar] [CrossRef]
  36. Hjelmquist, H. Embryo sac development of Cotoneaster. Bot. Not. 1962, 115, 209–236. [Google Scholar]
  37. Mandryk, V.Y. Seed reproduction features in Cotoneaster melanocarpus Fisch. ex Blytt. (Rosaceae) from Ukrainian Carpathian flora. Ukr. Botan. Zhurn. 2006, 19, 66–72. (In Ukrainian) [Google Scholar]
  38. Mandryk, V.Y. Embryologic investigation of Cotoneaster melanocarpus Fisch. ex Blytt (Rosaceae). Ukr. Botan. Zhurn. 1994, 51, 86–93. (In Ukrainian) [Google Scholar]
  39. Flinck, K.E.; Hylmö, B. Cotoneaster harrysmithii, a new species from western China. Bot. Not. 1962, 112, 17–64. [Google Scholar]
  40. Kšiňan, S.; Ďurišová, L.; Elias, P. Genome size estimation of Cotoneaster species (Rosaceae) from the Western Carpathians. Biologia 2021, 76, 1–12. [Google Scholar] [CrossRef]
  41. Doležel, J.; Bartoš, J.; Voglmayr, H.; Greilhuber, J. Nuclear DNA content and genome size of trout and human. Cytometry 2003, 51, 127–128. [Google Scholar] [CrossRef]
  42. Leitch, I.J.; Bennett, M.D. Genome downsizing in polyploid plants. Biol. J. Linn. Soc. 2004, 82, 651–663. [Google Scholar] [CrossRef] [Green Version]
  43. Liu, H.-M.; Ekrt, L.; Koutecky, P.; Pellicer, J.; Hidalgo, O.; Marquardt, J.; Pustahija, F.; Ebihara, A.; Siljak-Yakovlev, S.; Gibby, M.; et al. Polyploidy does not control all: Lineage-specific average chromosome length constrains genome size evolution in ferns. J. Syst. Evol. 2019, 57, 418–430. [Google Scholar] [CrossRef] [Green Version]
  44. Zenil-Ferguson, R.; Ponciano, J.M.; Burleigh, J.G. Evaluating the role of genome downsizing and size thresholds from genome size distributions in angiosperms. Am. J. Bot. 2016, 103, 1175–1186. [Google Scholar] [CrossRef] [Green Version]
  45. Dobeš, C.; Lückl, A.; Hülber, K.; Paule, J. Prospects and limits of the flow cytometric seed screen–insights from Potentilla sensu lato (Potentilleae, Rosaceae). New Phytol. 2013, 198, 605–616. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Talent, N.; Dickinson, T.A. The potential for ploidy level increases and decreases in Crataegus (Rosaceae, Spiraeoideae, tribe Pyreae). Can. J. Bot. 2007, 85, 570–584. [Google Scholar] [CrossRef] [Green Version]
  47. Kolář, F.; Čertner, M.; Suda, J.; Schönswetter, P.; Husband, B.C. Mixed-ploidy species: Progress and opportunities in polyploid research. Trends Plant Sci. 2017, 22, 1041–1055. [Google Scholar] [CrossRef] [PubMed]
  48. Goad, D.M.; Baxter, I.; Kellog, E.A.; Olsen, K.M. Hybridization, polyploidy and clonality influence geographic patterns of diversity and salt tolerance in the model halophyte seashore paspalum (Paspalum vaginatum). Mol. Ecol. 2021, 30, 148–161. [Google Scholar] [CrossRef] [PubMed]
  49. Dickinson, T.A.; Yan, B.X.; Han, S.; Zarrei, M. Niche shifts, hybridization, polyploidy and geographic parthenogenesis in western North American hawthorns (Crataegus subg. Sanguineae, Rosaceae). Agronomy 2021, 11, 2133. [Google Scholar] [CrossRef]
  50. Cosendai, A.C.; Hörandl, E. Cytotype stability, facultative apomixis and geographical parthenogenesis in Ranunculus kuepferi (Ranunculaceae). Ann. Bot. 2010, 105, 457–470. [Google Scholar] [CrossRef] [Green Version]
  51. Kolář, F.; Štech, M.; Trávníček, P.; Rauchová, J.; Urfus, T.; Vít, P.; Kubešová, M.; Suda, J. Towards resolving the Knautia arvensis agg. (Dipsacaceae) puzzle: Primary and secondary contact zones and ploidy segregation at landscape and microgeographic scales. Ann. Bot. 2009, 103, 963–974. [Google Scholar] [CrossRef] [Green Version]
  52. Ritter-Studnička, H. Reliktgesellschaften auf Dolomitböden in Bosnien und der Hercegovina. Vegetatio 1967, 15, 190–212. [Google Scholar] [CrossRef]
  53. Ennos, R.A.; Whitlock, R.; Fay, M.F.; Jones, B.; Neaves, L.E.; Payne, R.; Taylor, I.; de Vere, N.; Hollingsworth, P.M. Process-based species action plans: An approach to conserve contemporary evolutionary processes that sustain diversity in taxonomically complex groups. Bot. J. Linn. Soc. 2012, 168, 194–203. [Google Scholar] [CrossRef] [Green Version]
  54. Burgess, M.B.; Cushman, K.R.; Doucette, E.T.; Talent, N.; Frye, C.T.; Campbell, C.S. Effects of apomixis and polyploidy on diversification and geographic distribution in Amelanchier (Rosaceae). Am. J. Bot. 2014, 101, 1375–1387. [Google Scholar] [CrossRef]
  55. Talent, N.; Dickinson, T.A. Endosperm formation in aposporous Crataegus (Rosaceae, Spiraeoideae, tribe Pyreae): Parallels to Ranunculaceae and Poaceae. New Phytol. 2007, 173, 231–249. [Google Scholar] [CrossRef]
  56. Šarhanová, P.; Vašut, R.J.; Dančák, M.; Bureš, P.; Trávníček, B. New insights into the variability of reproduction modes in European populations of Rubus subgen. Rubus: How sexual are polyploid brambles? Sex. Plant Reprod. 2012, 25, 319–335. [Google Scholar] [CrossRef]
  57. Sutherland, B.L.; Galloway, L.F. Variation in heteroploid reproduction and gene flow across a polyploid complex: One size does not fit all. Ecol. Evol. 2021, 11, 9676–9688. [Google Scholar] [CrossRef]
  58. Köhler, K.; Scheid, O.M.; Erilova, A. The impact of the triploid block on the origin and evolution of polyploid plants. Trends Genet. 2010, 26, 142–148. [Google Scholar] [CrossRef] [PubMed]
  59. Husband, B.C. The role of triploid hybrids in the evolutionary dynamics of mixed-ploidy populations. Biol. J. Linn. Soc. 2004, 82, 537–546. [Google Scholar] [CrossRef]
  60. Scott, R.J.; Spielman, M.; Bailey, J.; Dickinson, H.G. Parent-of-origin effects on seed development in Arabidopsis thaliana. Development 1998, 125, 3329–3341. [Google Scholar] [CrossRef]
  61. Vinkenoog, R.; Bushell, C.; Spielman, M.; Adams, S.; Dickinson, H.G.; Scott, R.J. Genomic imprinting and endosperm development in flowering plants. Mol. Biotechnol. 2003, 25, 149–184. [Google Scholar] [CrossRef]
  62. Ludwig, S.; Robertson, A.; Rich, T.C.G.; Djordjevic, M.; Cerovic, R.; Houston, L.; Harris, S.A.; Hiscock, S.J. Breeding systems, hybridization and continuing evolution in Avon Gorge. Sorbus. Ann. Bot. 2013, 111, 563–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Hamston, T.J.; de Vere, N.; King, R.A.; Pellicer, J.; Fay, M.F.; Cresswell, J.E.; Stevens, J.R. Apomixis and hybridization drives reticulate evolution and phyletic differentiation in Sorbus L.: Implications for conservation. Front. Plant Sci. 2018, 9, 1796. [Google Scholar] [CrossRef]
  64. Robertson, A.; Rich, T.C.G.; Allen, A.M.; Houston, L.; Roberts, C.; Bridle, J.R.; Harris, S.A.; Hiscock, S.J. Hybridization and polyploidy as drivers of continuing evolution and speciation in Sorbus. Mol. Ecol. 2010, 19, 1675–1690. [Google Scholar] [CrossRef]
  65. Šarhanová, P.; Sharbel, T.F.; Sochor, M.; Vašut, R.J.; Dančáak, M.; Trávníček, B. Hybridization drives evolution of apomicts in Rubus subgenus Rubus: Evidence from microsatellite markers. Ann. Bot. 2017, 120, 317–328. [Google Scholar] [CrossRef]
  66. Karunarathne, P.; Reutemann, A.V.; Schedler, M.; Glücksberg, M.A.; Martínez, E.J.; Honfi, A.I.; Hojsgaard, D.H. Sexual modulation in a polyploid grass: A reproductive contest between environmentally inducible sexual and genetically dominant apomictic pathways. Sci. Rep. 2020, 10, 8319. [Google Scholar] [CrossRef]
  67. Paule, J.; Sharbel, T.F.; Dobeš, C. Apomictic and sexual lineages of the Potentilla argentea L. group (Rosaceae): Cytotype and molecular genetic differentiation. Taxon 2011, 60, 721–732. [Google Scholar] [CrossRef]
  68. Paun, O.; Greilhuber, J.; Temsch, E.M.; Hörandl, E. Patterns, sources and ecological implications of clonal diversity in apomictic Ranunculus carpaticola (Ranunculus auricomus complex, Ranunculaceae). Mol. Ecol. 2006, 15, 897–910. [Google Scholar] [CrossRef]
  69. Lepers-Andrzejewski, S.; Siljak-Yakovlev, S.; Brown, S.C.; Wong, M.; Dron, M. Diversity and dynamics of plant genome size: An example of polysomaty from a cytogenetic study of Tahitian vanilla (Vanilla × tahitensis, Orchidaceae). Am. J. Bot. 2011, 98, 986–997. [Google Scholar] [CrossRef] [PubMed]
  70. Bourge, M.; Brown, S.C.; Siljak-Yakovlev, S. Flow cytometry as tool in plant sciences, with emphasis on genome size and ploidy level assessment. Genet. Appl. 2018, 2, 1–12. [Google Scholar] [CrossRef] [Green Version]
  71. Suda, J.; Krahulcová, A.; Trávníček, P.; Krahulec, F. Ploidy level versus DNA ploidy level: An appeal for consistent terminology. Taxon 2006, 55, 447–450. [Google Scholar] [CrossRef]
  72. Greilhuber, J.; Doležel, J.; Lýsak, M.; Bennett, M.D. The origin, evolution and proposed stabilization of the terms ‘genome size’ and ‘C-value’ to describe nuclear DNA contents. Ann. Bot. 2005, 95, 255–260. [Google Scholar] [CrossRef]
  73. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 1–9. [Google Scholar]
  74. Feulner, M.; Weig, A.; Voss, T.; Schott, L.F.; Aas, G. Central European polyploids of Sorbus subgenus Aria (Rosaceae) recurrently evolved from diploids of central and south-eastern Europe: Evidence from microsatellite data. Bot. J. Linn. Soc. 2019, 191, 315–324. [Google Scholar] [CrossRef]
  75. Uhrinová, V.; Zozomová-Lihová, J.; Bernátová, D.; Paule, J.; Paule, L.; Gömöry, D. Origin and genetic differentiation of pink-flowered Sorbus hybrids in the Western Carpathians. Ann. Bot. 2017, 120, 271–284. [Google Scholar] [CrossRef]
  76. Beck, J.B.; Alexander, P.J.; Allphin, L.; Al-Shehbaz, I.A.; Rushworth, C.; Bailey, C.D.; Windham, M.D. Does hybridization drive the transition to asexuality in diploid Boechera? Evolution 2011, 66, 985–995. [Google Scholar] [CrossRef]
  77. Uozu, S.; Ikehashi, H.; Ohmido, N.; Ohtsubo, H.; Ohtsubo, E.; Fukui, K. Repetitive sequences: Cause for variation in genome size and chromosome morphology in the genus Oryza. Plant Mol. Biol. 1997, 35, 791–799. [Google Scholar] [CrossRef]
  78. Matzk, F.; Meister, A.; Schubert, I. An efficient screen for reproductive pathways using mature seeds of monocots and dicots. Plant J. 2000, 21, 97–108. [Google Scholar] [CrossRef] [PubMed]
  79. Tel-Zur, N.; Abbo, S.; Myslabodski, D.; Mizrahi, Y. Modified CTAB procedure for DNA isolation from epiphytic cacti of genera Hylocereus and Selenicereus (Cactaceae). Plant Mol. Biol. Rep. 1999, 17, 249–254. [Google Scholar] [CrossRef]
  80. Gianfranceschi, L.; Seglia, R.; Tarchini, R.; Komjanc, M.; Gessler, C. Simple sequence repeats for the genetic analysis of apple. Theor. Appl. Genet. 1999, 96, 1069–1076. [Google Scholar] [CrossRef]
  81. Oddou-Muratorio, S.; Aligon, C.; Decroocq, S.; Plomion, C. Microsatellite primers for Sorbus torminalis and related species. Mol. Ecol. Notes 2001, 1, 297–299. [Google Scholar] [CrossRef]
  82. Meirmans, P.G.; Van Tienderen, P.H. GENOTYPE and GENODIVE: Two programs for the analysis of genetic diversity of asexual organisms. Mol. Ecol. Resour. 2004, 4, 792–794. [Google Scholar] [CrossRef]
  83. Ellstrand, N.C.; Roose, M.L. Patterns of genotypic diversity in clonal plant species. Am. J. Bot. 1987, 74, 123–131. [Google Scholar] [CrossRef]
Plants 10 02798 g001 550
Figure 1. Absolute nuclear genome size (2C) of diploid, triploid, and tetraploid cytotypes (A) (boxes define 95th percentile interval; horizontal lines indicate mean value; whiskers indicate standard deviation: outliers are shown as circles), and geographical distribution of Cotoneaster integerrimus cytotypes and their frequency (%) per site/population (B). Site abbreviations (Site IDs) correspond to Table 2.
Figure 1. Absolute nuclear genome size (2C) of diploid, triploid, and tetraploid cytotypes (A) (boxes define 95th percentile interval; horizontal lines indicate mean value; whiskers indicate standard deviation: outliers are shown as circles), and geographical distribution of Cotoneaster integerrimus cytotypes and their frequency (%) per site/population (B). Site abbreviations (Site IDs) correspond to Table 2.
Plants 10 02798 g001
Plants 10 02798 g002 550
Figure 2. Flow cytometric histograms (log abscissa) of Cotoneaster integerrimus seed formation from sexual (S1–S6) or asexual (A1–A9) origin. The first fluorescence peaks correspond to the internal standard (ST, Oryza sativa L. ssp. japonica ‘Nipponbare’), the second ones to the embryo, and the third ones to the endosperm (blue: 2x mother, red: 3x mother, green: 4x mother).
Figure 2. Flow cytometric histograms (log abscissa) of Cotoneaster integerrimus seed formation from sexual (S1–S6) or asexual (A1–A9) origin. The first fluorescence peaks correspond to the internal standard (ST, Oryza sativa L. ssp. japonica ‘Nipponbare’), the second ones to the embryo, and the third ones to the endosperm (blue: 2x mother, red: 3x mother, green: 4x mother).
Plants 10 02798 g002
Plants 10 02798 g003 550
Figure 3. Principal coordinate analysis of the Jaccard distance matrix of multilocus genotypes of Cotoneaster integerrimus polyploids according to nuclear microsatellite data. Population abbreviations correspond to those in the Table 4. The numbers in the brackets denote the number of clones.
Figure 3. Principal coordinate analysis of the Jaccard distance matrix of multilocus genotypes of Cotoneaster integerrimus polyploids according to nuclear microsatellite data. Population abbreviations correspond to those in the Table 4. The numbers in the brackets denote the number of clones.
Plants 10 02798 g003
Table
Table 1. Absolute genome size values and ploidy levels of Cotoneaster integerrimus cytotypes.
Table 1. Absolute genome size values and ploidy levels of Cotoneaster integerrimus cytotypes.
DNA PloidyN2C DNACV %1Cx DNA
Mean ± SD (pg)Min–Max (pg)Mean ± SD (pg) *Mean (Mbp **)
2x201.280 ± 0.0461.23–1.362.470.636 ± 0.015 a622.00
3x71.856 ± 0.0911.63–1.905.040.608 ± 0.030 b594.62
4x1812.481 ± 0.0602.36–2.732.450.621 ± 0.015 b607.73
TOTAL208
* Mean ± SD (pg) followed by the same letter in superscript was not significantly different according to the Tukey’s HSD test; ** 1pg = 978 Mbp [41].
Table
Table 2. Geographic origin of Cotoneaster mother individuals covering three cytotypes and corresponding sample sizes (N) for genome size/ploidy level, reproductive mode, and microsatellite analysis.
Table 2. Geographic origin of Cotoneaster mother individuals covering three cytotypes and corresponding sample sizes (N) for genome size/ploidy level, reproductive mode, and microsatellite analysis.
Origin of SampleNumber of Samples (N)
Site No.Site IDLocalityEastNorthAltitude (m)DiploidsTriploidsTetraploidsTotal N per Site
1GoGornja Grkarica, Mt. Igman, Bosnia and Herzegovina (B&H)18.2883343.745281350-, -, --, -, -20, 98, 1920, 98, 19
2UmUmoljani, Mt. Bjelašnica, B&H18.2216743.66194130018, 108, 124, 3, 48, 27, 830, 138, 24
3VrVrbe, Čemerno, B&H18.5716743.178611090-, -, -1, -, -19, -, 1820, -, 18
4BoBosiljna, Mt. Čvrsnica, B&H17.4880643.511941319-, -, --, -, -20, 7, 1620, 7, 16
5RuRujište, Mt. Prenj, B&H17.9547243.461941010-, -, -1, -, 117, 14, 1818,14, 19
6BgBorova glava, Livno, B&H17.1213943.7786111831, 3, --, -, -22, 28, 2123, 31, 21
7SoSovićka vrata, Blidinje, B&H17.5047243.595281220-, -, -1, -, 116, 4, 1317, 4, 14
8PuPuhova ravan, Mt. Trebević, B&H18.4383343.828611450-, -, --, -, -20, 55, 2020, 55, 20
9CaČavljak, Mt. Ozren, B&H18.4383343.911941237-, -, --, -, -21, 75, 1921, 75, 19
10DeDevečanske stijene, Mt. Vlašić, B&H17.6213944.2788917521, 2, 1-, -, -18, 57, 1819, 59, 19
11RtMt. Rtanj, Serbia21.8891743.778611282-, -, -1 *, -, 18 *, -, 89, -, 9
12MuMt. Mučanj, Serbia20.0388943.545281530-, -, --, -, -7 *, 36, 77, 36, 7
13PrPremužićeva staza, Mt. Velebit, Croatia14.9872244.762221613-, -, -1 *, -, 13 *, 2, 34, 2, 4
14VoVošac, Mt. Biokovo, Croatia17.0547243.311941320-, -, -1 *, -, 18 *, 58, 89, 58, 9
15SiMt. Sinjajevina, Monte Negro19.3386142.961671773-, -, --, -, -5 *, -, 55, -, 5
16HeHechenberg-Ostgipfel von Kranebitten, Innsbruck, Austria11.3197247.279721390-, -, --, -, -5 *, -, 55, -, 5
17RoRotenfels, Rheinland-Pfalz, Germany7.81861149.83028232-, -, --, -, -5 *, 6, 55, 6, 5
18MoMonte Ventosola, Monte Sibillini, Umbria, Italy13.1702842.778331493-, -, --, -, -4 *, 8, 44, 8, 4
19SuMt. Suva Gora, Northern Macedonia21.2558341.844721184-, -, --, -, -7 *, -, 77, -, 7
20OlMt. Olympus, Greece22.3727840.077781935-, -, --, -, -5 *, -, 55, -, 5
21PaMt. Pangaion, Greece24.0733340.927781662-, -, --, -, -5 *, -, 55, -, 5
Total N per cytotype20, 113, 1310, 3, 9243, 475, 232273, 591, 254
* Ploidy level was inferred from microsatellite data.
Table
Table 3. Flow cytometric seed screening results and hypothesized pathways of seed formation and reproduction modes in analyzed Cotoneaster integerrimus populations.
Table 3. Flow cytometric seed screening results and hypothesized pathways of seed formation and reproduction modes in analyzed Cotoneaster integerrimus populations.
Cytometric ResultsHypothesized Pathways of Seed Formation
Maternal ploidy aEmbryopg ± SDEndosperm pg ± SDEmbryo ploidyEndosperm ploidyNumber of seedsEgg ploidyPolar nuclei ploidy/numberSperm ploidyNumber of sperms fecundating the endospermEndosperm maternal: paternal genome ratioEndosperm/embryo fluorescence (mean ratio)Seed origin
2x1.28 ± 0.031.95 ± 0.052x3x1031x1x/21x12m:1p1.52Sexual
2x1.91 ± 0.052.59 ± 0.063x4x8 b1x1x/22x11m:1p1.35Sexual
3x1.942.643x4x11x1x/22x11m:1p1.5Sexual
3x1.853.933x7x13x3x/21x16m:1p1.35Pseudogamous Apomixis
3x1.945.723x9x13x3x/23x12m:1p2.12Pseudogamous Apomixis
4x2.54 ± 0.073.84 ± 0.114x6x922x2x/22x12m:1p1.53Sexual
4x3.80 ± 0.106.43 ± 0.216x10x22 c4x4x/22x24m:1p2.5Sexual
4x1.90 ± 0.013.19 ± 0.053x5x2 d2x2x/11x14m:1p1.7Sexual
4x2.56 ± 0.16.38 ± 0.154x10x1484x4x/21x or 2x2 or 14m:1p2.5Pseudogamous Apomixis
4x2.54 ± 0.017.51 ± 0.374x12x1374x4x/22x22m:1p2.98Pseudogamous Apomixis
4x2.55 ± 0.047.04 ± 0.254x11x224x4x/23x12.7m:1p2.75Pseudogamous Apomixis
4x2.58 ± 0.125.83 ± 0.284x9x114x4x/21x18m:1p2.17Pseudogamous Apomixis
4x2.53 ± 0.055.08 ± 0.174x8x324x4x/2-04m:0p4:1Autonomous Apomixis
4x2.55 ± 0.084.46 ± 0.174x7x94x Not resolved
4x1.28 ± 0.003.84 ± 0.042x6x22x2x/21x or 2x2 or 1 Haploid Parthenogenesis
591
a Maternal ploidy inferred from cytometric leaf measurements. b Embryo ploidy was increased compared to mother ploidy. c Biii hybrids (explanation Section 2.2.3). d Embryo ploidy was decreased compared to mother ploidy.
Table
Table 4. Genetic diversity measures for the analyzed Cotoneaster integerrimus populations derived from nuclear microsatellites.
Table 4. Genetic diversity measures for the analyzed Cotoneaster integerrimus populations derived from nuclear microsatellites.
Locality (Site ID)Number of Individuals Sampled (N)Number of Multilocus Genotypes Detected (Ng)Total Number of Individuals Belonging to a Clone (Ni)Total Number of Unique Genotypes per Population (Un)Effective Number of Genotypes (Eff)Proportion of Clones Detected (Ng/N)Genotypic Diversity (Div)
Gornja Grkarica (Go)1931812.210.150.57
Umoljani (Um)2421419180.910.98
Vrbe (Vr)1871264.370.390.81
Bosiljna (Bo)1651453.550.310.76
Rujište (Ru)1941742.390.210.61
Borova glava (Bg)2122011.090.10.09
Sovićka vrata (So)1461144.070.420.81
Puhova ravan (Pu)20120110.050
Čavljak (Ca)1951641.570.260.38
Devečanske stijene (De)1961662.630.310.65
Mt. Rtanj (Rt)947430.440.75
Mt. Mučanj (Mu)72721.860.280.47
Premužićeva staza (Pr)43232.6610.83
Vošac (Vo)94732.070.330.58
Mt. Sinjajevina (Si)515110.200
Hechenberg-Ostgipfel (He)515110.200
Rotenfels (Ro)53332.270.600.70
Monte Ventosola (Mo)43232.660.750.83
Mt. Suva Gora (Su)73421.810.710.52
Mt. Olympus (Ol)54243.570.80.90
Mt. Pangaion (Pa)54243.570.80.90
2549219481
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bogunić, F.; Siljak-Yakovlev, S.; Mahmutović-Dizdarević, I.; Hajrudinović-Bogunić, A.; Bourge, M.; Brown, S.C.; Muratović, E. Genome Size, Cytotype Diversity and Reproductive Mode Variation of Cotoneaster integerrimus (Rosaceae) from the Balkans. Plants 2021, 10, 2798. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10122798

AMA Style

Bogunić F, Siljak-Yakovlev S, Mahmutović-Dizdarević I, Hajrudinović-Bogunić A, Bourge M, Brown SC, Muratović E. Genome Size, Cytotype Diversity and Reproductive Mode Variation of Cotoneaster integerrimus (Rosaceae) from the Balkans. Plants. 2021; 10(12):2798. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10122798

Chicago/Turabian Style

Bogunić, Faruk, Sonja Siljak-Yakovlev, Irma Mahmutović-Dizdarević, Alma Hajrudinović-Bogunić, Mickaël Bourge, Spencer C. Brown, and Edina Muratović. 2021. "Genome Size, Cytotype Diversity and Reproductive Mode Variation of Cotoneaster integerrimus (Rosaceae) from the Balkans" Plants 10, no. 12: 2798. https://0-doi-org.brum.beds.ac.uk/10.3390/plants10122798

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop