Translator Disclaimer
1 January 2012 The Evolution of Genome Size Variation in Drumstick Onions (Allium subgenus Melanocrommyum)
Author Affiliations +
Abstract

Allium subgenus Melanocrommyum is a species-rich group of perennial onions with uniform karyology, ecology, and breeding system. Thus, the cytological, ecological, and physiological factors often correlated with genome size should have negligible effect on genome size variation in Melanocrommyum. We measured DNA content in subgenus Melanocrommyum using flow cytometry based on propidium iodide staining and analyzed the evolution of genome size in a phylogenetic context. The observed 2C genome size variation in 160 accessions of 70 species of the subgenus was high, varying from 26.26–78.73 pg. The significant phylogenetic signal in genome size data suggests that distribution of genome size is in accordance with phylogenetic clades identified by the analysis of nuclear ITS sequences. Estimation of ancestral genome sizes using generalized least squares revealed lineages with increasing as well as decreasing DNA content. We found within-species genome size variation to be mostly below 2.5%. In species where intraspecific genome size differences were in a range of 6–9%, we suggest the existence of cryptic species, as previously inferred by molecular markers. Thus, genome size variation reflects incipient speciation or diversification in Allium subgenus Melanocrommyum. About two-fold differences in DNA content in several Melanocrommyum species indicate the occurrence of diploid and tetraploid cytotypes in these taxa, which for some species has been confirmed by chromosome counts.

Angiosperms vary tremendously in nuclear DNA content and much of this genome size variation is due to non-genic, repetitive DNA. It is well documented that in addition to polyploidization (Leitch and Bennett 1997; Soltis and Soltis 2000), transposon amplification is a major cause of genome expansion (Kalendar et al. 2000; Hawkins et al. 2008; Grover and Wendel 2010). The mechanisms responsible for genome size decrease include unequal crossing over, illegitimate recombination, a higher overall rate of deletions than insertions, and selection against transposable elements (Morgan 2001; Devos et al. 2002; Petrov 2002; Wendel et al. 2002; Ma et al. 2004; Bennetzen et al. 2005; Grover et al. 2007). The biological and evolutionary significance of wide variation in DNA content is a longstanding question in evolutionary biology. Most theories about genome size evolution emphasize the effects of DNA mass and volume on phenotype. Such effects include an increase in nucleus and cell size, prolonged duration of both mitosis and meiosis, and possible consequences of these changes such as slower cell growth and increased minimum generation time (Bennett 1972, 1987; Grime 1986; Wyman et al. 1997). While the above-mentioned theories focus on the adaptive importance of genome size variation, others point to the role of stochastic processes in genome size evolution (Lynch and Conery 2003). Oliver and colleagues (2007) propose a proportional model of evolution, implying that larger genomes evolve faster without imposing selection pressure. Petrov (2001) argues that while the mutations affecting the phenotype are subject to natural selection, the mutations with negligible selective effects are influenced mainly by genetic drift. Therefore, in the latter case the fixation of genome size will be a result of neutral stochastic events (Jakob et al. 2004).

Recent accumulation of C-value data (e.g. Plant DNA C-values database;  http://data.kew.org/cvalues/homepage.html) and the availability of phylogenies for many plant groups present an opportunity to study genome size evolution in angiosperms within a phylogenetic context. Few studies have followed this approach at the family or genus level (Bennetzen and Kellogg 1997; Cox et al. 1998; Jakob et al. 2004; Johnston et al. 2005; Leitch et al. 2007), while others have conducted broader analyses across all angiosperms (Leitch et al. 1998; Soltis et al. 2003) or the entire embryophyta (Leitch et al. 2005). However, rigorous analyses of genome size variation among closely related species in a phylogenetic framework are still missing for most plant groups. Several studies found a correlation of genome size with life history, breeding system, and ecology within several taxonomic groups (Albach and Greilhuber 2004; Weiss-Schneeweiss et al. 2006) and even within single species. Within-species genome size variation has been found in samples from geographically separated populations (Jakob et al. 2004; Schmuths et al. 2004; Smarda and Bures 2006), while in other cases variation in DNA content was correlated with ecological differences (Kalendar et al. 2000; Knight and Ackerly 2002; Jakob et al. 2004) or with differences in plant phenotype (Knight et al. 2005; Murray 2005; Beaulieu et al. 2007; Achigan-Dako et al. 2008). However, neither spatial pattern nor correlation of genome size with vegetation types or microclimatic conditions was found in a highly variable population of Festuca pallens Host (Šmarda et al. 2008).

In this study we analyze genome size evolution using phylogenetic and statistical approaches in a species-rich group of drum-stick onions (Allium L. subg. Melanocrommyum (Webb et Berthel.) Rouy). Subgenus Melanocrommyum is one of the largest in the genus Allium comprising approximately 160 perennial species with an extremely short developmental period (ephemeral growth form; Hanelt et al. 1992) and outcrossing mating system (R. M. Fritsch, pers. comm.). Melanocrommyum taxa are almost exclusively adapted to arid conditions, and possess uniform karyotypes (Fritsch and Astanova 1998). The members of the subgenus are diploid (base chromosome number x = 8) with a few exceptions. Occasionally, tetraploid individuals have been reported in species otherwise known to be diploid but they were mostly neglected due to their rare occurrence (A. atropurpureum Waldst. et Kit., Ohri et al. 1998; A. nigrum L., Shopova 1972; A. akaka Gmelin ex Roem. et Schult., Vakhtina 1969). The observed karyological and ecological uniformity (i.e. absence of adaptive factors that often are correlated with genome size) makes the group particularly interesting to study genome size variation and evolution. Thus, any variation within the subgenus is expected to be phylogenetically constrained and influenced by the number of speciation events and neutral population-level phenomena (e.g. genetic drift and population size). Studies concerning genome size variation across the genus Allium using Feulgen cytophotometry (Ohri et al. 1998; Ohri and Pistrick 2001) reported variation in genome size in Melanocrommyum even though only 15 species of the subgenus were studied. Interestingly, the authors found no general correlation between genome size and life strategy differences among species of the genus Allium (Ohri and Pistrick 2001). In the present investigation genome size variation among 70 taxa of subgenus Melanocrommyum was determined using flow cytometry, and analyzed in a phylogenetic framework using the generalized least squares (GLS) method (Pagel 1997, 1999). Here, we (1) examine the range of genome size variation in Melanocrommyum, (2) study the correlation between genome size variation and phylogenetic relationships, (3) determine if there is evidence for alterations in genome sizes (both increase and decrease) in phylogenetic groups, (4) provide evidence that diploid and tetraploid cytotypes occur within single species of subgenus Melanocrommyum, and (5) suggest that within-species genome size differences confirm the existence of cryptic species inferred by molecular markers and emphasize the taxonomic importance of genome size variation.

MATERIALS AND METHODS

Plant Material—Fresh leaf material of 160 accessions (70 taxa) representing all current taxonomic sections and phylogenetic groups of subgenus Melanocrommyum was obtained in early spring from the Allium Reference Collection of the IPK Gatersleben (Germany) (Table 1). For the majority of taxa two or more accessions were analyzed. Voucher specimens were deposited in the herbarium of the IPK Gatersleben (GAT).

Nuclear Genome Size Estimation—For flow cytometric measurements of the genome size, samples of fresh young leaves were chopped with a sharp razor blade together with Vicia faba L. cultivar ‘Tinova’ (2C = 26.21 pg) (Genebank Gatersleben accession number: FAB602) as an internal reference standard (Fig. 1) in a Petri dish containing 1 mL nuclei isolation buffer (Galbraith et al. 1983) supplemented with DNase-free RNase (50 µg/mL) and propidium iodide (50 µg/mL). The mixture was filtered using tubes with cell-strainer caps (35 µm; BD Biosciences, New Jersey). For species with genomes larger than 60 pg/2C, Allium cepa L. cultivar ‘Vsetatska’ (2C = 33.69 pg, Genebank Gatersleben accession number: ALL47), and with genomes smaller than 30 pg/2C, Secale cereale L. subsp. cereale (2C = 16.01 pg, Genebank Gatersleben accession number: R737) were used as internal reference standards. The relative fluorescence intensities of stained nuclei were measured using either a FACStarPLUS (BD Biosciences) flow sorter equipped with an argon ion laser INNOVA 90C (Coherent GmbH, Dieburg, Germany) or a FACSAria flow sorter (BD Biosciences). Usually, 10,000 nuclei per sample were analyzed. Accessions with the most contrasting values were measured several times. The absolute DNA amounts of samples were calculated based on the values of the G1 peak means.

Chromosome Counts—Chromosome numbers were determined from root tip meristems for 34 accessions of 14 species (Table 1). The growing root tips were pre-treated in ice-cold (0°C) tap water for 24 hr. Root tips were then fixed in acetic acid:ethanol mixture (1:3) for 24 hr. After washing in distilled water, and hydrolyzing in 0.2N HCl for 20 minutes at room temperature, roots were stained with a saturated acetocarmine solution (carmine boiled in 45% acetic acid and filtered to remove undissolved crystals) and squashed in 45% acetic acid.

Phylogenetic Analyses—Phylogenetic relationships of Melanocrommyum species were based on a study by Gurushidze et al. (2008), which used sequences of the nrDNA ITS to infer phylogenetic relationships. The analyses presented here were based on a simplified Bayesian tree (Fig. 2) in which only a single individual per species was retained, and which depicts the estimated ancestral genome sizes as well as the 2C DNA values of the measured taxa. In cases when the individuals of one species showed topological differences in the phylogenetic tree, the individuals with different topologies were kept.

C-value Statistics and Evolution of Genome Size—Nuclear DNA-values (2C) for each species were calculated as mean values from at least three accessions. The statistical analyses were done both with and without known tetraploid accessions to see if well-known genome reorganization after polyploidization (Wendel 2000; Leitch et al. 2008) would significantly change the results. We used Kolmogorov-Smirnov test for every phylogenetic clade to check if the genome sizes were normally distributed. The entire data set was normally distributed, but the data did not show normal distributions within each clade. Therefore, Kruskal-Wallis test was used to determine whether DNA-value differences among phylogenetic clades were significant, defining clades resulting from the ITS phylogenetic analysis as grouping variable. Kruskal-Wallis is a nonparametric test that checks the null hypothesis that the different samples in the comparison were drawn from the same distribution or from distributions with the same median (Hill and Lewicki 2006). Both above mentioned statistical tests were conducted with SPSS 12.0 (SPSS Inc., Chicago, Illinois). However, statistical analysis of genome size differences is problematic because conventional statistical methods consider species as independent samples, although they are evolutionarily related. A more suitable approach for analyzing data that may not be statistically independent is generalized least squares (Pagel 1997, 1999). The method is implemented in the computer program BayesContinuous (within the package Bayestraits,  www.evolution.rdg.ac.uk) and can be used to test for trait evolution, correlated trait evolution, and to perform regressions for continuously varying traits (Pagel and Meade 2007). Useful feature of BayesContinuous are the three scaling parameters kappa (κ;), delta (δ), and lambda (λ), that can be estimated for a given data set and phylogeny. These scaling parameters allow tests of the mode, tempo, and phylogenetic associations of trait evolution (Pagel 2004). The kappa (κ;) parameter differentially scales individual phylogenetic branch lengths and can be used to test for a punctuated versus gradual mode of trait evolution. If κ; = 0, trait evolution is independent of the length of the branch and is consistent with a punctuated mode of evolution, while κ; > 1 indicates that longer branches contribute more to trait evolution (gradual mode). Parameter delta (δ) scales the overall path in the tree (distance from the root to the species, as well as the shared path lengths) and can be used to test for adaptive radiation. If δ > 1, then longer path lengths have contributed disproportionately to trait evolution, which is interpreted as accelerated evolution over time (species-specific adaptation), whereas an estimate of δ < 1 is the signature of an adaptive radiation. The parameter λ is particularly important, as it tests for phylogenetic signal in the comparative data, i.e. assesses the contribution of phylogeny to the covariance among species on a given trait. A value of λ > 0 indicates strong phylogenetic signal, whereas λ = 0 suggests a star phylogeny (values of taxa are independent). The parameters κ;, δ, and λ were estimated using the ‘random-walk’ model of evolution (model A). Determining the most appropriate model as well as evaluation of the above mentioned parameters were conducted using likelihood ratio statistics (Huelsenbeck and Rannala 1997).

Ancestral genome sizes within Melanocrommyum were assessed using the generalized least squares (GLS) method implemented in the computer program Compare version 4.6b,  http://compare.bio.indiana.edu/ (Martins 2004). GLS tests the evolution of continuous characters along the phylogenetic tree from extant taxa and predicts each ancestral state as the weighted average of other taxa in the phylogeny (Martins and Hansen 1997). As both programs require a completely resolved binary tree, polytomies were resolved introducing short branches according to the authors' recommendation. The predicted genome sizes with standard errors (SE) for major lineages were plotted on the phylogenetic tree (Fig. 2).

RESULTS

Genome Size Variation in Subgenus Melanocrommyum —The 2C DNA values of all measured species of subgenus Melanocrommyum varied from 26.26 pg in diploid A. fetisowii Regel to 78.73 pg in tetraploid A. cyrilli Tenore (Table 1). Several examples of the histograms of relative fluorescence intensity are shown in Fig. 1. As chromosome numbers were not available for approximately one-third of the measured taxa, we report here only the 2C values and not the monoploid genome size (Cx), as the latter corresponds to the DNA content of the chromosome base number (Greilhuber et al. 2005).

TABLE 1.

Accessions, chromosome numbers (when available), and genome sizes of the Melanocrommyum taxa used in this study. N refers to the number of analyzed individuals from each species; SD refers to standard deviation calculated from measurements of all individuals of each species. Chromosome counts are either from the present study, or taken from the literature (Ohri et al. 1998; Friesen et al. 2006). 1Measurements were based on one or two individuals (thus SD-s are not available); 2two distinct 2C values available within one taxon (for details see text); 3chromosome counts from present study. Herbarium specimens for each taxon are deposited at GAT (see Gurushidze et al. 2008).

t01a_96.gif

Continued

t01b_96.gif

Usually, genome size measurements from different accessions of the same species were similar, as the coefficient of variation was below 2.5%. However, we found conspicuous within-species genome size variability, which could be classified as (i) relatively small (6–9% of the measured genome size) and (ii) large (two-fold) differences in DNA content. The relatively small within-species genome size variation was detected in the following species: 6% of the whole genome size in A. shelkovnikovii Grossh. and A. cardiostemon Fischer et Meyer, 7.4% in A. rosenorum R. M. Fritsch, 7.8% in A. fetisowii Regel, 8% in A. darwasicum Regel, and 8.8% in A. alexeianum Regel. To test if the observed differences indeed reflected intraspecific genome size variations rather than technical artifacts, we co-chopped the leaves of accessions 6298 (38.28 pg/2C) and 1869 (41.11 pg/2C) of A. rosenorum and measured the fluorescent intensity of the nuclei mixture of these two accessions. A similar procedure was undertaken for accessions 6553 (37.25 pg/2C) and 6134 (33.94 pg/2C) of A. darwasicum. This experiment resulted in two clearly distinguishable individual peaks (Fig. 3), indicating that the genome sizes of the two corresponding accessions of the same species are unambiguously different from each other. Roughly two-fold differences in genome size were observed in two species: 74.71 pg and 39.39 pg DNA in A. akaka, and 59.81 pg and 32.49 pg DNA in A. protensum Wendelbo. To date, A. akaka has mainly been reported to be diploid (2n = 16; Pedersen and Wendelbo 1966; Pogosian 1983; Özhatay 1986), while one study noted tetraploid and even hexaploid individuals (Vakhtina 1969). We counted chromosomes of the A. akaka accession with doubled 2C DNA content and confirmed the tetraploid status of this accession (2n = 4x = 32; Table 1). Although chromosome counts for A. protensum were available only for one accession (2n = 16, Pedersen and Wendelbo 1966), the two-fold difference in DNA amount also suggests the occurrence of both diploid and tetraploid cytotypes in this species.

In A. saposhnikovii E. Nikitina (60.15 pg) and A. sewerzowii Regel (56.45 pg) the measured genome sizes were relatively large and corresponded to a doubled amount compared to closely related taxa (A. tashkenticum Khassanov et R. M. Fritsch 0.44 pg and A. severtzovioides R. M. Fritsch 32.09 pg). However, to date, both species (A. saposhnikovii and A. sewerzowii) have been reported to be diploid (Pogosian and Seisums 1992; Fritsch and Astanova 1998).

Genome Size Variation Among Phylogenetic Clades—The phylogenetic tree based on ITS sequences of 110 species of subgenus Melanocrommyum consists of a “basal” grade (eight species) and a core clade consisting of the remaining species. In total, 12 well-supported groups (A-L) were circumscribed according to the ITS phylogeny (Fig. 2) (Gurushidze et al. 2008). The variation of genome size among the defined phylogenetic clades is represented in a box plot (Fig. 4). Table 2 shows the mean genome sizes, number of the species, maximum and minimum genome size values, and standard deviations for each phylogenetic group. The genome sizes in the “basal” grade were generally lower than in the core clade (33.65 ± 5.72 pg DNA and 42.01 ± 8.73 pg DNA, respectively). The basal-most group A possessed the lowest (29.70 ± 4.11 pg) and group C, belonging to the core clade, the highest (50.21 ± 5.71 pg) average DNA values (Fig. 4). The latter group unites morphologically similar circum-Mediterranean taxa, all belonging to section Melanocrommyum, and also contains the tetraploid species A. cyrilli.

The Kruskal-Wallis tests showed that genome size differences among the phylogenetic clades were significant. The ranking of the groups differed from one clade to another (Kruskal-Wallis χ2 = 40.19, df = 11, p < 0.001).

Genome Size Evolution and Ancestral Genome Size Estimation —The presence of a phylogenetic signal was suggested by the generalized least squares (GLS), which predicted the maximum likelihood value of λ to be significantly greater than 0. The maximum likelihood estimates of the parameters λ, κ;, δ, and respective interpretations are given in Table 3. The results of GLS (values of κ; and δ) supported a punctuated mode of genome size evolution in Melanocrommyum, and suggests that most of the genome size variation evolved later rather than earlier in the phylogeny of the subgenus (species-specific adaptations).

FIG. 1.

Histograms of flow-cytometric genome size measurements in taxa of Allium subgenus Melanocrommyum. Examples of relative fluorescence intensity of propidium-iodide stained nuclei of the species A. taeniopetalum Popov & Vved., A. robustum Kar. et Kir., A. bakhtiaricum Regel, A. pseudobodeanum R. M. Fritsch et Matin measured with the internal standard Vicia faba.

f01_96.jpg

Ancestral genome size estimation showed that genome size both increased and decreased during evolution along the lineages leading to extant species. In the “basal” grade most of the taxa show relatively low DNA contents, and GLS predicted for the root of the subgenus a 2C DNA content of 35.68 ± 6.21 pg, which is smaller than genome sizes of most contemporary taxa. However, it is obvious that both gain and loss could have occurred in closely related lineages. For example, in clades F, I, and K genome size increased in one sister lineage, while it decreased in another (Fig. 2).

FIG. 2.

Phylogenetic tree resulting from a Bayesian analysis of ITS sequences of subgenus Melanocrommyum (from Gurushidze et al. 2008). Gray double-lines indicate branches not present in the strict consensus (MP) tree. The numbers next to the species names depict 2C DNA values and standard deviations of respective species (wherever available), and the numbers along the branches show estimated ancestral genome sizes with standard errors (SE). Species for which genome sizes are not available, are depicted by #, the well-supported groups are indicated by uppercase letters on the right side (clades A - L).

f02_96.jpg

FIG. 3.

A histogram of relative fluorescent intensity of co-measured A. rosenorum individuals with deviating genome sizes. The histogram shows two clearly distinguishable individual peaks (38.28 pg and 41.11 pg), indicating real genome size differences between individuals.

f03_96.jpg

DISCUSSION

Genome Size Variation and Ploidy Level —The high genome size values (2C DNA content ranging from 26.26–78.73 pg) in Melanocrommyum are in agreement with previous reports (Ohri et al. 1998; Ohri and Pistrick 2001). The observation that plants with large genomes have short periods of growth followed by long periods of dormancy (e.g. geophytes in the orders Asparagales and Liliales; Chase et al. 2005) holds also in Melanocrommyum. Although the underlying mechanisms are unknown, the subgenus Melanocrommyum is another example of plants having an extremely short developmental period and high amounts of DNA. Genome size differences across species are more than two-fold. Although it is clear, that all known tetraploids possess doubled DNA amount compared to confirmed diploid species (Table 1), insufficient information about the ploidy of many taxa with high amounts of DNA does not allow for conclusions as to whether the genome sizes among species are always correlated to ploidy level. For example, the high 2C DNA values in A. saposhnikovii (60.15 pg) and A. sewerzowii (56.45 pg) could indicate that tetraploid cytotypes occur in these species, taking into account that the closely related species have only half of those genome sizes (30–34 pg). However, the tetraploidy of these individuals first needs to be confirmed by chromosome counts. In contrast, within A. protensum where DNA content differs by two-fold between different accessions, the occurrence of diploid and tetraploid cytotypes is highly likely, and in A. akaka the detected doubled DNA amount corresponds to a confirmed tetraploid cytotype.

FIG. 4.

Box plot displaying the distribution of 2C DNA values (Y axis) across the defined phylogenetic clades (X - axis), resulting from the ITS analyses (phylogenetic relationships among clades A - L are depicted on Fig. 2). A solid box represents the median, whiskers standard deviation, and Arabic numerals the number of the species in the respective clade.

f04_96.jpg

TABLE 2.

Genome size amount (pg) in phylogenetic clades.

t02_96.gif

Taxonomic Relevance of Genome Size Variation—Pronounced within-species differences in genome size are often reported (e.g. Price et al. 1981; Graham et al. 1994; Bennett and Leitch 1995) and correlated with differences in the habitat (e.g. Kalendar et al. 2000; Jakob et al. 2004) or in plant phenotype (Knight et al. 2005; Murray 2005; Beaulieu et al. 2007). However, as noted by Obermayer and Greilhuber (2005), intraspecific variation in genome size must be taken with utmost caution, as some earlier studies probably reported differences caused by technical artifacts. Nevertheless, intraspecific genome size variations were also found to indicate microevolutionary differentiation and could be taxonomically significant (Murray 2005; Achigan-Dako et al. 2008). This was demonstrated in Juncus biglumis L. (Schönswetter et al. 2007), where genome size differences and molecular markers revealed congruent groups.

TABLE 3.

Maximum likelihood parameter estimates for λ, δ, and κ; calculated using the generalized least squares (Pagel 1997, 1999).

t03_96.gif

Although the aim of this paper was not to analyze intraspecific genome size variation in Melanocrommyum, in some taxa within-species differences were observed. In A. rosenorum and A. darwasicum, morphologically cryptic lineages were phylogenetically clearly separated by molecular data (ITS sequences and fingerprint markers; Gurushidze et al. 2008), and they also differed in genome sizes. Two species showing intraspecific genome size variation (A. cardiostemon and A. fetisowii) are not polyphyletic according to ITS sequences. However, the genetic differentiation at the ITS locus of these species is high (0.5% and 2.9% genetic distance, respectively). Thus, genetic divergence coupled with genome size differences within A. rosenorum and A. darwasicum indicates separately evolving lineages within morphologically indistinguishable (cryptic) species, while in A. cardiostemon and A. fetisowii genome size differences could suggest the process of ongoing microevolutionary differentiation. These examples underline the potential of genome size variation for species delimitation. There are several other examples (e.g. Hordeum spontaneum K. Koch, Kalendar et al. 2000; Dactylis glomerata L., Reeves et al. 1998; Juncus biglumis L., Schönswetter et al. 2007) where intraspecific genome size variation might be indicative of incipient speciation, and thus have taxonomic relevance.

Another explanation for the congruent differences of genome size and ITS sequences within the above-mentioned species could be ancient interspecific hybridization. Thorough population genetic analyses with sampling multiple individuals and analyzing additional unlinked loci would be required to distinguish between ongoing microevolutionary diversification and ancient hybridization in these cases.

In A. shelkovnikovii (6% difference in 2C DNA content) and A. alexeianum (8.8%), ITS data provided no indication for intraspecific differentiation. Nevertheless, the range of variation of genome size between different accessions seems high compared to the “normal” range of variation (up to 3% of measured genome size for most taxa). Although it is possible that the genetic differentiation within these species is too young for the resolution of the molecular markers used (Gurushidze et al. 2008, 2010), the hypothesis that differences in genome size in A. shelkovnikovii and A. alexeianum are the result of ancient introgression, is also plausible. Although we obtained no indication for the occurrence of B chromosomes in the karyotyped species, Fritsch and Astanova (1998) reported the presence of B chromosomes in one species of Melanocrommyum, thus the hypothesis that differences in genome size in A. shelkovnikovii and A. alexeianum are due to B chromosomes cannot be ruled out.

Genome Size Changes in Relation to Phylogeny—Both analyses, the generalized least squares and Kruskal-Wallis test, support the assumption that distribution of genome size is not random but rather dependent on phylogenetic relationships. The outcomes of the statistical analyses are comparable for both data sets, including and excluding polyploids. The test for phylogenetic signal was positive using GLS, suggesting that genome size variation among species is adequately predicted by our phylogenetic tree. Based on comparison of the genome size between the “basal” grade and core clade, and estimated ancestral and contemporary genome size, the data suggest that the ancestral genome size in subgenus Melanocrommyum was relatively small. However, genome size changes do not always show a trend of increasing size, as both DNA gain and loss were detected even within closely related lineages. This is in agreement with other evolutionary reconstructions showing that both increases and decreases have occurred in different taxonomic groups (e.g. Wendel et al. 2002; Jakob et al. 2004).

According to GLS analyses, the evolution of genome size in subgenus Melanocrommyum follows the punctuated rather than gradual mode. In other taxa, the hypothesis of punctuated genome size evolution received support in Orobanche L. (Weiss-Schneeweiss et al. 2006) and Liliaceae (Leitch et al. 2007), whereas a gradual mode of evolution was inferred in Hieracium L. (Chrtek et al. 2009). Intraspecific stability of genome size (Greilhuber 1998; Gregory 2004) and reports of rapid genome reorganization after speciation (Rieseberg et al. 1995; Song et al. 1995; Özkan et al. 2001) promoted the idea that a punctuated model with most changes occurring immediately after speciation is more realistic rather than a gradual change in genome size (Albach and Greilhuber 2004; for review see Gregory 2004). However, to prove this hypothesis (punctuated mode more plausible than gradual), statistical analyses regarding the mode of the genome size evolution in a phylogenetic framework are needed in many other eukaryotic groups.

Considering the ecological uniformity of the entire subgenus, a pattern of non-adaptive evolution (δ > 1) was rather expected. However, high variability of genome size in Melanocrommyum prompts to invoke non-adaptive theories of genome size evolution. Oliver et al. (2007) have proposed a simple proportional model of genome size evolution without invoking selection pressure. The model implies that the rate of genome size evolution is simply proportional to genome size. The authors showed that the fastest rates of evolution occur in the largest genomes, and this trend is evident across 20 major eukaryote taxonomic groups. The fixation of genome size changes (particularly the ones caused by neutral or effectively neutral length mutations) might also be neutral, primarily determined by stochastic events such as drift, especially in small populations or after founder events. Unfortunately, population genetic studies in non-model groups of organisms like Allium are largely missing, which makes this challenging hypothesis currently untestable.

ACKNOWLEDGMENTS.

The authors thank Christina Koch and Ina Faustmann for expert technical assistance. We are grateful to Reinhard Fritsch and Enoch Achigan-Dako for fruitful discussions on an early version of the manuscript. We also thank Anne Bruneau and two anonymous reviewers for their comments and suggestions to improve the manuscript. The first author was partly funded by VolkswagenStiftung (Hannover, Germany) in the frame of “PharmAll” project.

LITERATURE CITED

  1. E. G. Achigan-Dako , J. Fuchs , A. Ahanchede , and F. R. Blattner . 2008. Flow cytometric analysis in Lagenaria siceraria (Cucurbitaceae) indicates correlation of genome size with usage types and growing elevation. Plant Systematics and Evolution 276: 9–19. Google Scholar

  2. D. A. Albach and J. Greilhuber . 2004. Genome size variation and evolution in Veronica. Annals of Botany 94: 897–911. Google Scholar

  3. J. M. Beaulieu , T. A. Moles , I. J. Leitch , M. D. Bennett , J. B. Dickie , and C. A. Knight . 2007. Correlated evolution of genome size and seed mass. The New Phytologist 173: 422–437. Google Scholar

  4. M. D. Bennett 1972. Nuclear DNA content and minimum generation time in herbaceous plants. Proceedings of the Royal Society of London. Series B. Biological Sciences 181: 109–135. Google Scholar

  5. M. D. Bennett 1987. Variation in genomic form in plants and its ecological implications. The New Phytologist 106: 177–200. Google Scholar

  6. M. D. Bennett and I. J. Leitch . 1995. Nuclear DNA amounts in angiosperms. Annals of Botany 76: 113–176. Google Scholar

  7. J. L. Bennetzen and E. A. Kellogg . 1997. Do plants have a one-way ticket to genomic obesity? The Plant Cell 9: 1509–1514. Google Scholar

  8. J. L. Bennetzen , J. X. Ma , and K. M. Devos . 2005. Mechanisms of recent genome size variation in flowering plants. Annals of Botany 95: 127–132. Google Scholar

  9. M. W. Chase , L. Hanson , V. A. Albert , W. M. Whitten , and N. H. Williams . 2005. Life history evolution and genome size in subtribe Oncidiinae (Orchidaceae). Annals of Botany 95: 191–199. Google Scholar

  10. J. Chrtek , J. Zahradníček , K. Krak , and U. Fehrer . 2009. Genome size in Hieracium subgenus Hieracium (Asteraceae) is strongly correlated with major phylogenetic groups. American Journal of Botany 104: 161–178. Google Scholar

  11. A. V. Cox , G. J. Abdelnour , M. D. Bennett , and I. J. Leitch . 1998. Genome size and karyotype evolution in the slipper orchids (Cypripedioideae: Orchidaceae). American Journal of Botany 85: 681–687. Google Scholar

  12. K. M. Devos , J. K. M. Brown , and J. L. Bennetzen . 2002. Genome size reduction through illegitimate recombination counteracts genome expansion in Arabidopsis. Genome Research 12: 1075–1079. Google Scholar

  13. N. Friesen , R. M. Fritsch , and F. R. Blattner . 2006. Phylogeny and new intrageneric classification of Allium L. (Alliaceae) based on nuclear ribosomal DNA ITS sequences. Aliso 22: 372–395. Google Scholar

  14. R. M. Fritsch and S. B. Astanova . 1998. Uniform karyotypes in different sections of Allium L. subgen. Melanocrommyum (Webb & Berth.) Rouy from Central Asia. Feddes Repertorium 109: 539–549. Google Scholar

  15. D. W. Galbraith , R. K. Harkins , M. J. Maddox , M. N. Ayres , P. D. Sharma , and E. Firoozabady . 1983. Rapid flow cytometry analysis of the cell cycle in intact plant tissues. Science 220: 1049–1051. Google Scholar

  16. M. J. Graham , C. D. Nickell , and A. L. Rayburn . 1994. Relationship between genome size and maturity group in soybean. Theoretical and Applied Genetics 88: 429–432. Google Scholar

  17. T. R. Gregory 2004. Macroevolution, hierarchy theory, and the C-value enigma. Paleobiology 30: 179–202. Google Scholar

  18. J. Greilhuber 1998. Intraspecific variation in genome size: a critical reassessment. Annals of Botany 82: 27–35. Google Scholar

  19. J. Greilhuber , J. Doležel , M. A. Lysák , and M. D. Bennett . 2005. The origin, evolution and proposed stabilization of the terms ‘Genome Size’ and ‘C-Value’ to describe nuclear DNA contents. Annals of Botany 95: 255–260. Google Scholar

  20. J. P. Grime 1986. Prediction of terrestrial vegetation responses to nuclear winter conditions. The International Journal of Environmental Studies 28: 11–19. Google Scholar

  21. C. E. Grover and J. F. Wendel . 2010. Recent insights into mechanisms of genome size change in plants. Le Journal de Botanique , doi:  10.1155/2010/382732Google Scholar

  22. C. E. Grover , J. S. Hawkins , and J. F. Wendel . 2007. Tobacco genomes quickly go up in smoke. The New Phytologist 175: 599–602. Google Scholar

  23. M. Gurushidze , R. M. Fritsch , and F. R. Blattner . 2008. Phylogenetic analysis of Allium subgenus Melanocrommyum infers cryptic species and demands a new sectional classification. Molecular Phylogenetics and Evolution 49: 997–1007. Google Scholar

  24. M. Gurushidze , R. M. Fritsch , and F. R. Blattner . 2010. Species level phylogeny of Allium subgenus Melanocrommyum — incomplete lineage sorting, hybridization and trnF gene duplication. Taxon 59: 829–840. Google Scholar

  25. P. Hanelt , J. Schultze-Motel , R. M. Fritsch , J. Kruse , H. I. Maaβ , H. Ohle , and K. Pistrick . 1992. Infrageneric grouping of Allium - the Gatersleben approach. Pp. 107–123 in The Genus Allium - Taxonomic problems and genetic resources. Proceedings of international symposium, 11–13 June 1991, eds. P. Hanelt , K. Hammer , and H. Knüpffer . Gatersleben: IPK. Google Scholar

  26. J. S. Hawkins , G. Hu , R. A. Rapp , J. L. Grafenberg , and J. F. Wendel . 2008. Phylogenetic determination of the pace of transposable element proliferation in plants: copia and LINElike elements in Gossypium. Genome 51: 11–18. Google Scholar

  27. T. Hill and P. Lewicki . 2006. Statistics: methods and applications. A comprehensive reference for science, industry, and data mining. Ed 1. Tulsa: StatSoft. Google Scholar

  28. J. P. Huelsenbeck and B. Rannala . 1997. Phylogenetic methods come of age: testing hypotheses in an evolutionary context. Science 276: 227–232. Google Scholar

  29. S. S. Jakob , A. Meister , and F. R. Blattner . 2004. The considerable genome size variation in Hordeum species (Poaceae) is linked to phylogeny, life form, ecology, and speciation rates. Molecular Biology and Evolution 21: 860–869. Google Scholar

  30. J. S. Johnston , A. E. Pepper , A. E. Hall , Z. J. Chen , G. Hodnett , J. Drabek , R. Lopez , and H. J. Price . 2005. Evolution of genome size in Brassicaceae. Annals of Botany 95: 229–235. Google Scholar

  31. R. Kalendar , J. Tanskanen , S. Immonen , E. Nevo , and A. H. Schulman . 2000. Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence. Proceedings of the National Academy of Sciences USA 97: 6603–6607. Google Scholar

  32. C. A. Knight and D. D. Ackerly . 2002. Variation in nuclear DNA content across environmental gradients: a quantile regression analysis. Ecology Letters 5: 66–76. Google Scholar

  33. C. A. Knight , N. A. Molinari , and D. A. Petrov . 2005. The large genome constraint hypothesis: evolution, ecology and phenotype. Annals of Botany 95: 177–190. Google Scholar

  34. I. J. Leitch and M. D. Bennett . 1997. Polyploidy in angiosperms. Trends in Plant Science 2: 470–476. Google Scholar

  35. I. J. Leitch , M. W. Chase , and M. D. Bennett . 1998. Phylogenetic analysis of DNA C-values provides evidence for a small ancestral genome size in flowering plants. Annals of Botany 82: 85–94. Google Scholar

  36. I. J. Leitch , D. E. Soltis , P. S. Soltis , and M. D. Bennett . 2005. Evolution of DNA amounts across land plants (Embryophyta). Annals of Botany 95: 207–217. Google Scholar

  37. I. J. Leitch , J. M. Beaulieu , K. Cheung , L. Hanson , M. A. Lysak , and M. F. Fay . 2007. Punctuated genome size evolution in Liliaceae. Journal of Evolutionary Biology 20: 2296–2308. Google Scholar

  38. I. J. Leitch , L. Hanson , K. Y. Lim , A. Kovarik , M. W. Chase , J. J. Clarkson , and A. R. Leitch . 2008. The ups and downs of genome size evolution in polyploid species of Nicotiana (Solanaceae). Annals of Botany 101: 805–814. Google Scholar

  39. M. Lynch and J. S. Conery . 2003. The origins of genome complexity. Science 302: 1401–1404. Google Scholar

  40. J. Ma , K. M. Devos , and J. L. Bennetzen . 2004. Analyses of LTRretrotransposon structures reveal recent and rapid genomic DNA loss in rice. Genome Research 14: 860–869. Google Scholar

  41. E. P. Martins 2004. Compare, version 4.6b. Computer program for the statistical analysis of comparative data.  http://compare.bio.indiana.edu/. Bloomington, Indiana: Department of Biology, Indiana University. Google Scholar

  42. E. P. Martins and T. F. Hansen . 1997. Phylogenies and the comparative method: a general approach to incorporating phylogenetic information into the analysis of interspecific data. American Naturalist 149: 646–667. Google Scholar

  43. M. T. Morgan 2001. Transposable element number in mixed mating populations. Genetical Research 77: 261–275. Google Scholar

  44. B. G. Murray 2005. When does intraspecific C-value variation become taxonomically significant? Annals of Botany 95: 119–125. Google Scholar

  45. R. Obermayer and J. Greilhuber . 2005. Does genome size in Dasypyrum villosum vary with fruit colour? Heredity 95: 91–95. Google Scholar

  46. D. Ohri , R. M. Fritsch , and P. Hanelt . 1998. Evolution of genome size in Allium (Alliaceae). Plant Systematics and Evolution 210: 57–86. Google Scholar

  47. D. Ohri and K. Pistrick . 2001. Phenology and genome size variation in Allium L. - a tight correlation? Plant Biology 3: 654–660. Google Scholar

  48. M. J. Oliver , D. A. Petrov , D. D. Ackerly , P. G. Falkowski , and O. M. Schofield . 2007. The mode and tempo of genome size evolution in eukaryotes. Genome Research 17: 594–601. Google Scholar

  49. N. Özhatay 1986. Kuzey Anadolu'da yeti§en Allium L. türleri ve bunların kromozom sayıları. Doğa Turk Biyoloji Dergisi 10: 452–458. Google Scholar

  50. H. Özkan , A. Levy , and M. Feldman . 2001. Allopolyploidy-induced rapid genome evolution in the wheat (Aegilops-Triticum) group. The Plant Cell 13: 1735–1747. Google Scholar

  51. M. Pagel 1997. Inferring evolutionary processes from phylogenies. Zoologica Scripta 26: 331–348. Google Scholar

  52. M. Pagel 1999. Inferring the historical patterns of biological evolution. Nature 401: 877–884. Google Scholar

  53. M. Pagel 2004. Continuous manual. Available at  http://www.evolution.rdg.ac.uk/BayesTraits.htmlGoogle Scholar

  54. M. Pagel and A. Meade . 2007. BayesTraits. Computer program and documentation available at  http://www.evolution.rdg.ac.uk/BayesTraits.htmlGoogle Scholar

  55. K. Pedersen and P. Wendelbo . 1966. Chromosome numbers of some SW Asian Allium species. Blyttia 24: 307–313. Google Scholar

  56. D. A. Petrov 2001. Evolution of genome size: new approaches to an old problem. Trends in Genetics 17: 23–28. Google Scholar

  57. D. A. Petrov 2002. DNA loss and evolution of genome size in Drosophila. Genetica 115: 81–91. Google Scholar

  58. A. I. Pogosian 1983. Chromosome numbers of some species of the genus Allium (Alliaceae) distribution in Armenia and Iran. Botanicheskii Zhurnal 68: 652–660. Google Scholar

  59. A. I. Pogosian and A. G. Seisums . 1992. Chromosome numbers of some Allium species (Alliaceae) from Afghanistan, Turkey and Central Asian states. Botanicheskii Zhurnal 77: 103–104. Google Scholar

  60. H. J. Price , K. L. Chambers , and K. Bachmann . 1981. Genome size variation in diploid Microseris bigelovii (Asteraceae). Botanical Gazette (Chicago, Ill.) 142: 156–159. Google Scholar

  61. G. Reeves , D. Francis , M. S. Davies , H. J. Rogers , and T. Hodkinson . 1998. Genome size is negatively correlated with altitude in natural populations of Dactylis glomerata. Annals of Botany 82: 99–105. Google Scholar

  62. L. H. Rieseberg , C. van Vossen , and A. M. Desrochers . 1995. Genomic reorganization accompanies hybrid speciation in wild sunflowers. Nature 375: 313–316. Google Scholar

  63. H. Schmuths , A. Meister , R. Horres , and K. Bachmann . 2004. Genome size variation among accessions of Arabidopsis thaliana. Annals of Botany 93: 317–321. Google Scholar

  64. P. Schönswetter , J. Suda , M. Popp , H. Weiss-Schneeweiss , and C. Brochmann . 2007. Circumpolar phylogeography of Juncus biglumis (Juncaceae) inferred from AFLP fingerprints, cpDNA sequences, nuclear DNA content and chromosome numbers. Molecular Phylogenetics and Evolution 42: 92–103. Google Scholar

  65. M. Shopova 1972. Citoloshki istrazhuvajna na rodot Allium od florata na Makedonijha. Godišen Zbornik Biologia 24: 83–103. Google Scholar

  66. P. Šmarda and P. Bureš . 2006. Intraspecific DNA content variability in Festuca pallens on different geographical scales and ploidy levels. Annals of Botany 98: 665–678. Google Scholar

  67. P. šmarda , P. Bureš , L. Horová , and O. Rotreklová . 2008. Intrapopulation genome size dynamics in Festuca pallens. Annals of Botany 102: 599–607. Google Scholar

  68. D. E. Soltis and P. S. Soltis . 2000. Contributions of plant molecular systematics to studies of molecular evolution. Plant Molecular Biology 42: 45–75. Google Scholar

  69. D. E. Soltis , P. S. Soltis , M. D. Bennett , and I. J. Leitch . 2003. Evolution of genome size in the angiosperms. American Journal of Botany 90: 1596–1603. Google Scholar

  70. K. Song , P. Lu , K. Tang , and T. C. Osborn . 1995. Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proceedings of the National Academy of Sciences USA 92: 7719–7723. Google Scholar

  71. L. I. Vakhtina 1969. Sravnitel'no-kariologicheskoe issledovanie nekotorykh vidov luka sekcii Molium Don. Botanicheskii Zhurnal 54: 143–153. Google Scholar

  72. H. Weiss-Schneeweiss , J. Greilhuber , and G. M. Schneeweiss . 2006. Genome size evolution in holoparasitic Orobanche (Orobanchaceae) and related genera. American Journal of Botany 93: 148–156. Google Scholar

  73. J. F. Wendel 2000. Genome evolution in polyploids. Plant Molecular Biology 42: 225–249. Google Scholar

  74. J. F. Wendel , R. C. Cronn , S. Johnston , and H. J. Price . 2002. Feast and famine in plant genomes. Genetica 115: 37–47. Google Scholar

  75. J. Wyman , S. Laliberte , and M. F. Tremblay . 1997. Nuclear DNA content variation in seeds from 22 half-sib families of jack pine (Pinus banksiana, Pinaceae). American Journal of Botany 84: 1351–1361. Google Scholar

© Copyright 2012 by the American Society of Plant Taxonomists
Maia Gurushidze, Jörg Fuchs, and Frank R. Blattner "The Evolution of Genome Size Variation in Drumstick Onions (Allium subgenus Melanocrommyum)," Systematic Botany 37(1), 96-104, (1 January 2012). https://doi.org/10.1600/036364412X616675
Published: 1 January 2012
JOURNAL ARTICLE
9 PAGES


SHARE
ARTICLE IMPACT
Back to Top