- Research article
- Open Access
Genome size expansion and the relationship between nuclear DNA content and spore size in the Asplenium monanthes fern complex (Aspleniaceae)
BMC Plant Biology volume 13, Article number: 219 (2013)
Homosporous ferns are distinctive amongst the land plant lineages for their high chromosome numbers and enigmatic genomes. Genome size measurements are an under exploited tool in homosporous ferns and show great potential to provide an overview of the mechanisms that define genome evolution in these ferns. The aim of this study is to investigate the evolution of genome size and the relationship between genome size and spore size within the apomictic Asplenium monanthes fern complex and related lineages.
Comparative analyses to test for a relationship between spore size and genome size show that they are not correlated. The data do however provide evidence for marked genome size variation between species in this group. These results indicate that Asplenium monanthes has undergone a two-fold expansion in genome size.
Our findings challenge the widely held assumption that spore size can be used to infer ploidy levels within apomictic fern complexes. We argue that the observed genome size variation is likely to have arisen via increases in both chromosome number due to polyploidy and chromosome size due to amplification of repetitive DNA (e.g. transposable elements, especially retrotransposons). However, to date the latter has not been considered to be an important process of genome evolution within homosporous ferns. We infer that genome evolution, at least in some homosporous fern lineages, is a more dynamic process than existing studies would suggest.
Homosporous ferns are renowned for their high chromosome numbers, e.g. Ophioglossum reticulatum has the highest chromosome number (2n = ca.1400) so far reported for any eukaryote . Moreover, the mean chromosome number for homosporous ferns (n = 57.05) is far higher than any other plant group (including heterosporous ferns, n = 13.6; and angiosperms, n = 15.99) [2, 3]. The reason for this disparity between land plant groups remains enigmatic and is a major focus of ongoing research in this field (see [3–5]).
The application of novel genome-wide analytical methods, including genome size analysis, is providing significant insight into the processes that shape homosporous fern genomes [3, 4, 6]. Although our knowledge of genome sizes in ferns is limited (< 1% of species have been analysed), available data suggest that patterns of genome size evolution, which include (i) polyploidisation; (ii) paleopolyploidsation; and, (iii) changes in chromosome size, are not operating uniformly across all fern lineages [6–10].
In contrast to the extreme diversity of genome sizes encountered in angiosperms, which range c. 2,400-fold , genome sizes in ferns (monilophytes) are less variable, ranging just c. 94-fold [8, 9]. In fact, if we focus within homosporous ferns, this variation in nuclear DNA content only spans c. 25-fold, from 1C = 2.95 pg in Athyrium filix-femina to 1C = 72.68 pg in Psilotum nudum var. rubra. While some of this diversity arises from polyploidy (e.g. Ophioglossum petiolatum, 2n = 32x = c. 960, and 1C = 65.55 pg, see ), genome size changes can also arise within the same ploidy level in some genera. For example, the 1.5-fold range of genome sizes encountered in Davallia have taken place at the diploid level (2n = 80) with the different genome sizes between species reflected in contrasting chromosome sizes . Such diversity may be attributed to arise through different balances between the amplification of repetitive DNA such as transposable elements (TE) (especially retrotransposons) leading to genome and chromosome size increases and DNA elimination, as frequently observed in angiosperms [13, 14].
Available cytogenetic data indicates that apart from a few examples , most homosporous ferns are characterized by possessing small and rather conserved chromosome sizes with little evidence of retrotransposon activity [3, 6, 16, 17]. These studies would indicate that genome size variation in homosporous ferns is largely driven by polyploidy. This hypothesis is supported by the relatively small variation (compared with angiosperms) reported in the monoploid genome size (1Cx-value; 1Cx = 2.95 pg - 21.02 pg) and the correlation between 2C-value (DNA content of the whole chromosome complement) and chromosome numbers . The monoploid genome size (1Cx-value) is the DNA content of one un-replicated chromosome set, whilst the holoploid genome size (1C-value) is the DNA content per haploid chromosome complement (see ).
In ferns, studies on closely related species in Dryopteris[19, 20] and Polypodium have shown that genome size can be a powerful marker for taxonomic delimitations. However, as far as we are aware, no study to date has investigated the evolution of genome size in a closely related group of fern species, or determined whether it is correlated with breeding system or any morphological traits. Studies in angiosperms have demonstrated that genome size is correlated with several ecological and morphological traits such as seed mass and stomatal density (e.g. [22–25]). Similar studies would be highly informative in ferns, as traits such as spore size and stomatal cell size are often used to infer changes in ploidy levels among closely related species [26–28]. This is based on the implicit assumption that species with higher ploidy levels (and hence larger genomes) will have larger spores, although this has never been systematically tested for genome size within a phylogenetically well defined group.
In this study we investigate the evolution of genome size and spore size within the Asplenium monanthes complex (Figure 1). The complex includes a group of closely related species whose phylogenetic relationships have been recovered, showing evidence of reticulate evolution and multiple apomictic lineages . In addition, polyploidy is known to occur based on previously reported chromosome counts in some taxa (see Methods and ).
Our work had three aims: First, we investigated the validity of using C-value measurements obtained from silica dried material. Here, we report the extent of differences in C-value estimates obtained from measurements of fresh and silica dried material (stored for 6 months and/or 2 years). Second, we tested the widely held assumption that DNA amount and spore size are correlated. For this, we conducted comparative analyses between DNA content and spore length for multiple taxa within the A.monanthes complex. Third, we examined genome size evolution and determined to what extent genome size variation reflected changes in chromosome size or ploidy levels within this group of ferns. To do this, genome size estimations were made for taxa with known ploidy levels (i.e. karyologically determined in previous studies) to infer DNA ploidies for those species without reported chromosome counts. The data were then analysed within the phylogenetic framework of Dyer et al.  to provide insights into genome size dynamics.
Shifts in genome size associated with the preservation of plant material
Since fresh material was not available for all the specimens studied, the impact of silica drying on the relative fluorescence of nuclei and hence genome size estimates was investigated. Genome size estimates were very similar between fresh leaf material and material that had been stored in silica for 6 months. This was shown for A.formosum where genome size estimates for RD28 and RD33 (2C = 13.46 pg and 13.19 pg respectively) were very similar to those obtained from the 6 month old silica dried samples of A.formosum RD28, RD33, RD157 and RD158 (2C = 13.95 pg, 13.82 pg, 13.93 pg and 13.87 pg respectively) (Table 1). We therefore inferred no meaningful shift in the genome size measurements for specimens stored in silica for six months and thus genome size estimates obtained from four 6 month old silica samples (A.formosum, RD157 and RD158; spec.nov.1, RD162; and spec.nov.2, RD163), were used in all analyses as no fresh material was available.
In contrast, samples processed using 2 year old silica material showed that 2C-values were on average, 18.6% higher (ranging from 7.83–30.09%), compared with estimates made from fresh material (Table 1). Furthermore, the quality of the measurements made on 2 year old silica samples was notably lower compared with both fresh samples and samples 6 months in age, based on the mean CV% of the peaks in the flow histograms (i.e. fresh samples = 3.37; silica samples 6 month old = 3.25; Silica samples 2 years old = 5.65). Data from 2 year old material were therefore excluded from further analysis.
Variation in genome size and spore size
Within the A.formosum and A.resiliens clades, 2C-values were highly conserved (13.19-13.93 pg and 13.70-14.50 pg respectively), as were the spore lengths (27.85-34.56 μm and 42.26-45.07 μm respectively) (Figure 1). All accessions of A.resiliens formed a well supported monophyletic lineage (clade A, Figure 1), with accessions of A.aff.heterochroum occupying a sister position (clade B, Figure 1). Of the two specimens investigated of A.aff.heterochroum, one (RD75, 2C = 14.38 pg) had a similar 2C DNA content to A.resiliens (mean 2C of four individuals = 13.99 pg), while the other was distinctly higher (RD9a, 2C = 18.67 pg). In contrast to this genome size variability, both specimens showed very similar spore lengths (37.43-38.62 μm), which were smaller compared with the mean value for the four A.resiliens specimens analysed (43.29 μm).
Within the A.monanthes clade, the smallest genome size was found in the diploid spec.nov.1, which has been shown to be sister to the rest of the clade (Figure 1). Diploid sexual, spec.nov.2 showed a distinctly larger genome size (22.20 pg). The specimens of A.hallbergii and A.aff.hallbergii, which form a sister clade to the main clades of A.monanthes (MO1 and MO2), showed very little variation in nuclear DNA contents (2C = 26.00-26.23 pg), but a relatively large variation in spore length (41.33-48.69 μm). Within the A.monanthes MO1 lineage, three sub-clades with similar, but distinct 2C value ranges could be identified: MO1 (I), 29.68-30.72 pg; MO1 (II), 27.96-28.03 pg; and MO1 (III), 26.45-27.41 pg. Spore size across these three sub-clades varied from 38.09 μm to 48.54 μm. The nuclear DNA content of individuals within the A.monanthes (MO2) clade was found to be similar (2C = 26.09-27.43 pg), with the exception of a single specimen (RD96, designated MO2 II), which had a significantly smaller genome size (2C = 18.58 pg) (Figure 1 and Table 2). Spore size varied from 38.94-48.54 μm, and the value for specimen RD96 also fell within this range (42.11 μm).
Asplenium fibrillosum, A.soleirolioides, and A.polyphyllum, which together comprise the A.castaneum clade, were observed to have the highest 2C values in the complex (37.78 pg for A.fibrillosum, 29.95 pg for A.polyphyllum, and 50.29 pg for A.soleirolioides) (Table 2 and Table 1), but since these estimates were made from material that had been stored on silica for 2 years their reliability is questioned (as noted above). Interestingly, the mean spore length of the A.fibrillosum specimen was also very large (47.90 μm), although those for A.soleirolioides and A.polyphyllum were distinctly smaller (39.37 μm and 35.42 μm respectively).
Correlation between DNA amount and spore size
Regression analysis of the raw data showed a strong positive correlation (p = 4.57×10-6, r2 = 0.542, slope = 10.986, 95% CI’s = 8.665-13.929) between DNA amount and spore length (Figure 2A). However, phylogenetically independent contrasts analysis indicated no significant relationship between contrasts (Grafen branch transformation, p = 0.101) (Figure 2B and Table 3). This indicated that divergences in DNA amount are not associated with divergences in spore length, and therefore there is no correlated evolution between DNA amount and spore length.
Genome size and ploidy level
Monoploid genome size (1Cx-values) and ploidy levels were not inferred for A.fibrillosum, A.soleirolioides and A.polyphyllum, due to the shifts observed in the 2C DNA contents of these species, as a consequence of their storage for 2 years in silica (see above). Ploidy levels and 1Cx-values were inferred for all other subclades analysed (Table 4 and Figure 3).
Within the A.resiliens clade, different ploidy levels were inferred for the two A.aff.heterochroum specimens. One specimen of A.aff.heterochroum (RD75; 2C = 14.38 pg) had a similar 2C-value to A.resiliens (mean 2C = 13.99 pg), a species reported to be triploid, suggesting that RD75 could also be a triploid. In contrast, the other A.aff.heterochroum specimen analysed (RD9a) had a significantly higher 2C-value (18.67 pg) and a 1Cx-value approximately four times that inferred for RD75, which lead to the inference that this specimen was a tetraploid.
The suggestion that most A.monanthes specimens are triploid is based on the observation that, with the exception of one report of a tetraploid count, all counts for this species have been 2n = 3x = 108. Under this assumption, the mean 1Cx-values for MO1 sub-clades I, II and III are 10.06 pg, 9.33 g, and 8.95 pg, respectively (Table 4). A triploid level is also assumed for MO2 (I), based on the 2C-values, resulting in a mean 1Cx-value of 9.02 pg. The A.monanthes MO2 (II) sub-clade, which comprised a single specimen (RD96), had a significantly lower 2C-value of 18.58 pg, and we infer that this specimen is likely to be a diploid cytotype. This would result in a 1Cx value of 9.29 pg, which is in broad agreement with the monoploid genome sizes calculated for the remaining A.monanthes sub-clades (Table 4). Nevertheless, it is noted that the mean spore size of RD96 (42.11 μm) conflicts with this finding, suggesting that this specimen is a polyploid.
In A.hallbergii, the similarity of its mean 2C-value (26.11 pg) and mean spore size (45.24 μm) to that of the triploid A.monanthes lineages (MO1 and MO2.I) suggests this species is also a triploid. As for the two new species, spec.nov.1 and spec.nov.2, they differ in 2C-values (17.67 pg and 22.20 pg respectively) but have similarly small spore sizes. They are currently inferred to be diploid although further data are needed to confirm or refute this.
Overall, the results show that monoploid genome size varies 2.4-fold (4.66-11.10 pg), with little variation in 1Cx-values within clades, but considerable variation between some of them (Table 4 and Figure 3). For example, compared with the small 1Cx-values in the A.resiliens clade (4.66-4.79 pg), the mean 1Cx-value for diploid A.formosum (6.81 pg) was markedly higher. A larger monoploid genome size was also noted in the A.monanthes clade compared with the A.resiliens clade. Indeed, the largest monoploid genome sizes were encountered within the A.monanthes clade (8.70-11.10 pg), with spec.nov.2 having the highest 1Cx-value (11.10 pg) of all, which is more than double the mean monoploid genome size of the A.resiliens clade.
Considerations for the use of silica dried material for genome size estimation
The question as to whether it is possible to use silica preserved material to measure genome size in absolute units has recently become the focus of intense debate [30–32]. In order to avoid ambiguous interpretations, some authors only advocate the use of fresh material to estimate genome size , and restrict the use of silica-dried material to estimating DNA ploidy levels (e.g. ). However, recent studies have started to explore the use of silica dried material for genome size measurements (see ). This would help overcome some of the constraints imposed by the need for fresh material, especially when the quality of measurements is not compromised (i.e. <5% of variation between fresh and silica-dried tissues).
In the current analysis some of the accessions were only available as silica dried material (either 6 months or 2 years old), so the reliability of the genome size estimates obtained from such material was investigated. We were able to compare fresh and silica-dried tissues belonging to the same accession, for three separate species (Table 1). We found a significant increase in fluorescence intensity in all silica-dried samples of 2 years in age, relative to the fresh material, which resulted in genome size estimates up to c. 30% larger compared with those estimated on fresh material (Table 1). In addition, as previously reported, %CV was also seen to increase, and in most cases these exceeded acceptable values (i.e. >5%). We therefore conclude that for Asplenium, the use of 2 year old silica dried material is unsuitable for obtaining accurate genome size estimations.
In contrast, samples of two accessions of A.formosum stored for 6 months in silica gave very similar results to those obtained from fresh material (RD28 and RD33; Table 1). In addition, there was no evidence of any reduction in the quality of the flow histogram with the CVs for both fresh and silica samples being similar. Nevertheless, while these results suggest that the short term storage of Asplenium samples in silica is suitable for genome size estimations, the extent to which these findings can be extrapolated to other plant genera, including other ferns remains unknown.
Relationship between DNA amount and spore size
Numerous studies in angiosperms have shown a positive correlation between genome size and a number of morphological traits, including seed mass, cell size and stomatal density (e.g. [23, 24, 35, 36]). Beyond the angiosperms, such studies are relatively rare (e.g. [19, 21]), and while spore size is often used as a proxy for inferring ploidy level (e.g. [26–28]), there have been no empirical studies to date that test these inferences using genome size. In this study we investigated the relationship between genome size and spore size to determine the extent to which such an assumption is valid.
We found a significant and positive correlation between nuclear DNA amount and spore length using the raw data. However, when evolutionary relationships were considered, using phylogenetically independent contrasts (PICs), no significant correlation was found. The discrepancy between analyses is interesting, and although the significance of the raw data should not be discounted it does highlight the importance of using PICs to determine the evolutionary association of statistically non-independent traits .
The insignificant relationship observed between DNA amount and spore size calls into question the use of spore size for inferences of ploidy level within homosporous ferns. Previous authors have pointed out that this inference should be restricted to very close relatives, and the distinction between diploids and their autoploid offspring [26, 27]. Our findings are consistent with these suggestions. In addition, our study sample is mainly comprised of apomictic accessions, and due to a low sample number of sexually reproducing species we are unable to test for the effect of reproductive mode on the relationship between traits.
Considerable variation was found in spore length, both within and between specimens of the apomictic lineages. In part this variation may be a consequence of the non-globulose shape of the monolete spores found in this complex, making precise measurements difficult. However, the variation may also be due to the variable nature of Döpp Manton sporogenesis in triploid apomicts of homosporous ferns [38, 39]. The variation in spore size within and between these taxa is reflected by the inability to determine ploidy level (beyond the difference between diploid and polyploids) from spore size measurements in this complex .
Overall, our findings challenge the utility of spore size for inferring ploidy level within ferns. Nevertheless, further investigations are needed to test our findings for spore size in a greater number of lineages, and also determine the effect, if any, of reproductive mode on the relationship between DNA amount and spore size.
The evolution of genome size in the Asplenium monanthescomplex
The high chromosome numbers, and conserved chromosome sizes reported for many homosporous ferns has contributed to the hypothesis that the evolution of fern genomes is less dynamic than the evolution of angiosperm genomes [3, 4, 9]. This has been suggested to be due to a higher retention rate of chromosomes and the possible suppression of transposable elements (TEs) in homosporous ferns [5, 6]. The inferred constancy of chromosome size is based on physical measurements , the low number of reported retrotransposons , and the reported correlation between chromosome number and genome size . Bainard et al.  suggest that these observations show that genome size expansion is mainly driven by polyploidy, and that stepwise increases observed in 1Cx-values across monilophytes might be indicative of paleopolyploidy.
Indeed, in the Asplenium monanthes complex studied here some of the 1C-value variation within clades clearly arises from polyploidisation, with ploidy levels ranging from 2x to 4x (Table 4), as inferred by comparing the genome size data with previously published chromosome counts (available for some of the species). This is consistent with the hypothesis that genomes size variation in homosporous ferns is driven by polyploidisation. However, contrary to this hypothesis, we also report variation in the monoploid genome size (1Cx-value), indicating chromosome size variation between species (Figure 3). Our study provides evidence that both processes of genome evolution are occurring in this complex, both are discussed below.
Polyploidy is known to occur in several species of this complex based on previously reported chromosome counts (see Methods). The genome size data obtained in this study proved useful in identifying putative cryptic species (e.g. A.aff.heterochroum, and MO2,II) and in inferring ploidy levels for taxa without reported chromosome numbers. Inferences of ploidy were carried out with caution as even closely related taxa with different ploidy levels can display similar genome sizes, potentially leading to incorrect ploidy determinations . A.formosum and the apomict A.resiliens clearly illustrated this situation as both taxa have very similar 2C-values (mean 2C = 13.61 pg and 13.99 pg respectively) (Figures 1 and 3). Without additional information it might be tempting to assume they had the same ploidy level. Yet complementary information from published chromosome counts strongly supported the presence of different ploidy levels, with A.formosum being a sexual diploid and A.resiliens a triploid apomict. Different ploidy levels were also identified among A.aff.heterochroum specimens (RD75 = 3x and RD9a = 4x), when compared to the genome size of the sister species A.resiliens.
In the A.monanthes complex moderate variation in 1C-values was noted between the accessions identified as A.monanthes (Figure 1). However, in all but one case (i.e. RD96), these specimens were inferred to be triploid. It is suggested that such variation has most likely arisen from differences in the DNA contents of the progenitor species that gave rise to the triploids, a claim supported by the different 1C-values reported for the putative diploid progenitor species, spec.nov.1 and spec.nov.2 . The exception noted in specimen RD96, which belonged to the MO2.II lineage of A.monanthes, was inferred to be a diploid apomict, due to a significantly smaller 2C-value (18.58 pg) than the mean 2C of the remaining triploid taxa (27.64 pg). This suggestion is consistent with the occurrence of diploid apomicts in other ferns, including taxa of the Dryopteris affinis complex [19, 41] and the Pteris cretica complex .
Chromosome size expansion/reduction is indicated in this complex by the observed differences in the monoploid genome size of the study species. This finding would suggest that chromosome size may not be as conserved as widely reported in homosporous ferns (e.g. ). In our study, the 1Cx-values of specimens within the A.resiliens clade (1Cx = 4.66-4.79 pg) are similar to more distantly related Asplenium species including A.trichomanes ssp. quadrivalens (1Cx = 4.53 pg; a species within A.trichomanes complex) [6, 8]. This indicates that A.formosum and species within the A.monanthes clade have undergone up to a two-fold expansion in monoploid genome size.
An expansion of the chromosome/ monoploid genome size could arise via retrotransposon-driven changes, as suggested for Equisetum and Pstilotum[9, 38, 43, 44]. Retrotransposon proliferation and elimination is linked with genomic and environmental factors such as effective population size, environmental stress, hybridisation and polyploidy [13, 45–49]. Given that the apomictic A.monanthes clade and related lineages show strong patterns of reticulate evolution, it is possible that these processes may be acting as triggers for retrotransposon activity leading to the range of 1Cx-values observed. Indeed, given the extent of hybridization and reticulate evolution reported in homosporous ferns in general , it seems likely that retrotransposon driven changes in genome size are probably more widespread across ferns but may have been largely overlooked due to the low level of sampling.
Our findings indicate that the evolution of genome size and spore size are not correlated within the A.monanthes complex. These findings challenge the utility of spore size for inferences of ploidy level within ferns. However, the prevalence of apomixis within this complex may be the cause of these findings, and the effect, if any, of reproductive mode on the relationship between DNA amount and spore size is currently unclear.
Our study also provides important insight into the dynamism of the fern genome. Previous studies have suggested that chromosome size expansion plays only a minor role in the evolution of the fern genome; with most genome size variation generated by polyploidisation. These findings were based on a broad taxonomic sample but included a very low coverage (<1%) of homosporous fern species. In our analysis of genome size in a polyploid species complex, we have found evidence to suggest that genome size variation is not explained by polyploidy alone, but also by mechanisms inducing changes in the amount of DNA per chromosome (chromosome size), without altering the number of chromosomes per genome. This finding indicates the potential for retrotransposon-driven chromosome/genome size expansion within homosporous ferns. This would have large implications for our understanding of the evolution of homosporous fern genomes in general, and it highlights the need for a substantial increase in genome size sampling, in order to determine the full extent to which these processes operate across the diversity of ferns.
The Asplenium species included in the present study are listed in Table 2. They were selected based on the phylogenetic investigation of the complex by Dyer et al. and whether the species could be successfully cultivated (see Table 2). Spec.nov1 and spec.nov.2, which belong to the A.monanthes complex, are reported by Dyer et al.  to be sexual diploids (based on spore size and nuclear DNA sequence analysis) and considered to be putative progenitor species to the two distinct apomictic lineages of A.monanthes, referred to as MO1 and MO2 (see Figure 1).
Chromosome counts and derived ploidy levels are reported in the literature for: A.formosum, 2n = 2x = 72 [1, 51–54]; A.monanthes, mainly 2n = 3x = 108, although there is a single count of 2n = 4x = 144 [38, 55–60]; A.resiliens, 2n = 3x = 108 and A.heterochroum, 2n = 4x = 144/2n = 6x = 216, 2n = 5x = 180 [52, 57, 61–67].
Fresh material was unavailable for A.fibrillosum, A.polyphyllum, A.soleirolioides, spec.nov.1, spec.nov.2, and some accessions of A.formosum (see Table 2), so silica dried samples were used instead. In total, the study sample comprised 31 fresh samples cultivated from specimen spores, and seven silica-dried samples.
Flow cytometry: genome size and DNA ploidy
The nuclear DNA content of 38 specimens was measured following the one-step procedure described by Dolezel et al.. Individual sporophytes were prepared as follows: several pinnae (after removing the rachis) were co-chopped using new razor blades, together with the appropriate calibration standard (Pisum sativum ‘Ctirad’, 2C = 9.09 pg) , in a Petri dish containing 2 mL of ‘General purpose buffer’ (GPB) (Loureiro et al., ) supplemented with of 3% PVP-40 . The suspension of nuclei was then filtered through a 30 μm nylon mesh, stained with 100 μl of propidium iodide (Sigma; 1 mg · mL-1), and treated with 34 μl of 3 mg · mL-1 ribonuclease A (RNase A; Sigma). Both propidium iodide and RNase A were adjusted to a final concentration of 50 μg/mL. Samples were kept on ice for 30 min and 5,000 particles recorded using a Partec Cyflow SL3 flow cytometer (Partec GmbH) fitted with a 100 mW green solid state laser (532 nm, Cobolt Samba). Flow histograms were analysed with the FlowMax software (v. 2.4, Partec GmbH). Three sporophytes were measured separately for each specimen, and three replicates of each were processed. Measurements obtained from fresh and silica dried materials of the same specimens were compared to determine the extent to which the preservation method influenced the relative fluorescence (i.e. nuclear DNA content) estimate.
We compared 1Cx-values within species (for species with karyologically determined ploidy levels) to check for cryptic ploidy levels. In order to infer DNA ploidy levels (i.e. ploidy levels inferred only from DNA amount cf. ) for the species without chromosome counts, we compared their holoploid genome size and mean spore length (see section below), with those of taxa within the same clades whose ploidy level had been karyologically determined. Given the comprehensive variety (quantity and geographic spread) of the previously reported karyological data, we consider our inferences of 1Cx-values for the remaining taxa (by means of the DNA content) to be legitimate.
To analyse the relationship between spore size and genome size we combined the spore size data from Dyer et al. , with additional spore measurements made for the remaining specimens sampled in this study. Spores from individual specimens were mounted onto slides using glycerine jelly. Each spore length was measured using AxioVision on a calibrated light-microscope (v4.8.2, Zeiss). An average of 25 spores were measured per specimen, and a mean spore length was calculated (Table 2). Special care was taken to identify putative abortive spores that were then excluded from the analyses. Individual box plots were compiled to show spore size variance within each specimen, as well as the interquartile range and the median.
The relationship between DNA amount and spore length
We investigated the relationship between DNA amount and spore length by comparing the results obtained from analysing the raw data with those obtained using phylogenetic independent contrasts (PIC),. A standard regression analysis assumes that the data points are independent, which is not the case in related species data. PIC’s take phylogenetic relationships into account by transforming the species data into statistically independent data. Different DNA amount values were used for specimens based on their mode of reproduction: apomictic ferns produce unreduced spores, and so 2C-values were used for analysis; sexually reproducing ferns produce reduced spores, and hence 1C-values were used for analysis. The raw data were not normally distributed; therefore, in order to linearize the data for PIC analysis , the mean measurements for DNA amount were log transformed.
Regression analysis of raw data
We used a linear model (LM) regression analysis to test for the correlation between traits in all sampled specimens. We then fitted a line of best fit using a standardised major axis (SMA) to obtain a slope estimate and r2 value (Beaulieu et al., ; Connolly et al., ). The SMA reduces the residuals in both the dependent and independent variables (rather than just the dependent variable, as in the LM model), and is therefore useful here as it is unknown which variable is which. An SMA was fitted using the (S)MATR package in R [72, 73].
Regression analysis incorporating phylogenetic data
To incorporate phylogenetic information into the regression analysis using PIC, we first reconstructed a phylogenetic tree of the sampled taxa. Bayesian inference (BY) was performed on a combined matrix of three plastid regions for the species listed in Table 2 (see , for voucher and accession numbers), using substitution models determined in jModeltest according to BIC criterion . Sequence data were incomplete for some of the taxa sampled in this study (i.e. not all three plastid regions were present). On these occasions, sequence data from very closely related taxa (based on analysis of individual plastid regions) were used as substitutes in phylogenetic reconstruction (see Table 2). Analysis was carried out in MrBayes 3.1 [75, 76], with Markov Chain Monte Carlo (MCMC) run for 5 million generations and sampled every 500 generations to approximate the posterior probabilities of trees. Two analyses were run simultaneously, and a conservative burn-in phase of 25% was implemented to disregard trees prior to convergence on the maximum likelihood. Remaining trees were then compiled to give 7,500 trees for each run, from which a 50% majority rule consensus was calculated.
Phylogenetically independent contrasts (PIC)
In order to determine whether PIC analysis was appropriate, we tested for ‘phylogenetic signal’ (= κ value), i.e., trait similarity among closely related species . We used the ‘Analysis of Traits’ model in Phylocom , using the ‘Picante package’ in R , to asses phylogenetic signal for a series of different branch length transformations. Significant phylogenetic signal was shown for all branch length transformation (Table 3), supporting phylogenetic regression by the PIC method [37, 80].
PICs were calculated using the PDAP: PDTREE module in Mesquite v.2.75 [81, 82]. This method uses branch lengths to standardise contrasts between closely related taxa and is able to deal with the soft polytomies present in our phylogeny (Garland and Díaz-Uriarte ). We had several zero-length terminal branches, so it was necessary to transform branch lengths in Mesquite, in order to generate the PICs.
To test whether branch lengths had adequately standardised the contrasts, we performed regression analysis on the PICs of both characters against their standard deviation. Logarithmic, Pagel  and Nee  transformation of branch lengths resulted in an insignificant relationship for spore length, but was significant for DNA content (see Table 3) indicating that the contrasts had not been significantly standardised. However, the branch transformation methods of Grafen , showed no significant relationship between contrasts and their standard deviations for both characters (DNA content, p = 0.289; Spore length, p = 0.090) (see Table 3), indicating the contrasts had been adequately standardised by branch lengths. The contrasts were then standardised by dividing them by their respective standard deviations. The sign of the DNA amount contrasts were made positive, and the spore length contrasts were compared in the same direction across the node . Regression analysis of standardised contrasts (forced through the origin) was performed to test for the correlated evolution of traits in R [37, 81]. The slope estimate and r2 value was obtained using SMA analysis in R (as above), and forced through the origin.
Ilia J Leitch and Harald Schneider are joint senior authorship.
Coefficient of variation
Phylogenetically independent contrasts
Standardised major axis
Standardised major axis tests and routines
Bayesian information criterion
Markov Chain Monte Carlo.
Ghatak J: Biosystematic survey of pteridophytes from shevaroy hills, South India. Nucl. 1977, 20: 105-108.
Klekowski EJ, Baker HG: Evolutionary significance of polyploidy in the Pteridophyta. Science (80- ). 1966, 153: 305-307. 10.1126/science.153.3733.305.
Nakazato T, Barker MS, Rieseberg LH, Gastony GJ: Evolution of the nuclear genome of ferns and lycophytes. Biol Evol Ferns Lycpohytes. Edited by: Ranker TA, Haufler CH. Cambridge, UK: Cambridge Unversity Press; 2008, 175-198.
Barker MS, Wolf PG: Unfurling fern biology in the genomics age. Bioscience. 2010, 60: 177-185. 10.1525/bio.2010.60.3.4.
Barker MS: Karyotype and genome evolution in pteridophytes. Plant genome Divers vol 2, Phys Struct Behav Evol plant genomes. Edited by: Leitch IJ, Greilhuber J, Doležel J, Wendel JF. 2013, Wien: Springer, 245-253.
Bainard JD, Henry TA, Bainard LD, Newmaster SG: DNA content variation in monilophytes and lycophytes: large genomes that are not endopolyploid. Chromosom Res. 2011, 19: 763-775. 10.1007/s10577-011-9228-1.
Bennett M, Leitch IJ: Nuclear DNA amounts in Pteridophytes. Ann Bot. 2001, 87: 335-345. 10.1006/anbo.2000.1339.
Pteridophyte DNA C-values database (release 4.0, Dec. 2010). http://www.kew.org/cvalues/homepage.html.
Leitch IJ, Leitch AR: Genome size diversity and evolution in land plants. Plant genome Divers vol 2, Phys Struct Behav Evol plant genomes. Edited by: Leitch IJ, Greilhuber J, Doležel J, Wendel JF. 2013, Wien: Springer, 307-322.
Obermayer R, Leitch IJ, Hanson L, Bennett MD: Nuclear DNA C-values in 30 species double the familial representation in pteridophytes. Ann Bot. 2002, 90: 209-217. 10.1093/aob/mcf167.
Pellicer J, Fay MF, Leitch IJ: The largest eukaryotic genome of them all?. Bot J Linn Soc. 2010, 164: 10-15. 10.1111/j.1095-8339.2010.01072.x.
Grime J, Hodson J, Hunt R: Comparative Plant Ecology: a Functional Approach to Common British Species. London: Unwin Hyman; 1988.
Grover CE, Wendel JF: Recent insights into mechanisms of genome size change in plants. J Bot. 2010, 1-8.
Leitch AR, Leitch IJ: Ecological and genetic factors linked to contrasting genome dynamics in seed plants. New Phytol. 2012, 194: 629-646. 10.1111/j.1469-8137.2012.04105.x.
Britton D: Chromosome studies on ferns. Am J Bot. 1953, 40: 575-583. 10.2307/2438442.
Wagner W, Wagner F: Polyploidy in pteridophytes. Polyploidy, Biol Relev Volume 13. Edited by: Lewis WH. New York: Plenum Press; 1980, 199-214.
Brandes A, Heslop-Harrison JS, Kamm A, Kubis S, Doudrick RL, Schmidt T: Comparative analysis of the chromosomal and genomic organization of Ty1-copia-like retrotransposons in pteridophytes, gymnosperms and angiosperms. Plant Mol Biol. 1997, 33: 11-21. 10.1023/A:1005797222148.
Greilhuber J, Doležel J, Lysák MA, Bennett MD: 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. 10.1093/aob/mci019.
Ekrt L, Trávníček P, Jarolímová V, Vít P, Urfus T: Genome size and morphology of the Dryopteris affinis group. Preslia. 2009, 81: 261-280.
Ekrt L, Holubová R, Trávnícek P, Suda J: Species boundaries and frequency of hybridization in the dryopteris carthusiana (dryopteridaceae) complex: a taxonomic puzzle resolved using genome size data. Am J Bot. 2010, 97: 1208-1219. 10.3732/ajb.0900206.
Bures P, Tichy L, Wang Y-F, Bartos J: Occurrence of polypodium × mantoniae and new localities for P.Interjectum in the Czech republic confirmed using flow cytometry. Preslia. 2003, 75: 293-310.
Zedek F, Smerda J, Smarda P, Bureš P: Correlated evolution of LTR retrotransposons and genome size in the genus Eleocharis. BMC Plant Biol. 2010, 10: 265-274. 10.1186/1471-2229-10-265.
Beaulieu JM, Moles AT, Leitch IJ, Bennett MD, Dickie JB, Knight CA: Correlated evolution of genome size and seed mass. New Phytol. 2007, 173: 422-437. 10.1111/j.1469-8137.2006.01919.x.
Beaulieu JM, Leitch IJ, Patel S, Pendharkar A, Knight CA: Genome size is a strong predictor of cell size and stomatal density in angiosperms. New Phytol. 2008, 179: 975-986. 10.1111/j.1469-8137.2008.02528.x.
Greilhuber J, Leitch IJ: Genome size and the phenotype. Plant genome Divers vol 2, Phys Struct Behav Evol plant genomes. Edited by: Leitch IJ, Greilhuber J, Doležel J, Wendel JF. 2013, Wien: Springer, 323-344.
Moran RC: The Asplenium trichomanes complex in the United States and adjacent Canada. Am Fern J. 1982, 72: 5-11. 10.2307/1547078.
Barrington DS, Paris CA, Ranker TA: Systematic inferences from spore and stomate size in the ferns. Am Fern J. 1986, 76: 149-159. 10.2307/1547723.
Beck JB, Windham MD, Yatskievych G, Pryer KM: A diploids-first approach to species delimitation and interpreting polyploid evolution in the fern genus Astrolepis (Pteridaceae). Syst Bot. 2010, 35: 223-234. 10.1600/036364410791638388.
Dyer RJ, Savolainen V, Schneider H: Apomixis and reticulate evolution in the Asplenium monanthes fern complex. Ann Bot. 2012, 110: 1515-1529. 10.1093/aob/mcs202.
Bainard JD, Husband BC, Baldwin SJ, Fazekas AJ, Gregory TR, Newmaster SG, Kron P: The effects of rapid desiccation on estimates of plant genome size. Chromosom Res. 2011, 19: 825-842. 10.1007/s10577-011-9232-5.
Cires E, Cuesta C, Peredo EL, Revilla MÁ, Prieto JAF: Genome size variation and morphological differentiation within Ranunculus parnassifolius group (Ranunculaceae) from calcareous screes in the Northwest of Spain. Plant Syst Evol. 2009, 281: 193-208. 10.1007/s00606-009-0201-9.
Sánchez-Jiménez I, Hidalgo O, Canela MÁ, Siljak-Yakovlev S, Šolić ME, Vallès J, Garnatje T: Genome size and chromosome number in Echinops (Asteraceae, Cardueae) in the Aegean and Balkan regions: technical aspects of nuclear DNA amount assessment and genome evolution in a phylogenetic frame. Plant Syst Evol. 2012, 298: 1085-1099. 10.1007/s00606-012-0618-4.
Doležel J, Greilhuber J, Suda J: Estimation of nuclear DNA content in plants using flow cytometry. Nat Protoc. 2007, 2: 2233-2244. 10.1038/nprot.2007.310.
Suda J, Trávnícek P: Estimation of relative nuclear DNA content in dehydrated plant tissues by flow cytometry. Curr Protoc Cytom. 2006, 38: 7.30.1-7.30.14.
Hodgson JG, Sharafi M, Jalili A, Díaz S, Montserrat-Martí G, Palmer C, Cerabolini B, Pierce S, Hamzehee B, Asri Y, Jamzad Z, Wilson P, Raven JA, Band SR, Basconcelo S, Bogard A, Carter G, Charles M, Castro-Díez P, Cornelissen JHC, Funes G, Jones G, Khoshnevis M, Pérez-Harguindeguy N, Pérez-Rontomé MC, Shirvany FA, Vendramini F, Yazdani S, Abbas-Azimi R, Boustani S, et al: Stomatal vs. genome size in angiosperms: the somatic tail wagging the genomic dog?. Ann Bot. 2010, 105: 573-584. 10.1093/aob/mcq011.
Knight CA, Beaulieu JM: Genome size scaling through phenotype space. Ann Bot. 2008, 101: 759-766. 10.1093/aob/mcm321.
Garland TJ, Harvey PH, Ives AR: Procedures for the analysis of comparative data using phylogenetically independent contrasts. Syst Biol. 1992, 41: 18-32.
Manton I: Problems of Cytology and Evolution in the Pteridophyta. Cambridge: The syndics of the Cambridge University Press; 1950.
Döpp W: Die Apogamie bei Aspidium remotum Al.Br. Planta. 1932, 17: 86-152. 10.1007/BF01909625.
Suda J, Krahulcová A, Trávnícek P, Krahulec F: Ploidy level versus DNA ploidy level: an appeal for consistent terminology. Taxon. 2006, 55: 447-450. 10.2307/25065591.
Schneller J, Krattinger K: Genetic composition of Swiss and Austrian members of the apogamous Dryopteris affinis complex (Dryopteridaceae, Polypodiopsida) based on ISSR markers. Plant Syst Evol. 2010, 286: 1-6. 10.1007/s00606-009-0264-7.
Huang Y, Hsu SM, Hsieh TS, Chou HU, Chiou WL: Three Pteris species (Pteridaceae: Pteridophyta) reproduce by apogamy. Bot Stud. 2011, 52: 79-87.
Brownsey PJ, Lovis J: Chromosome numbers for the New Zealand species of Psilotum and Tmesipteris, and the phylogenetic relationships of the Psilotales. New Zeal J Bot. 1987, 25: 439-454. 10.1080/0028825X.1987.10413360.
Guillon J-M: Molecular phylogeny of horsetails (Equisetum) including chloroplast atpB sequences. J Plant Res. 2007, 120: 569-574. 10.1007/s10265-007-0088-x.
Grandbastien M-A, Audeon C, Bonnivard E, Casacuberta J, Chalhoub B, Costa A-P, Le Q, Melayah D, Petit M, Poncet C, Tam S, Van Sluys M, Mhiri C: Stress activation and genomic impact of Tnt1 retrotransposons in Solanaceae. Cytogenet Genome Res. 2005, 110: 229-241. 10.1159/000084957.
Petit M, Guidat C, Daniel J, Denis E, Montoriol E, Bui Q, Lim K, Kovarik A, Leitch AR, Grandbastien M-A, Mhiri C: Mobilization of retrotransposons in synthetic allotetraploid tobacco. New Phytol. 2010, 186: 135-147. 10.1111/j.1469-8137.2009.03140.x.
Kalendar R, Tanskanen J, Immonen S, Nevo E, Schulman AH: Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence. Proc Natl Acad Sci U S A. 2000, 97: 6603-6607. 10.1073/pnas.110587497.
Leitch IJ, Bennett M: Genome downsizing in polyploid plants. Biol J Linn Soc. 2004, 82: 651-663. 10.1111/j.1095-8312.2004.00349.x.
Lockton S, Ross-Ibarra J, Gaut BS: Demography and weak selection drive patterns of transposable element diversity in natural populations of Arabidopsis lyrata. Proc Natl Acad Sci U S A. 2008, 105: 13965-13970. 10.1073/pnas.0804671105.
Lovis J: Evolutionary Patterns and Processes in Ferns. Advances in Botanical Research. Volume 4. Edited by: Preston RD, Woolhouse HW. 1977, 402.
Manton I: Cytological information of the ferns of West Tropical Africa. Ferns Fern Allies West Trop Africa. Edited by: Alston AHG. 1959, Cambridge, 75-81.
Walker TG: A cytotaxonomic survey of the pteridophytes of Jamaica. Transacions R Soc Edinburgh. 1966, 66: 169-237. 10.1017/S0080456800023516.
Ammal LS, Bahavanandan KV: Cytological studies on the genus Asplenium Linn. Indian Fern J. 1991, 8: 69-73.
Guillén RH, Daviña JR: Estudios cromosómicos en especies de Asplenium (Aspleniaceae) de la Argentina. Darwiniana. 2005, 43: 44-51.
Manton I, Vida G: Cytology of the fern flora of Tristan da Cunha. Proc R Soc London, Ser B. 1968, 170: 361-371. 10.1098/rspb.1968.0045.
Manton I, Vida G, Gibby M: Cytology of the fern flora of Madeira. Bull Br Museum (Natural Hist Bot. 1986, 15: 123-161.
Wagner W, Farrar D, McAlpin B: Pteridology of the Highlands Biological Station Area, Southern Appalachians. J Elisha Mitchell Sci Soc. 1970, 86: 1-24.
Tryon R, Voeller B, Tryon A, Riba R: Fern biology in Mexico (A class field program). Bioscience. 1973, 23: 28-33. 10.2307/1296364.
Lovis J, Rasbach H, Rasbach K, Reichstein T: Asplenium azoricum and other ferns of the A. trichomanes group from the Azores. Am Fern J. 1977, 67: 81-93. 10.2307/1546770.
Smith AR, Mickel JT: Chromosome counts for Mexican ferns. Brittonia. 1977, 29: 391-398. 10.2307/2806481.
Morzenti V, Wagner W: Southeastern American blackstem spleenworts of the Asplenium heterochroum-resiliens complex. Assoc Southeast Biol Bull. 1962, 9: 40-41.
Morzenti V: Morphological and cytological data on southeastern United States species of the Asplenium heterochroum-resiliens complex. Am Fern J. 1966, 56: 167-177. 10.2307/1545935.
Wagner W: A biosystematic study of United States ferns - preliminary abstract. Am Fern J. 1963, 53: 1-16. 10.2307/1546805.
Wagner W: Two new species of ferns from the United states. Am Fern J. 1966, 56: 3-17. 10.2307/1546108.
Wagner W, Wagner F: Pteridophytes of the mountain lake aras, Gilles Co., Virginia: Biosystematic Studies 1964–65. Castanea. 1966, 31: 121-140.
Windham MD: The ferns of Elden Mountain, Arizona. Am Fern J. 1983, 73: 85-93. 10.2307/1546855.
Haufler CH, Soltis D: Genetic evidence suggests that homosporous ferns with high chromosome numbers are diploid. Proc Natl Acad Sci U S A. 1986, 83: 4389-4393. 10.1073/pnas.83.12.4389.
Doležel J, Greilhuber J, Lucretti S: Plant genome size estimation by flow cytometry: Inter-laboratory comparison. Ann Bot. 1998, 82: 17-26.
Loureiro J, Rodriguez E, Doležel J, Santos C, Article T: Two new nuclear isolation buffers for plant DNA flow cytometry: A test with 37 species. Ann Bot. 2007, 100: 875-888. 10.1093/aob/mcm152.
Quader S, Isvaran K, Hale RE, Miner BG, Seavy NE: Nonlinear relationships and phylogenetically independent contrasts. J Evol Biol. 2004, 17: 709-715. 10.1111/j.1420-9101.2004.00697.x.
Connolly JA, Oliver MJ, Beaulieu JM, Knight CA, Tomanek L, Moline MA: Correlated evolution of genome size and cell volume in diatoms (Bacillariophyceae). J Phycol. 2008, 44: 124-131. 10.1111/j.1529-8817.2007.00452.x.
SMATR: Standardised major axis tests and Routines, ver 2.0. http://www.bio.mq.edu.au/ecology/SMATR/.
Warton DI, Duursma RA, Falster DS, Taskinen S: SMATR 3- an R package for estimation and inference about allometric lines. Methods Ecol Evol. 2011, 3: 257-259.
Posada D: jModelTest: Phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Blomberg SP, Garland TJ, Ives AR: Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution (N Y). 2003, 57: 717-745.
Webb CO, Ackerly DD, Kembel SW: Phylocom: software for the analysis of phylogenetic community structure and trait evolution. Bioinformatics. 2008, 24: 2098-2100. 10.1093/bioinformatics/btn358.
Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, Blomberg SP, Webb CO: Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010, 26: 1463-1464. 10.1093/bioinformatics/btq166.
Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15. 10.1086/284325.
PDAP Package of Mesquite. http://mesquiteproject.org/pdap_mesquite.
Mesquite: A modular system for evolutionary analysis. Version 2.5. http://mesquiteproject.org.
Garland TJ, Díaz-Uriarte R: Polytomies and phylogenetically independent contrasts: examination of the bounded degrees of freedom approach. Syst Biol. 1999, 48: 547-558. 10.1080/106351599260139.
Pagel MD: A method for the analysis of comparative data. J Theor Biol. 1992, 156: 431-442. 10.1016/S0022-5193(05)80637-X.
Read A, Nee S: Inference from binary comparative data. J Theor Biol. 1995, 173: 99-108. 10.1006/jtbi.1995.0047.
Grafen A: The phylogenetic regression. Philos Trans R Soc Lond B Biol Sci. 1989, 326: 119-157. 10.1098/rstb.1989.0106.
We would like to thank Charlie Marsh for his assistance with statistical analysis. The work was supported by the Natural History Museum in London, the Royal Botanic Gardens at Kew, the Chinese Academy of Science, the Royal Society, and the European Research Council.
The authors declare that there are no competing interests concerning the content and publication of this manuscript.
RJD designed the study, carried out data collection and statistical analyses, and drafted the manuscript. JP participated in the design of the study, and helped with the acquisition of genome size data and the drafting of the manuscript. IJL participated in the design of the study, and helped with statistical analyses and the drafting of the manuscript. VS and HS participated in the conception of the study and drafting of the manuscript. All authors have read and approved the final manuscript.
About this article
Cite this article
Dyer, R.J., Pellicer, J., Savolainen, V. et al. Genome size expansion and the relationship between nuclear DNA content and spore size in the Asplenium monanthes fern complex (Aspleniaceae). BMC Plant Biol 13, 219 (2013) doi:10.1186/1471-2229-13-219
- Genome size
- Chromosome size
- DNA content
- Asplenium monanthes
- Spore size