Abstract

Two birch species coexist in Iceland, dwarf birch Betula nana and tree birch B. pubescens. Both species are variable morphologically, which has been shown to be due to introgressive hybridization via interspecific hybrids. The aim of this study was to examine if the introgression could be related to genome size. We characterized 42 plants from Bifröst woodland morphologically and cytogenetically. The population consisted of diploid B. nana (38%), tetraploid B. pubescens (55%), and triploid hybrids (7%). Genome size was measured from 12 plants, using Feulgen DNA image densitometry (FDM) on spring leaf buds and flow cytometry (FCM) with dormant winter twigs. The use of winter twigs for FCM is novel. The average 1C-values for diploid, triploid, and tetraploid plants were 448, 666, and 882 Mbp, respectively. Monoploid genome sizes were found to be statistically constant among ploidy levels. This stability is in contrast to the different taxonomic positions of the di- and tetraploids and also contrasts with the frequent occurrence of genome downsizing in polyploids.

1. Introduction

Two species of Betula L. coexist in Iceland [1]: the diploid dwarf birch B. nana L. and the tetraploid tree birch B. pubescens Ehrh. Betula nana, a circumpolar species [2], grows up to one metre in height with procumbent, wide-spreading stems and a strong branching habit. Betula pubescens is a European species [3], represented by the subspecies pubescens Ehrh., which may grow up to 25 m tall with single or many stems (monocormic or polycormic type), and subspecies tortuosa (Ledeb.) Nyman is a shrub or low tree found in the mountain regions of northern Europe and is therefore called “mountain birch”. Betula pubescens in Iceland has sometimes also been referred to as mountain birch due to its low stature and scrub-like growth form [4]. Elkington [5] compared birch from Iceland and Scotland and suggested that the morphological variability and the shrub-like appearance of Icelandic birch were due to genetic introgression from B. nana into B. pubescens via triploid interspecific hybrids. This was confirmed by both crossing experiments and population-based studies [1, 6]. Furthermore, these studies revealed that the introgression occurred in both directions, that is, not only from the diploid to the tetraploid species but also from the tetraploid to the diploid birch species, resulting in a high degree of variation in both B. nana and B. pubescens.

Due to this bidirectional introgression between the two Betula species, taxonomic delineation in Icelandic birch has been difficult, especially with the tree birch species B. pubescens. Morphological variation in birch plants from natural woodlands throughout Iceland was examined using both qualitative species-specific botanical characters and quantitative morphological variables [7]. The study revealed that Icelandic birch populations consisted mostly of introgressant types. Whereas about 11% of all plants examined were “morphologically pure” B. nana, not a single plant had all of the species-specific characters of B. pubescens. Indeed only 1% of the plants examined had more than 80% of the B. pubescens characters. But when the plants were sorted according to chromosome number, which was determined by direct chromosome counting of mitotic metaphases, three groups of birch emerged: diploid (2 = 28), triploid (2 = 42), and tetraploid (2 = 56). The diploid (2 = 28) group consisted mainly of B. nana and B. nana-like plants, although some diploid plants showed intermediate “hybrid” morphology. The tetraploid group consisted of B. pubescens-like plants and a significant proportion of plants that had intermediate morphology. The triploid group (about 10% of all plants examined) showed no “typical” morphology, although there was a trend towards greater similarity to B. nana. The multivariate data analysis of morphological variables was able to separate diploid B. nana from tetraploid B. pubescens most of the time, but triploid hybrids were difficult to predict;-only half of them were assigned correctly.

The three genome groups in Icelandic birch were distinguishable statistically but their morphological variation overlapped in most of the distribution, implying bidirectional introgression [7, 8]. The question is whether this introgressive hybridization has affected genome size and if so, to what extent. There might be large-scale genomic recombination and shuffling following the hybridization and back-crossing, leading to genome size changes and significant variation within and between genome groups. Genome downsizing has been shown to be common in angiosperms [9], whereby genome size does not increase in direct proportion with ploidy. Genome size expansion due to transposition has also been documented, but there seems to be a large variation in the nature and extent of opposing forces to genome obesity [10, 11], making the direction of genome size change in a given species difficult to predict. On the other hand, there may be no deviation at all in genome size proportional to ploidy if there is sufficient (strong) selection pressure against genome-wide (neutral) introgression, especially in hybrid zones. For the hybrid zones to be stable, they require some form of divergent selection (intrinsic and/or extrinsic) to maintain differentiation (species integrity) despite extensive gene flow. Divergent selection has been suggested to contribute to heterogeneous genomic divergence and ecological speciation [12] and to play an important role in hybrid zone maintenance [13]. Our morphological [7, 8] and molecular [14] studies have supported the suggestion that Iceland is a birch hybrid zone, harbouring genetic variation which may be advantageous in subarctic regions.

The aim of the present study was therefore to get the first set of genome size data on Icelandic birch, including samples from all three ploidy groups (diploid, triploid, and tetraploid) from one of the most extensive birch woodlands in Iceland. Two approaches were used to measure genome size in these samples: (A) the Feulgen DNA image densitometry (FDM) method for squashed preparations of meristematic cells in newly bursting shoot buds and (B) the flow cytometry (FCM) method for cells from winter-dormant twigs. The genome size data generated in this study will undoubtedly contribute to the understanding of the process and mechanism of introgressive hybridization in Icelandic birch.

2. Materials and Methods

2.1. Plant Materials

Plant materials in this study came from the Bifröst woodland (64.76°N/21.59°W, 62 m.a.s.l., 11°C mean July temperature), by Lake Hredavatn in Borgarnes district, West Iceland. This is part of the largest continuous woodland area still in existence today. In this woodland, like most other woodlands in Iceland, both B. nana and B. pubescens grow most often side by side, although B. nana tends to be more dominant on wet ground, whereas B. pubescens occupies drier land. A total of 42 plants were selected within an area of about one km2 and with a distance of at least 20 m between plants. The plants were marked in the field to enable a long-term study, whereby samples could be collected from the same plants in different visits over a number of years.

2.2. Morphological and Cytogenetic Characterization

For the morphological analysis, 30 leaves were collected from each plant. These were late leaves from long shoots, normally from the third to fifth positions from the shoot tip. Discrete, species-specific morphological characters were used to determine the morphological variation of B. nana and B. pubescens. Eight characters including growth form, growth habit, and characteristics of leaf shape were examined qualitatively and scored as described previously [8], whereby each character was given two or three possible scores from zero (B. nana) to one or two (B. pubescens). Scores for leaf characters of each plant were assigned using average values from 30 leaves that had been pressed, sealed in plastic, and scanned. The scores of all characters were then combined for each plant into a single value, called a morphology index. This was assigned to place B. nana at the lowest ranks and B. pubescens at the highest.

The same sets of leaves were analysed quantitatively. The quantitative analysis of leaf characters was performed for nine variables as in [7], using the leaf morphology analysis program WinFolia (Regent Instruments, Quebec, Canada). The measured characters were analysed using the multivariate analysis of variance (MANOVA) and the linear discriminant analysis (LDA) described in [15].

Samples of bursting leaf buds in early spring (May-June) were collected in the field for chromosome analysis. The chromosome number of each plant was determined by direct chromosome counting, from metaphases prepared from shoot-tip meristems using the protoplast dropping method of Anamthawat-Jónsson [16]. The chromosome preparations were stained with the fluorescent dye DAPI (4, 6-diaminophenylindole) and visualized under a Nikon epifluorescence microscope (Eclipse 800) using 1000 magnification. The images were captured with a Nikon digital camera (DXM 1200F) using 12.5 million pixels. Chromosomes were counted from 10–20 metaphases from each plant.

2.3. Genome Size Determination Using Feulgen DNA Image Densitometry

Seven samples (individuals) of birch from the Bifröst woodland (Table 1) were collected in the field in May 2007 and were airmailed as twigs with early sprouting buds under a moist condition from Reykjavik to Vienna. Buds were detached and flooded with 4% phosphate-buffered formaldehyde (pH 7) in embryo cups and cut open. Bud meristems and leaf primordia were dissected out and fixed for 1.5 h at 20°C synchronously with embryonic roots and shoots from dry seeds of Pisum sativum “Kleine Rheinländerin” for internal standardization of the Feulgen reaction. Extensive washes with acetic-methanol (1 : 3) were performed, to remove formaldehyde from the tissue. The material was then rinsed extensively and finally stored in 96% ethanol at 20°C until use.

The Feulgen reaction was done as described in [17], that is, hydrolysis in 5N HCl for 90 min. at 20°C, then Schiff reagent for 90 min., with appropriate washes between all steps. Squashes were made in 45% acetic acid with birch and pea tissue particles placed side by side under separate cover slips on one slide. The slides were frozen over a cold plate and the cover slips were removed. The slides were air-dried and were kept dust- and light-protected for a maximum of a few days.

Measurement of DNA content was performed with video-based DNA image densitometry using the CIRES, version 3.1 (Kontron, Munich), using a 63 oil objective (Zeiss), green and neutral filters, the green channel of the video camera, and “local” background determination (around each nucleus separately). Three slides per plant were measured except on one occasion where two slides were examined, and usually 100 unreplicated nuclei per slide of birch and pea each were measured for DNA content determination. For each slide the 1C DNA content of birch was calculated using the birch/pea mean 2C DNA content ratio and 4.42 pg (1C) for pea [18]. The pg/Mbp conversion factor used is 1 pg = 978 Mbp [19]. Genome sizes are presented as 1C-values (holoploid genome size) and 1Cx-values (monoploid genome size), as proposed in [20, 21]. We use the term “DNA ploidy” [22] to indicate that ploidy level was inferred from DNA amounts instead of chromosome counts.

2.4. Genome Size Determination Using Flow Cytometry

In contrast to the Feulgen method using formaldehyde fixation, the FCM approach was shown to be unsuitable for sprouting leaf buds of birch because of an abundant amount of secondary metabolites. Icelandic birch, like the mountain birch of Scandinavia, is known to produce excessive amounts of primary and secondary metabolites in its leaves, especially the herbivory-related phenolic compounds [23]. We therefore looked for different materials for the FCM in this project.

Samples of young, winter-dormant twigs were collected from nine plants in the Bifröst woodland (Table 1) in November 2009. These plants were analysed morphologically and cytogenetically, except the two samples collected at this site from plants having the growth habit of B. nana (Bn-1 and Bn-2). Four of the karyotyped plants were analysed by FDM in 2007. The twig samples were sent by airmail to Vienna for the FCM study. The sample preparation for FCM began six days after the samples were collected in the field in Iceland.

About 1-2 cm long pieces of the youngest twig internodes were decorticated (Figure 1) and the cortex discarded, while the remaining twig was kept moist. Using a sharp razor blade, the decorticated twig segments were cochopped [24] with Solanum pseudocapsicum leaves (1.29 pg DNA/1C [25] for internal standardization) in Otto’s buffer I [26, 27]. The resulting nuclei isolate was filtered through a 30  m mesh (Saatilene) and treated for 30 min. with RNase-A (0.15 mg/mL) at 37°C in a water bath. Thereafter, propidium iodide (PI, final concentration 50  g/mL) dissolved in Otto’s buffer II [26] was added to stain the nuclei at neutral pH for at least 1 h in a refrigerator. Decorticated winter twigs yielded enough nuclei for suitable histograms (Figure 2). Good separation of the 2C peaks of Betula and Solanum from debris is shown by the SSC versus fluorescence analysis in Figure 2 (insert). Tests for the interference of secondary compounds with nuclei fluorescence were made with chopped birch twigs and standard leaves alone and in combination. The measurements were made at constant instrument settings [27].

Genome size measurements were done using FloMax software on a CyFlow ML flow cytometer (Partec, Münster, Germany) equipped with a green laser (100 mW, 532 nm, Cobolt Samba, Cobolt AB, Stockholm, Sweden). The C-values were calculated for every run from the ratio of the mean fluorescence intensities of the standard and the birch G1 peaks. Usually the CV% of single peaks was lower than 3% and three runs (5000 particles/run) per preparation were sufficient. If the CV% was higher (up to maximal 5%), an additional one or two runs were done. Statistical tests (ANOVA and Wilcoxon pair test) were made using SPSS 10.

3. Results

3.1. Ploidy and Morphological Characterization

The plants from Bifröst were characterized morphologically and cytogenetically (Figure 3) and were localized onto the quantitative introgression map of this woodland (Figure 4). Of the 42 plants examined, 16 plants (38%) were confirmed with accurate chromosome counts as being diploid (2 = 2 = 28), 23 plants (55%) as tetraploid (2 = 4 = 56), and 3 plants (7%) as triploid with 42 chromosomes (see examples in Figures 3(d)3(f)). The percentages of different ploidy groups within the Bifröst population are similar to those found in all woodlands combined [7]. No aneuploid cells or aneuploid plants were detected. Ten of these plants were used in the genome size study (see Table 1): three diploid Betula nana, two triploid and five tetraploid B. pubescens plants.

The plants from Bifröst woodland represent the overall features of Icelandic birch in general. The plants are mostly morphologically introgressed birch although they are remarkably stable at the ploidy level, that is, the population consists only of diploid, triploid, and tetraploid plants, no aneuploids. But to what extent has Icelandic birch been modified, via introgressive hybridization? Botanically [28] the two Betula species differ mainly in that the dwarf birch species B. nana has orbicular or obovate-orbicular sessile leaves which are regularly crenate and are rounded at both leaf tip and base (similar to the leaves in Figure 3(a)), whereas the leaf shape of B. pubescens is usually cordate with dentate margins and a distinct petiole (similar to that of the tetraploid birch in Figure 3(c)). In the morphology index scale from zero for B. nana to 13 for B. pubescens (see Materials & Methods), given null introgression, the Icelandic diploid, triploid, and tetraploid groups from Bifröst were found to have average scores of 1.4 (higher than 0), 4.0 (relative intermediate), and 8.4 (much lower than 13), respectively. These are about the same average scores as those obtained from all major birch woodlands in Iceland combined [7]. The bidirectional introgression has clearly brought the two species closer to each other. Apart from these mean values, however, there is a considerable variation within ploidy groups. The diploid B. nana plants from Bifröst, for example, have a morphology index of 2 (Figure 3(a)), but 1 and 3 in the plants used in our genome size study (Table 1). The triploid plant in Figure 3(b) has a morphology index of 4, like the other two triploids examined in this study (Table 1). The tetraploid B. pubescens plant in Figure 3(c) has the typical index of 9 for this species in Icelandic woodlands, whereas those used in the genome size study show 6–10 in their morphology scale (Table 1).

The majority of the plants that were recorded as B. nana and B. pubescens in the field were found to be diploid and tetraploid, respectively. Triploid individuals, on the other hand, resembled either species as frequently as they appeared as morphological intermediates. But when quantitative species-specific characteristics of leaf shape were analysed, the plants from Bifröst woodland were separated into two groups only, that is, diploid and tetraploid groups (Figure 4). The triploids were not statistically distinguishable from B. pubescens. The same analyses applied to triploid plants from woodlands throughout Iceland: the triploids did not form an intermediate group [7].

3.2. Genome Size in Icelandic Birch
3.2.1. Feulgen DNA Image Densitometry (FDM)

Formaldehyde fixation and the Feulgen procedure resulted in delicately magenta-stained and often isolated birch nuclei, and a colourless transparent cytoplasm, in which solid bodies of polymerized condensed tannins were abundantly present in many but not all cells (compare with Greilhuber [29]). These tannin bodies were localized within the vacuoles. None or very little mitotic activity was observed in the slides. Mitotic nuclei, if found, conformed to the corresponding interphase 2C and 4C distributions. Somewhat higher CVs within some slides could be attributed to some background granulation. In the cells of the resting pea embryos used as standard, no tannins were found. Their nuclei were nonreplicating and conformed well to ploidy classes. Only 2C nuclei were used for calculation. Within-slide CVs in birch were mostly below 6%, with one exception of 7%, while in pea they were 1.5% to 3%.

Of seven birch samples (individual plants), three were “DNA diploid” with 441.4 Mbp on average, one was “DNA triploid” with 670.8 Mbp, and three were “DNA tetraploid” with 879.2 Mbp (1C) on average (Tables 1 and 2). The 1Cx-value was fairly constant between individuals at 441.4 Mbp (CV = 1.6%, = 7, nonsignificant difference between individuals upon two-sided monofactorial ANOVA, = 1.09, = 0.418), and the variances were homogeneous (Levene test, = 0.177). Likewise, different ploidy levels did not differ significantly in their Cx-values ( = 0.340, = 0.730).

3.2.2. Flow Cytometry Using Decorticated Twigs

As far as we know, the use of woody plant twigs in FCM for otherwise recalcitrant material is new. The test of interference of phenolics showed unexpectedly somewhat increased fluorescence and improved peak CVs in both birch and standard, if chopped in combination, compared with lone isolates of both. In five tests, the increase in nuclei fluorescence was 1.043-fold in the standard and 1.038-fold in Betula pubescens. This is at present difficult to explain and needs further investigation, but the tests seem to indicate a compensating effect of the joint preparation which justifies the internal standardization.

The results from nine accessions examined clearly disclosed the respective ploidy levels: two were “DNA diploid” with 457.6 Mbp on average, two were “DNA triploid” with 663.8 Mbp, and five were “DNA tetraploid” with 883.9 Mbp (1C) on average (Tables 1 and 2). The ploidy status derived from FCM again corresponded unambiguously with the chromosome number obtained cytogenetically. The two samples that had the growth habit of B. nana in the field (no leaves in mid-winter for botanical identification) were correctly identified as belonging to the diploid species (Table 1). The Cx-values of these three ploidy levels were also subjected to ANOVA, but the differences were found to be nonsignificant ( = 3.048, = 0.122). However, a weak but real downsizing cannot be rejected and this might be detected in the future using a larger sample size. The mean Cx-value was 0.966-fold in tetraploids compared with diploids, and the intermediate Cx-value in triploids was very close to the theoretical expectation.

3.2.3. Flow Cytometry Compared with Feulgen Densitometry

From the four individuals (one triploid and three tetraploid) measured using both techniques (FDM and FCM), no technique-dependent trend in C-values can be deduced (Table 3). The nonparametric Wilcoxon test for paired samples revealed a two-sided error significance of 0.715, and so both methods yielded equivalent results. All FCM-derived Cx-values compared with FDM-derived ones showed that the individuals in both groups were not statistically different (df1,14   = 0.894, = 0.360). The Cx-values of the three ploidy levels obtained with either method (i.e., all 2 , 3 , and 4 data compared) were also not different (df2,13   = 0.966, = 0.406). The overall mean 1Cx-value from all three ploidy groups is 444 Mbp (Figure 5).

In summary, based on the average 1C-values obtained, the 1C genomes of diploid Betula nana, triploid interspecific hybrid and tetraploid B. pubescens were found to be 448, 666 and 882 Mbp DNA, respectively (Table 2). The mean C-value of the tetraploid species B. pubescens was therefore very close to being twice (1.97-fold) that of the diploid species B. nana (Table 2). The mean C-value of the triploid hybrids was intermediate between the two species, very close to the theoretical expectation. The morphological assessment of these plants, on the other hand, did not reveal a linear relationship among the three ploidy groups (Table 1, Figure 4).

4. Discussion

Genome sizes of Icelandic birch obtained in this study correlate perfectly with the ploidy level—the average 1C genome size of the tetraploid Betula pubescens (882 Mbp) is about twice the size of the genome of the diploid species B. nana (448 Mbp), and the triploid group has the average size of 666 Mbp, which is mid-way between the two species.

The discovery that the triploid hybrid genome is intermediate in its size has a great biological significance. Firstly, it confirms that the triploids are interspecific hybrids of these two coexisting birch species, most probably containing a complete genome from each parental species. Chromosomal mapping in Icelandic birch using the 18S-25S ribosomal markers [1, 30] has also shown that the triploid genome is a hybrid containing one set of chromosomes from each parental species. The triploid genome has five nonhomologous 18S-25S rDNA sites, two from B. nana and three from B. pubescens [1]. Secondly, the present study indicates that triploid birch in Iceland is probably of recent origin. As the genome size of the triploid group is an almost exact intermediate between the currently introgressed parental species, it must have been formed through recent hybridization. Our palynological study [31] has indeed shown that hybridization in Icelandic birch occurred after the first colonization of birch species in Iceland in the early Holocene and has been going on ever since. Thirdly, there seems to be no major genome/chromosome rearrangements accompanying the interspecific hybridization that could have resulted in significant deviation from the intermediate hybrid genome size observed in the present study. Although we need to measure genome size in a much larger sample size, these first results are statistically significant and hence meaningful. It is impressive that even in the face of extensive gene flow, based on morphological assessment, genome size has proven to be a stable character. There must have been some form of selection pressure that has maintained the species integrity despite introgression. Some elements in the arctic and subarctic environments may have a role in maintaining the genomic stability, by for example, selecting against aneuploidy.

Genome downsizing cannot be statistically proven in the polyploid genomes of these plants. This could be considered a rare situation given that genome downsizing is common in angiosperms [9], whereby polyploids are expected to have proportional larger C-values than their diploid progenitors. The present study shows 1.97 increase in genome size on average, from diploid to tetraploid levels, whereas the analysis of observed genome size changes in angiosperms [9] reveals at the most 1.2 increase for eudicots. The average Cx-values among different ploidies in the present study are also constant statistically. The reason for not seeing genome downsizing in our study could be that the diploid species B. nana is probably not the diploid progenitor of the tetraploid B. pubescens. The two species are taxonomically and geographically differentiated in a broad sense. In geographically overlapping zones, such as Iceland, the two species hybridize but this secondary contact does not give rise to a new species. Taxonomically, B. nana belongs to the subsection Nanae, whereas B. pubescens belongs to the subsection Albae with the other European tree birch species B. pendula [28, 32]. The ribosomal mapping on chromosomes of B. nana and B. pubescens [30] supports the two species having different genome origins. The phylogenetic studies using the nuclear ADH gene and DNA sequences of the chloroplast matK gene [33] as well as AFLP markers [34] indicate that the tetraploid species B. pubescens has a hybrid origin, that is, is an allopolyploid. However, no conclusive evidence has been found as to the identity of its diploid progenitors. Until the hybrid genomes of the tetraploid species B. pubescens have been correctly identified, it is not possible to test if genome downsizing has actually occurred in Betula.

The Plant DNA C-values Database [35] reports 1C-values of only two Betula species, an eastern N-American tree birch species B. populifolia with 0.2 pg [36] and B. pubescens (European tree birch species examined in the present study) with 0.75 pg [37], previously given as 0.7 pg [38]. Both species were measured using Feulgen densitometry after acetic alcohol fixation, without chromosome counts reported. For both species, a diploid chromosome number 2 = 28 is given in the C-values Database. However, all chromosome counts for B. pubescens according to Index to Plant Chromosome Numbers (IPCN, on 16 October 2009, http://mobot.mobot.org/W3T/Search/ipcn.html) and many earlier sources report the tetraploid number 2 = 56. Apart from this, the differing 1C-values of 0.7 or 0.75 pg for B. pubescens [37, 38] come from using a recalibrated standard value of Senecio vulgaris in Mowforth’s thesis [37] (compare Bennett and Smith [39, 40] and see footnotes b and i in the latter). Assuming that the values from [37, 38] actually refer to B. pubescens with 2 = 56, their 1C-value of 0.7 pg or 0.75 pg is 0.78- to 0.83-fold the present 1C-value of 0.902 pg for the tetraploid level. Both these lower Feulgen values are however compatible with our data, if possible tannin interference in their samples, which one could expect in birch after acetic alcohol fixation (see Greilhuber [29]) is assumed. The 1C-value of 0.2 pg for B. populifolia [36] appears too low for a diploid birch.

The congruence of FDM- with FCM-values in the present study and their unbiased correlation with ploidy levels agrees with the visual aspect of good, regularly stained formaldehyde-fixed Feulgen-stained nuclei, even though ample solidified tannins were present in many cells of the bud tissues. FDM is clearly a valid method for the purpose. Regarding FCM, the suitability of the chopping method for young decorticated and, in our case, dormant twigs is a useful addition to the flow cytometric practice. First of all, twigs of trees are, it seems, an organ of woody plants that are always available but hitherto never or rarely employed (see Supplemental Appendix in [41]). It is worth noting that in some cases stems of shrubby Fabaceae were used instead of leaves to improve histogram quality [42]. Secondly, this approach extends the application to the cold season, when leaves are not available and buds may be unsuitable because of secondary compounds such as resins. In birch, twig tissue has, in contrast to buds, no disturbing resins. The use of young and even dormant twigs for genome size measurement as reported in this study should be applicable for studying a broad range of tree species.

Acknowledgment

The authors thank I. J. Leitch for communicating details on methods used for Feulgen DNA measurement in Betula pubescens from the thesis of M. A. G. Mowforth.