Original Research ARTICLE
Polyploidy in the Conifer Genus Juniperus: An Unexpectedly High Rate
- 1Laboratoire Biodiversité et Génomique Fonctionnelle, Faculté des Sciences, Université Saint-Joseph, Campus Sciences et Technologies, Beirut, Lebanon
- 2Ecologie Systématique Evolution, Univ. Paris-Sud, CNRS, AgroParisTech, Université Paris-Saclay, Orsay, France
- 3Royal Botanic Gardens Kew, Richmond, United Kingdom
- 4Laboratori de Botànica, Facultat de Farmàcia, Universitat de Barcelona, Unitat Associada CSIC, Barcelona, Spain
- 5Biology Department, Sorbonne Université, Paris, France
- 6Biology Department, Baylor University, Waco, TX, United States
Recent research suggests that the frequency of polyploidy may have been underestimated in gymnosperms. One notable example is in the conifer genus Juniperus, where there are already a few reports of polyploids although data are still missing for most species. In this study, we evaluated the extent of polyploidy in Juniperus by conducting the first comprehensive screen across nearly all of the genus. Genome size data from fresh material, together with chromosome counts, were used to demonstrate that genome sizes estimated from dried material could be used as reliable proxies to uncover the extent of ploidy diversity across the genus. Our analysis revealed that 16 Juniperus taxa were polyploid, with tetraploids and one hexaploid being reported. Furthermore, by analyzing the genome size and chromosome data within a phylogenetic framework we provide the first evidence of possible lineage-specific polyploidizations within the genus. Genome downsizing following polyploidization is moderate, suggesting limited genome restructuring. This study highlights the importance of polyploidy in Juniperus, making it the first conifer genus and only the second genus in gymnosperms where polyploidy is frequent. In this sense, Juniperus represents an interesting model for investigating the genomic and ecological consequences of polyploidy in conifers.
Polyploidy or whole genome duplication (WGD) is the heritable condition of possessing more than two complete sets of chromosomes (Comai, 2005). Typically, polyploidy arises either as a result of genome duplication within a species (i.e., autopolyploidy), or from hybridization between two different species followed by chromosome doubling (allopolyploidy) (Stebbins, 1947; Comai, 2005). Most of our understanding of the consequences of polyploidy in plants has come from the study of angiosperms, where it has been shown that polyploidization generally causes a dramatic change in genomic structure, dynamics and expression, and cell organization (Tayalé and Parisod, 2013; Van de Peer et al., 2017; Wendel et al., 2018). Indeed, polyploidy is considered to have played a major role in angiosperm evolution (Blanc and Wolfe, 2004; Chen, 2007; Otto, 2007; Soltis and Soltis, 2009).
While polyploidy has been reported to occur across all major taxonomic land plant groups (Barker et al., 2016), it has been estimated to be very frequent in angiosperms with 50–80% of species being polyploid (Masterson, 1994; Otto and Whitton, 2000) and possibly all angiosperms contain at least one WGD in their ancestry (Van de Peer et al., 2017). In contrast, only 5% of all gymnosperms are reported to be polyploid based on chromosome counts (Khoshoo, 1959; Ahuja, 2005; Husband et al., 2013; Rice et al., 2015). Nevertheless, recent analyses of transcriptomic and genomic data (e.g., Li et al., 2015; Guan et al., 2016; Roodt et al., 2017) have suggested that the evolution of gymnosperms was accompanied by several ancient WGD events, including two within conifers, one at the base of Pinaceae (c. 200–342 million years ago) and one at the base of the cupressophytes (including Cupressaceae but excluding Araucaceae) (c. 210–275 million years ago). This highlights the importance of polyploidy in the very early evolution of conifers in contrast to the extreme rarity of this phenomenon among extant species [estimated to be 1.5% based on chromosome counts (Khoshoo, 1959; Husband et al., 2013; Rice et al., 2015)]. The one notable exception to the low frequency of polyploidy in extant gymnosperms is in Ephedra, which belongs to the non-coniferous lineage Gnetales. Here, polyploidy has been reported in over 65% of extant Ephedra species (Ickert-Bond et al., 2015). In this genus no evidence for any ancient WGDs has been detected in its ancestry (Li et al., 2015).
Conifers comprise the largest group of extant gymnosperms (Christenhusz et al., 2011), and from a phylogenetic perspective, they are divided into two major clades—(i) the Pinaceae and (ii) cupressophytes as they include Cupressaceae which is the most species-rich family (Lu et al., 2014; Ran et al., 2018). Within extant conifers, chromosome counts of all studied wild stands of all genera of Pinaceae are reported to be diploid (2n = 2x = 24) (Hizume, 1988; Murray, 2013) despite an exceptional genome size variation in some genera, such as Pinus L. (34.5–72.0 pg/2C) (Bogunic et al., 2003; Murray et al., 2012).
Similarly, in Cupressaceae, among ca. 91 species studied for their chromosome number to date (Hair, 1968; Murray, 2013), nearly all are diploid (2n = 2x = 22), with just three natural polyploids reported: Sequoia sempervirens is hexaploid with 2n = 6x = 66 (Ahuja and Neale, 2002; Scott et al., 2016), while Fitzroya cupressoides (Molina) I. M. Johnst. (alerce) and Juniperus thurifera L. are tetraploid with 2n = 4x = 44 (Hair, 1968; Romo et al., 2013; Vallès et al., 2015). It is also notable that within Juniperus, the study of just three species revealed each had polyploid cytotypes in some populations (Sax and Sax, 1933; Nagano et al., 2007). These findings raise the question of whether polyploidy may be common in this genus and hence whether it has played a more significant role in the evolution of Cupressaceae than previously recognized in gymnosperms as a whole, and in conifers in particular.
In this study, we focused on exploring the prevalence of polyploidy in wild populations of Juniperus. With 115 taxa (75 species with 40 varieties; Adams (2014), also see Table 1 for species and varieties), Juniperus is the most diverse genus in Cupressaceae and the second most diverse in all conifers after Pinus (Farjon, 2010; Romo et al., 2013). Juniperus has been shown to be a well-supported monophyletic genus (Mao et al., 2010; Adams and Schwarzbach, 2013; Adams, 2014), that can be divided into three monophyletic sections: (i) section Caryocedrus, with one species in the Mediterranean; (ii) sect. Juniperus, with 14 species, 12 in East Asia and the Mediterranean, and one with a circumboreal distribution (Juniperus communis L.) and one [J. jackii (Rehder) R. P. Adams] endemic to North America; and (iii) sect. Sabina, with ~60 species distributed in southwestern North America, Asia and the Mediterranean region, with outlier species in Africa and the Canary Islands. The few polyploids in wild populations noted above have all been reported to occur in species belonging to sect. Sabina. Both diploid and tetraploid cytotypes have been found in some populations of J. chinensis L. (Sax and Sax, 1933; Hall et al., 1973; Zonneveld, 2012) and in some populations of J. sabina L. (Siljak-Yakovlev et al., 2010; Farhat et al., 2019). Few sporadic triploid and tetraploid cytotypes have also been found in some ornamental cultivars. Juniperus thurifera is the only species reported to be exclusively tetraploid (2n = 4x = 44 and 40 pg/2C) (Romo et al., 2013; Vallès et al., 2015). More recently, Bou Dagher-Kharrat et al. (2013) showed that J. foetidissima Willd. had a very large genome (59.74 pg/2C), c. 3-fold larger than confirmed diploid Juniperus species which range from 19.02 to 26.40 pg/2C (Bennett and Leitch, 2012). The exceptional genome size of J. foetidissima, suggests this species may be hexaploid (Bou Dagher-Kharrat et al., 2013) but cytogenetic studies are needed to confirm this since genome size alone may be misleading as it can be highly variable between species of the same genus that have the same ploidy level (Ledig, 1998; Morse et al., 2009; Abdel Samad et al., 2014).
Table 1. List of the Juniperus taxa studied with indication of data collection, type of material, genome size data, and chromosome numbers.
Altogether, these observations suggest that Juniperus may have undergone an unusual evolutionary trajectory, involving polyploidization more frequently than encountered in other conifers. This paper takes a first step toward addressing these gaps in our data to fully understand the role that polyploidization has played in the evolutionary history of Juniperus. The objective was to assess variation in genome size across the whole genus and use these data as a proxy to estimate ploidy levels. Using classical cytogenetics approaches, we also determined the ploidy level of J. foetidissima, which has the biggest genome in this genus. Finally, we used phylogenetically-informed trait evolution modeling approaches to reconstruct ancestral genome sizes for the three main clades of Juniperus and identify the occurrence of polyploidization events during the evolution of Juniperus.
Materials and Methods
The origins of the studied accessions are presented in Table 1. We used Robert P. Adams's worldwide collection of Juniperus leaf material, dried in silica gel and kept frozen at −20°C. This material has been stored for years (the oldest sample was collected in 1985). To address its suitability for genome size analysis and ploidy screening, we carried out measurements on both dry and fresh material for a sub-sample of 12 species which were selected to cover as much of the genus diversity at the taxonomic (representatives of sections Juniperus and Sabina), morphological (needles-like and scale leaves) and cytogenetic (species with different ploidy levels) levels. Fresh leave material was obtained from plants growing in the living collections of the Royal Botanic Gardens, Kew, UK.
Genome Size Assessments by Flow Cytometry
Nuclear DNA contents of about 3,000 stained nuclei were estimated for each sample with a CyFlowSL Partec flow cytometer (Partec GmbH) following the one-step protocol of Doležel et al. (2007) with minor modifications as described in Clark et al. (2016). We selected Allium cepa L., 2C = 34.89 pg (Doležel et al., 1998; Clark et al., 2016) and the “CyStain PI Absolute P kit” buffer (Sysmex UK) as the most appropriate internal calibration standard and nuclei isolation buffer for ploidy screening in Juniperus.
We compiled published Juniperus chromosome numbers from the Chromosome Counts Database (CCDB; Rice et al., 2015). New chromosome counts were made for J. foetidissima and J. excelsa using 3 years old plants cultivated from seed of natural origin (from Turkey), and following Vallès et al. (2015) for protoplast preparation and Chromomycin A3 (CMA, Serva) staining.
Analyses of Genome Size and Chromosome Number Evolution
Trait evolution was modeled on the phylogenic tree of Adams (2014), pruned to the set of species and varieties with genome size data and made ultrametric with R v.3.2.2 (Team, 2016). However, five taxa with genome size estimates were not represented in the phylogeny and so they were discarded from these analyses [Juniperus communis var. kelleyi R. P. Adams, J. deltoides var. spilinanus (Yalt., Elicin and Terz.) Terz, J. durangensis var. topiensis R. P. Adams and S. Gonzalez, J. poblana var. decurrens R. P. Adams, J. semiglobosa var. talassica (Lipsky) Silba)]. The inference of ancestral genome size values was based on monoploid GS (1Cx-values) sensu Greilhuber et al. (2005). Ancestral 1Cx-values were reconstructed under ML using the “fastAnc” command and mapped onto the phylogeny with the “contMap” command of the Phytools package of R (Revell, 2012).
We used ChromEvol v.2 (Glick and Mayrose, 2014) to infer ancestral haploid (n) chromosome numbers in Juniperus. This program implements a series of likelihood models to infer duplication events, chromosome gains/losses and demi-duplications at ancestral nodes. The model that best fitted the data set was chosen under the Akaike information criterion (AIC) using default parameters.
Genome Size Diversity
Genome sizes were assessed for 111 Juniperus species and varieties (Table 1), representing 96.5% of taxonomic diversity. Low differences were found between values obtained with dried and fresh material for the 12 species analyzed using both types of leaf material. Differences varied around zero with six positive (minimum = 0.6%, maximum = 9.8% and mean difference = 3.1%) and six negative percentages (minimum = −0.42%, maximum = −3.16% and mean difference = −2.15%). Overall, the genome size estimates for Juniperus ranged 3.2-fold (from 21.81 to 70.58 pg/2C) but they were seen to be distributed into three non-overlapping classes (Figure 1), class A: 21.81–30.3 pg/2C, B: 46.29–50.7 pg/2C, and C: 70.58 pg/2C.
Figure 1. Genome size (2C-values, represented by black dots) classes in Juniperus and their unequivocal relationship with the chromosome number. Class A represents the range of genome sizes for all diploid species confirmed by published chromosome numbers. Class B represents the range of genome sizes for all tetraploid species confirmed by published chromosome number. Class C represents the genome size of the only hexaploid species so far reported (i.e., J. foetidissima). (A) Chromosomes of the diploid J. excelsa (our data); (B) chromosomes of the tetraploid J. thurifera (reproduced from Vallès et al., 2015), and (C) Chromosomes of J. foetidissima confirming its hexaploid status (our data). Monoploid genome size (1Cx-values, represented by yellow dots) of the three sections were also illustrated.
Ploidy Levels Inferred From Genome Size Data
We gathered chromosome number data from the CCDB for 41 Juniperus species and varieties (Table 1). In addition, we made the first chromosome counts for J. excels—a diploid with 2n = 22, and J. foetidissima—a hexaploid with 2n = 66 (Figures 1A,C, respectively). Ploidy levels based on chromosome numbers agreed with those inferred from genome size for all but two taxa, suggesting a strong correlation between genome size, ploidy level and chromosome number. Genome size values of class A corresponded to diploids with 2n = 2x = 22, class B to tetraploids with 2n = 4x = 44 and class C to hexaploids with 2n = 6x = 66 (Table 1; Figure 1). The two exceptions were J. seravschanica Kom. and J. chinensis var. sargentii A. Henry, which were both reported to be diploid in the CCDB but had genome size estimates indicating the samples analyzed here were tetraploid. We thus considered these taxa to have two cytotypes, as previously established for J. chinensis and J. sabina (Table 1).
Evolution of Chromosome Numbers
The best-fitting model in ChromEvol to explain the evolution of chromosome numbers in Juniperus was the CONST_RATE model (Supplementary Table S1), suggesting that polyploidy is the predominant mode of chromosome evolution in Juniperus. The ancestor of the whole genus was inferred to be diploid, with n = 11. It is noted that the polyploids were exclusively restricted to sect. Sabina (Figure 2). Three lineage-specific polyploidization events leading to tetraploidy were detected in the ancestors of the clades giving rise to (i) J. recurva, J. rushforthiana, J. indica, (ii) J. preswalskii, J. tibetica, J. morrisonicola, J. squamata, and (iii) J. thurifera, J. foetidissima (Figure 2). A further gain of 22 chromosomes was inferred in the lineage giving rise to the hexaploid J. foetidissima. Six species-specific or within-species polyploidization events (i.e., cytotypes) were found in J. coxii, J. sevaschanica, J. chinensis, J. chinensis var. procumbens, J. chinensis var. sargentii and J. sabina, all of which contained both diploid and tetraploid individuals (Figure 2).
Figure 2. Ancestral state reconstruction of genome size (1Cx/pg) and chromosome number (n) on a phylogenetic tree of Juniperus reconstructed using Bayesian approaches (Adams, 2014). An illustration of the leaf shape for each section is represent by: (a) J. drupacea (sect. Caryocedrus); (b) J. communis (sect. Juniperus); (c) J. excelsa (sect. Sabina).
Evolution of Genome Size
Beside the genome size variation explained by chromosome number difference, a small variation at the 1Cx-level was detected between ploidy levels. In addition, the distribution of 1Cx-values across Juniperus presented in Figures 1, 2 showed an ancestral genome size of 12.37 pg for the whole genus and overall larger values in species belonging to sect. Sabina (mean 1Cx 12.7 pg, ancestral 1Cx 12.64 pg) compared with those of sect. Caryocedrus (mean 1Cx 11.74 pg, ancestral 1Cx 12.15 pg) and sect. Juniperus (mean 1Cx 11.38 pg, ancestral 1Cx 11.59 pg). Nevertheless, decreases in 1Cx-values were observed in several taxa from sect. Sabina, including some –but not all– polyploids. Polyploid taxa showed limited 1Cx variation relative to the value inferred for their most recent ancestors, with a maximum 1Cx downsizing of 5.70% for J. squamata var. meyeri, and a maximum 1Cx upsizing of 1.71% in J. rushfortiana (Supplementary Table S2).
Reliability of Genome Size Estimates From Desiccated Leaf Material of Juniperus
Over the years considerable attention has focused on exploring the suitability of dried plant material for genome size and ploidy level analysis, especially given the challenges of collecting and analyzing fresh material from plants growing in remote locations. Dried material has certainly shown to be suitable for ploidy level analysis in many vascular plants (Suda and Trávníček, 2006; Schönswetter et al., 2007; Suda et al., 2007; Popp et al., 2008; Krejčíková et al., 2013; Wang and Yang, 2016). Nevertheless, the quality of data generated by flow cytometry using dried material has been shown to differ between species, buffers (Bainard et al., 2011) and type of desiccation used (Šmarda et al., 2005; Šmarda and Stančík, 2006; Suda and Trávníček, 2006) and it is now generally accepted that while desiccated material is suitable for ploidy level analysis, it is usually not reliable enough for accurate genome size estimations.
In contrast to these previous studies, our analyses of Juniperus showed that leaves dried in silica gel and stored continuously at −20°C are suitable for genome size estimations using flow cytometry, giving reasonable data quality (i.e., mean %CV = 3.9, S.D. = 0.96). This was supported by comparisons of 2C-values estimated for the same species from dried and fresh material where low differences between the two variances were found in the 12 species analyzed. We are thus confident that the genome size data generated from the desiccated material analyzed here are reliable and hence suitable for exploring genome size [but there might be a slight shift in “absolute” genome sizes (9.8% at maximum)] and ploidy diversity and evolution across Juniperus. Our results broadly agree with Bainard et al. (2011) who found that leaves desiccated immediately in the field using silica gel, was one of the most promising conservation methods, yielding reasonable quality flow cytometry peaks for some species.
Variability in Genome Size and Polyploidy in Juniperus
This study showed that junipers are characterized by possessing large genomes (mean genome size for diploid taxa = 25 pg/2C) with extensive variation between species (ranging 3.2-fold from 21.81 to 71.32 pg/2C). This large variation perfectly correspond to known ploidy levels (2x – 6x), while the variation in 1Cx is only 1.38-fold. The data considerably extend our knowledge of genome sizes in Juniperus which was previously based on data for just 19 species (Bennett and Leitch, 2012). They also show Juniperus now has the largest range in genome size so far reported for any conifer genus.
There are three main mechanisms which can lead to variation in genome size; (i) rapid loss or expansion of transposable and/or other repetitive elements, (ii) loss or gain of chromosomes (aneuploidy and dysploidy), and (iii) polyploidization, possibly followed by genome downsizing (Ramsey and Schemske, 1998; Leitch and Bennett, 2004; Greilhuber et al., 2005; Morse et al., 2009). While in Pinus the high variability in genome size (34.50–72.00 pg/2C; Murray et al., 2012) has been shown to be mainly driven by variation in copy numbers of repeats, such as retrotransposable elements (Morse et al., 2009; Kovach et al., 2010; Nystedt et al., 2013), in Juniperus, our data indicate that most of the variation in genome size is due to variation in ploidy levels. This does not exclude the occurrence of limited genome size variation within each ploidy level, but based on the data presented, it is relatively small, ranging just 1.4-fold in diploids (95 taxa) and 1.1-fold in tetraploids (15 taxa). The source of this variation is still unclear but likely to represent variation in repeat content since, to date, there have been no reports of aneuploidy in the genus (Murray, 2013).
Among the 111 taxa analyzed, just two (J. chinensis var. sargentii and J. seravschanica) showed a discrepancy between the chromosome number reported in the CCDB and the ploidy level estimated from the genome size data obtained here. This could be due to a technical error, such as misidentification of the species used for counting chromosomes and such an explanation is possible for J. seravschanica, where the synonym taxa J. macropoda Boiss. has been used to determine the ploidy level (Rice et al., 2015). Nevertheless, these exceptions could also be explained by the existence of intra-specific variability in ploidy levels (= cytotype diversity), a well-documented phenomenon encountered in many land plant lineages, especially in angiosperms and ferns (Husband et al., 2013). In contrast, cytotype diversity is rarely reported in gymnosperms, with Ephedra being the only genus where it occurs extensively (>50% of species have >1 cytotype—Ickert-Bond et al., 2015). Prior to this study, natural intraspecific variation in ploidy level in Juniperus had only been reported in a few species including in J. chinensis (2x, 4x) (Sax and Sax, 1933; Hall et al., 1973) and J. sabina (2x, 4x) (Siljak-Yakovlev et al., 2010; Farhat et al., 2019).
In view of these previous studies, the results presented here are striking—revealing a much higher frequency of polyploidy in Juniperus than hitherto detected, with 15% of taxa being tetraploid, and the discovery of an hexaploid (J. foetidissima), which is only the second hexaploid to be found in conifers. In addition, the use of ChromEvol to infer the evolution of chromosome numbers across the phylogeny of Juniperus suggests that there have been an unexpectedly high number of polyploidization events throughout its evolutionary history compared with other gymnosperm lineages (except Ephedra). Such a result suggests that mechanisms that promote polyploidization and/or the evolutionary success of polyploid species have occurred at a much higher frequency in Juniperus than in other conifers, and even in gymnosperms in general, apart from Ephedra. It is also worth noting that only one individual was analyzed for most taxa in this study. It is therefore possible that our data underestimate the importance of polyploidization in Juniperus as additional intraspecific ploidy diversity may well be uncovered when more individuals are analyzed, as already seen in J. sabina and J. chinensis.
Genome Size Evolution and Ploidy Levels of Juniper Ancestors
Studies exploring the evolution of genome size diversity across different land plant groups, have uncovered contrasting dynamics in genome size fluctuations throughout their evolution (Bainard and Villarreal, 2013; Clark et al., 2016; Soltis et al., 2018). Now that genome size data are available for almost every recognized taxa of Juniperus and that ploidy levels can be inferred given the robust relationship with genome size (Figure 1), the reconstruction of the ancestral genome size within this genus and inferred ancestral ploidy level is highly instructive. Indeed, apart from Pinus (Grotkopp et al., 2004), our study is the first to reconstruct ancestral genome size within a species-rich genus for any gymnosperm. Our analysis revealed that the ancestral ploidy level for Juniperus was diploid with an estimated genome size of 12.37 pg/1C, which fits within the range of 9–12.38 pg/1C inferred by Burleigh et al. (2012), based on a sampling including only two Juniperus species amongst 165 gymnosperm species.
Within the genus, we found evidence suggesting that fluctuations in genome size, both upsizing and downsizing, independent of polyploidy, have taken place during evolution, as also found in Pinus (Grotkopp et al., 2004) and across other gymnosperm lineages as well (Burleigh et al., 2012). However, while, in most other gymnosperm genera the shifts in genome size are likely to be driven by changes in the abundance of repetitive DNA (Nystedt et al., 2013; De La Torre et al., 2014), in Juniperus the large shifts in genome size are associated with polyploidization events, with a minimum of 10 such events predicted from our analyses (Figure 2). Whether the occurrence and frequency of polyploidy, which was seen to be restricted to sect. Sabina, contributes to the higher number of species in this section (c. 60 species) compared with the other two sections of Juniperus (sect. Juniperus = c. 13 species, sect. Caryocedrus = one species) is unclear, although previous studies pointing to higher diversification rates in some angiosperm lineages following polyploidy suggest this is possible (Wood et al., 2009; Landis et al., 2018).
Concerning the origin of the hexaploid, J. foetidissima, there are several possible pathways. It could have arisen from a triploid ancestor following one step. If so, then there are two possible routes; (i) fertilization between two unreduced triploid gametes of a triploid ancestor, or (ii) somatic doubling of a triploid, giving rise directly to the hexaploid. Alternatively, it could have arisen following two WGD events (two steps) as envisaged for the hexaploid Sequoia sempervirens (Scott et al., 2016). The first step being a WGD event either via autopolyploidy or allopolyploidy leading to the formation of a tetraploid with n = 2x, followed by hybridization with a diploid (n = x) leading to a triploid. The second step involves a WGD giving rise to a hexaploid. The reports of sporadic triploid Juniperus individuals indicate that triploids can indeed form (Hall et al., 1973). However, yet another possibility is that the origin of J. foetidissima does not involve a triploid, but instead arose from hybridization between an unreduced gamete from a tetraploid (4x) with either (a) a reduced gamete from another tetraploid (2x) or (b) an unreduced gamete from a diploid (2x). Currently, there is no information about the genomic makeup of J. foetidissima to know whether it is an auto- or allo-polyploid, or its mode of origin.
Why Is Polyploidy More Common in Juniperus Than Other Conifers?
The success of hexaploid Sequoia sempervirens and polyploid Ephedra species (4x – 8x), has been partially attributed to their capacity for vegetative propagation (Scott et al., 2016; Wu et al., 2016) and this may also contribute to the survival of polyploid Juniperus species as there is evidence that they too have the capacity for vegetative propagation [e.g., in J. sabina and J. communis (Houle and Babeux, 1994; Ronnenberg, 2005; Wesche et al., 2005; Tylkowski, 2010)]. Furthermore, the extreme longevity has been suggested to be another factor contributing to the success of polyploidy in S. sempervirens (Scott et al., 2016), and since Juniperus has been classified as long-lived (Ward, 1982; Gauquelin et al., 2012) this may also help the survival of polyploids, enabling them to become established.
Here we propose a novel hypothesis that may also contribute to higher frequency of polyploidy revealed in Juniperus—this is the high frequency of sympatry between juniper species. In contrast to most of the conifers, the geographical ranges of Juniperus species overlap considerably which opens up lots of opportunities for natural hybridization between species. For example, in Spain, hybrids between J. thurifera × J. sabina and J. thurifera × J. phoenicea and J. sabina × J. phoenicea in sympatry have been described (Rojo and Díaz, 2006, 2009; Rojo and Uribe-Echebarría, 2008). More recently, Adams et al. (2016) suggested that an ancient hybridization between J. thurifera and J. sabina gave rise to J. sabina var. balkanensis. Juniper hybrids are also common in North America between closely related species in areas of sympatry [e.g., between J. virginiana L. and J. horizontalis Moench, J. osteosperma Hook and J. occidentalis Torr. Little, J. virginiana var. silicicola, and J. bermudiana (Vasek, 1966; Palma-Otal et al., 1983; Adams and Kistler, 1991; Adams and Wingate, 2008; Adams, 2014)].
Even though the sympatry is a sine qua non condition for natural hybridization, there are few cases of conifers occurring in sympatry that do hybridize without giving rise to polyploids: e.g., Pinus taeda and P. echinata (Edwards-Burke et al., 1997). Furthermore, induced hybridization like for Cedrus species (Fady et al., 2003) produced only homoploids. Cases of unreduced gamete production were documented in Cupressaceae (Pichot and El Maâtaoui, 2000) and Ephedraceae (Wu et al., 2016). This ability to produce unreduced gametes may be the explanation for polyploidisation in Juniperus.
On the other hand, the genomic shock arising from hybridization can often be ameliorated by WGD and subsequent diploidization as it was shown in angiosperms (Hegarty et al., 2006). Given the high frequency of hybrid formation in Juniperus, and assuming that similar levels of genomic shock following hybridization also occur here, as in angiosperms, then it is possible to envisage that polyploidy may offer one potential solution to these genomic challenges, tipping the balance toward their survival in the wild. Clearly, studies are now needed at the molecular level to provide insights into whether our understanding of the genomic consequences of hybridization and polyploidization in angiosperms is also applicable to the growing list of gymnosperm polyploids.
Polyploidy or whole genome duplication is rare in conifers. The lack of studies on polyploidy within Juniperus prompted the present study, in which the ploidy level of 96.5% of the genus was screened in order to explore the extent of polyploidy across the genus. Silica gel-dried leaves of Juniperus were found to be highly suitable for genome size measurements using flow cytometry. This study uncovered a relatively high number of polyploidization events (at least 10) in Juniperus, compared to other conifers, and revealed that at least 15% of Juniperus taxa are tetraploids. In addition, we used both chromosome and genome size data to validate the presence of the only hexaploid in Juniperus (J. foetidissima) so far reported, and only the second hexaploid found in conifers (after Sequoia sempervirens). An analysis of the phylogenetic distribution of polyploids across Juniperus showed they were restricted to sect. Sabina and that three clades are exclusively made of polyploids (one including the hexaploid J. foetidissima), providing the first evidence of possible lineage-specific polyploidizations in the genus.
Overall, it seems clear that Juniperus is exceptional within conifers, and represents a second genus within gymnosperms where polyploidy is common. We propose that Juniperus should be considered to be a highly relevant model for studying polyploidization mechanisms and pathways in conifers, and comparisons with Ephedra will provide a comprehensive understanding of the evolutionary dynamics and consequences of polyploidy in gymnosperms.
MB designed the study. RA provided the Juniperus material. PF and OH carried out the flow cytometry measurements and analyzed the data. PF and SS-Y determined the chromosome numbers. PF wrote a first draft of the manuscript that was further critically reviewed by MB, RA, OH, SS-Y, IL, TR.
The authors thank the National Council for Scientific Research grant number CNRS-FS90—Lebanon, the Saint Joseph University Research Council (CR-USJ) FS-111 for supporting financially this work.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank the Royal Botanic Gardens Kew, London, UK for providing access to the flow cytometry facilities and living collections.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00676/full#supplementary-material
Abdel Samad, F., Baumel, A., Juin, M., Pavon, D., Siljak-Yakovlev, S., Médail, F., et al. (2014). Phylogenetic diversity and genome sizes of Astragalus (Fabaceae) in the Lebanon biogeographical crossroad. Plant Syst. Evol. 300, 819–830. doi: 10.1007/s00606-013-0921-8
Adams, R., Schwarzbach, A., and Tashev, A. (2016). Chloroplast capture by a new variety, Juniperus sabina var. balkanensis RP Adams and AN Tashev, from the Balkan peninsula: a putative stabilized relictual hybrid between J. sabina and ancestral J. thurifera. Phytologia 98, 100–111.
Ahuja, M. R., and Neale, D. B. (2002). Origins of polyploidy in coast redwood (Sequoia sempervirens (D. don) Endl. and relationship of coast redwood to other genera of Taxodiaceae. Silvae Genet. 51, 93–99.
Bainard, J. D., Husband, B. C, Baldwin, S., Fazekas, A., Gregory, T., Newmaster, S., et al. (2011). The effects of rapid desiccation on estimates of plant genome size. Chromosome Res. 19, 825–842. doi: 10.1007/s10577-011-9232-5
Bennett, M., and Leitch, I. (2012). Plant DNA C-Values Database (release 6.0, Dec. 2012). Available online at: http://www.kew.org/cvalues/
Bogunic, F., Muratovic, E., Brown, S., and Siljak-Yakovlev, S. (2003). Genome size and base composition of five Pinus species from the Balkan region. Plant Cell Rep. 22, 59–63. doi: 10.1007/s00299-003-0653-2
Bou Dagher-Kharrat, M., Abdel-Samad, N., Douaihy, B., Bourge, M., Fridlender, A., Siljak-Yakovlev, S., et al. (2013). Nuclear DNA C-values for biodiversity screening: case of the Lebanese flora. Plant Biosyst. 147, 1228–1237. doi: 10.1080/11263504.2013.861530
Burleigh, J. G., Barbazuk, W. B., Davis, J. M., Morse, A. M., and Soltis, P. S. (2012). Exploring diversification and genome size evolution in extant gymnosperms through phylogenetic synthesis. J. Bot. 2012:292857. doi: 10.1155/2012/292857
Chen, Z. J. (2007). Genetic and epigenetic mechanisms for gene expression and phenotypic variation in plant polyploids. Annu. Rev. Plant Biol. 58, 377–406. doi: 10.1146/annurev.arplant.58.032806.103835
Christenhusz, M. J., Reveal, J. L., Farjon, A., Gardner, M. F., Mill, R. R., and Chase, M. W. (2011). A new classification and linear sequence of extant gymnosperms. Phytotaxa 19, 55–70. doi: 10.11646/phytotaxa.19.1.3
Clark, J., Hidalgo, O., Pellicer, J., Liu, H., Marquardt, J., Robert, Y., et al. (2016). Genome evolution of ferns: evidence for relative stasis of genome size across the fern phylogeny. New Phytol. 210, 1072–1082. doi: 10.1111/nph.13833
Doležel, J., Greilhuber, J., Lucretti, S., Meister, A., Lysák, M., Nardi, L., et al. (1998). Plant genome size estimation by flow cytometry: inter-laboratory comparison. Ann. Bot. 82(Suppl_1), 17–26. doi: 10.1093/oxfordjournals.aob.a010312
Edwards-Burke, M. A., Hamrick, J. L., and Price, R. A. (1997). Frequency and direction of hybridization in sympatric populations of Pinus taeda and P. echinata (Pinaceae). 84, 879–886. doi: 10.2307/2446277
Fady, B., Lefèvre, F., Reynaud, M., Vendramin, G. G., Bou Dagher-Kharrat, M., Anzidei, M., et al. (2003). Gene flow among different taxonomic units: evidence from nuclear and cytoplasmic markers in Cedrus plantation forests. Theor. Appl. Genet. 107, 1132–1138. doi: 10.1007/s00122-003-1323-z
Farhat, P., Siljak-Yakovlev, S., Robert, A., Magda, B., and Robert, T. (2019). Genome size variation and polyploidy in the geographical range of Juniperus sabina L. (Cupressaceae). Bot. Lett. 68, 92–96. doi: 10.1080/00087114.2015.1024546
Gauquelin, T., Chondroyannis, P., Boukhdoud, N., Bouyssou, M., Brunel, C., Danneyrolles, V., et al. (2012). Le Genévrier thurifère, espèce partagée au Nord et au Sud de la Méditerranée. Forêt Méditerranéenne 33, 227–240.
Glick, L., and Mayrose, I. (2014). ChromEvol: assessing the pattern of chromosome number evolution and the inference of polyploidy along a phylogeny. Mol. Biol. Evol. 31, 1914–1922. doi: 10.1093/molbev/msu122
Greilhuber, J., DoleŽel, J., Lysák, M. A., and Bennett, M. D. (2005). The origin, evolution and proposed stabilization of the terms ‘genome size'and ‘C-value'to describe nuclear DNA contents. Ann. Bot. 95, 255–260. doi: 10.1093/aob/mci019
Grotkopp, E., Rejmánek, M., Sanderson, M. J., and Rost, T. L. (2004). Evolution of genome size in pines (Pinus) and its life-history correlates: supertree analyses. Evolution 58, 1705–1729. doi: 10.1111/j.0014-3820.2004.tb00456.x
Hegarty, M. J., Barker, G. L., Wilson, I. D., Abbott, R. J., Edwards, K. J., and Hiscock, S. J. (2006). Transcriptome shock after interspecific hybridization in Senecio is ameliorated by genome duplication. Curr. Biol. 16, 1652–1659. doi: 10.1016/j.cub.2006.06.071
Houle, G., and Babeux, P. (1994). Variations in rooting ability of cuttings and in seed characteristics of five populations of Juniperus communis var. depressa from subarctic Quebec. Can. J. Bot. 72, 493–498. doi: 10.1139/b94-066
Husband, B. C., Baldwin, S. J., and Suda, J. (2013). “The incidence of polyploidy in natural plant populations: major patterns and evolutionary processes,” in Plant Genome Diversity Volume 2: Physical Structure, Behaviour and Evolution of Plant Genomes, eds. J. Greilhuber, J. Dolezel, and J. F. Wendel (Vienna: Springer Vienna), 255–276. doi: 10.1007/978-3-7091-1160-4_16
Ickert-Bond, S., Pellicer, J., Souza, A., Metzgar, J., and Leitch, I. J. (2015). “Ephedra-the gymnosperm genus with the largest and most diverse genome sizes driven by a high frequency of recently derived polyploidy taxa and a lack of genome downsizing,” in Annual Meeting of the Botanical Society of America, Botany 2015, Abstract ID 862 (Edmonton).
Kovach, A., Wegrzyn, J. L., Parra, G., Holt, C., Bruening, G. E., Loopstra, C. A., et al. (2010). The Pinus taeda genome is characterized by diverse and highly diverged repetitive sequences. BMC Genomics 11, 1–14. doi: 10.1186/1471-2164-11-420
Krejčíková, J., Sudová, R., Lučanová, M., Trávníček, P., Urfus, T., Vít, P., et al. (2013). High ploidy diversity and distinct patterns of cytotype distribution in a widespread species of Oxalis in the Greater Cape Floristic Region. Ann. Bot. 111, 641–649. doi: 10.1093/aob/mct030
Landis, J. B., Soltis, D., Li, Z., Marx, H., Barker, M., Tank, D., et al. (2018). Impact of whole-genome duplication events on diversification rates in angiosperms. Am. J. Bot. 105, 348–363. doi: 10.1002/ajb2.1060
Li, Z., Baniaga, A., Sessa, E., Scascitelli, M., Graham, S., Rieseberg, L., et al. (2015). Early genome duplications in conifers and other seed plants. Sci. Adv. 1:e1501084. doi: 10.1126/sciadv.1501084
Lu, Y., Ran, J.-H., Guo, D.-M., Yang, Z.-Y., and Wang, X.-Q. (2014). Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes. PLoS ONE 9:e107679. doi: 10.1371/journal.pone.0107679
Mao, K., Hao, G., Liu, J., Adams, R., and Milne, R. (2010). Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals. New Phytol. 188, 254–272. doi: 10.1111/j.1469-8137.2010.03351.x
Morse, A. M., Peterson, D. G., Islam-Faridi, M. N., Smith, K. E., Magbanua, Z., Garcia, S. A., et al. (2009). Evolution of genome size and complexity in Pinus. PLoS ONE 4:e4332. doi: 10.1371/journal.pone.0004332
Murray, B., Leitch, I. J., and Bennett, M. D. (2012). Gymnosperm DNA C-Values Database (release 5.0, Dec. 2012). Available online at: http://www.kew.org/cvalues/
Murray, B. G. (2013). “Karyotype variation and evolution in gymnosperms,” in Plant Genome Diversity, Vol. 2, eds I. J. Greilhuber, J. Doležel, and J. Wendel (Vienna: Springer-Verlag), 231–243. doi: 10.1007/978-3-7091-1160-4_14
Nagano, K., Matoba, H., Yonemura, K., Matsuda, Y., Murata, T., and Hoshi, Y. (2007). Karyotype analysis of three Juniperus species using fluorescence in situ hybridization (FISH) with two ribosomal RNA genes. Cytologia 72, 37–42. doi: 10.1508/cytologia.72.37
Nystedt, B., Street, N. R., Wetterbom, A., Zuccolo, A., Lin, Y.-C., Scofield, D. G., et al. (2013). The Norway spruce genome sequence and conifer genome evolution. Nature 497:579. doi: 10.1038/nature12211
Palma-Otal, M., Moore, W., Adams, R., and Joswiak, G. (1983). Morphological, chemical, and biogeographical analyses of a hybrid zone involving Juniperus virginiana and J. horizontalis in Wisconsin. Can. J. Bot. 61, 2733–2746. doi: 10.1139/b83-301
Popp, M., Gizaw, A., Nemomissa, S., Suda, J., and Brochmann, C. (2008). Colonization and diversification in the African ‘sky islands' by Eurasian Lychnis L. (Caryophyllaceae). J. Biogeogr. 35, 1016–1029. doi: 10.1111/j.1365-2699.2008.01902.x
Ran, J.-H., Shen, T.-T., Wang, M.-M., and Wang, X.-Q. (2018). Phylogenomics resolves the deep phylogeny of seed plants and indicates partial convergent or homoplastic evolution between Gnetales and angiosperms. R. Soc. 285:20181012. doi: 10.1098/rspb.2018.1012
Rice, A., Glick, L., Abadi, S., Einhorn, M., Kopelman, N. M., Salman-Minkov, A., et al. (2015). The Chromosome Counts Database (CCDB)–a community resource of plant chromosome numbers. New Phytol. 206, 19–26. doi: 10.1111/nph.13191
Romo, A., Hidalgo, O., Boratynski, A., Sobierajska, K., Jasinska, A. K., Vallès, J., et al. (2013). Genome size and ploidy levels in highly fragmented habitats: the case of western Mediterranean Juniperus (Cupressaceae) with special emphasis on J. thurifera L. Tree Genet. Genomes 9, 587–599. doi: 10.1007/s11295-012-0581-9
Ronnenberg, K. (2005). Reproductive ecology of two common woody species, Juniperus sabina and Artemisia santolinifolia, in mountain steppes of southern Mongolia. Erforsch. Biol. Ress. Mongolei (Halle/Saale) 9, 207–223.
Roodt, D., Lohaus, R., Sterck, L., Swanepoel, R., Van de Peer, Y., and Mizrachi, E. (2017). Evidence for an ancient whole genome duplication in the cycad lineage. PLoS ONE 12:e0184454. doi: 10.1371/journal.pone.0184454
Schönswetter, P., Suda, J., Popp, M., Weiss-Schneeweiss, H., and Brochmann, C. (2007). Circumpolar phylogeography of Juncus biglumis (Juncaceae) inferred from AFLP fingerprints, cpDNA sequences, nuclear DNA content and chromosome numbers. Mol. Phylogenet. Evol. 42, 92–103. doi: 10.1016/j.ympev.2006.06.016
Scott, A. D., Stenz, N. W., Ingvarsson, P. K., and Baum, D. A. (2016). Whole genome duplication in coast redwood (Sequoia sempervirens) and its implications for explaining the rarity of polyploidy in conifers. New Phytol. 211, 186–193. doi: 10.1111/nph.13930
Siljak-Yakovlev, S., Pustahija, F., Šolić, E., Bogunić, F., Muratović, E., Bašić, N., et al. (2010). Towards a genome size and chromosome number database of Balkan flora: C-values in 343 taxa with novel values for 242. Adv. Sci. Lett. 3, 190–213. doi: 10.1166/asl.2010.1115
Šmarda, P., and Stančík, D. (2006). Ploidy level variability in South American fescues (Festuca L., Poaceae): use of flow cytometry in up to 5 1/2-year-old caryopses and herbarium specimens. Plant Biol. 8, 73–80. doi: 10.1055/s-2005-872821
Soltis, D., Soltis, P., Endress, P., Chase, M. W., Manchester, S., Judd, W., et al. (2018). Phylogeny and Evolution of the Angiosperms: Revised and Updated Edition. Chicago, IL: University of Chicago Press.
Suda, J., and Trávníček, P. (2006). Reliable DNA ploidy determination in dehydrated tissues of vascular plants by DAPI flow cytometry—new prospects for plant research. Cytometry Part A 69, 273–280. doi: 10.1002/cyto.a.20253
Suda, J., Weiss-Schneeweiss, H., Tribsch, A., Schneeweiss, G. M., Trávníček, P., and Schönswetter, P. (2007). Complex distribution patterns of di-, tetra-, and hexaploid cytotypes in the European high mountain plant Senecio carniolicus (Asteraceae). Am. J. Bot. 94, 1391–1401. doi: 10.3732/ajb.94.8.1391
Team, R. C. (2016). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/
Vallès, J., Garnatje, T., Robin, O., and Siljak-Yakovlev, S. (2015). Molecular cytogenetic studies in western Mediterranean Juniperus (Cupressaceae): a constant model of GC-rich chromosomal regions and rDNA loci with evidences for paleopolyploidy. Tree Genet. Genomes 11, 1–8. doi: 10.1007/s11295-015-0860-3
Wendel, J. F., Lisch, D., Hu, G., and Mason, A. (2018). The long and short of doubling down: polyploidy, epigenetics, and the temporal dynamics of genome fractionation. Curr. Opin. Genet. Dev. 49, 1–7. doi: 10.1016/j.gde.2018.01.004
Wesche, K., Ronnenberg, K., and Hensen, I. (2005). Lack of sexual reproduction within mountain steppe populations of the clonal shrub Juniperus sabina L. in semi-arid southern Mongolia. J. Arid Environ. 63, 390–405. doi: 10.1016/j.jaridenv.2005.03.014
Wood, T. E., Takebayashi, N., Barker, M. S., Mayrose, I., Greenspoon, P. B., and Rieseberg, L. H. (2009). The frequency of polyploid speciation in vascular plants. Proc. Natl. Acad. Sci. U.S.A. 106, 13875–13879. doi: 10.1073/pnas.0811575106
Wu, H., Ma, Z., Wang, M. M., Qin, A.-L., Ran, J. H., and Wang, X. Q. (2016). A high frequency of allopolyploid speciation in the gymnospermous genus Ephedra and its possible association with some biological and ecological features. Mol. Ecol. 25, 1192–1210. doi: 10.1111/mec.13538
Keywords: Juniperus, gymnosperms, conifers, polyploidy, genome size, flow cytometry
Citation: Farhat P, Hidalgo O, Robert T, Siljak-Yakovlev S, Leitch IJ, Adams RP and Bou Dagher-Kharrat M (2019) Polyploidy in the Conifer Genus Juniperus: An Unexpectedly High Rate. Front. Plant Sci. 10:676. doi: 10.3389/fpls.2019.00676
Received: 12 March 2019; Accepted: 06 May 2019;
Published: 22 May 2019.
Edited by:Michael R. McKain, University of Alabama, United States
Reviewed by:Dirk Carl Albach, University of Oldenburg, Germany
Petr Koutecký, University of South Bohemia, Czechia
Copyright © 2019 Farhat, Hidalgo, Robert, Siljak-Yakovlev, Leitch, Adams and Bou Dagher-Kharrat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Magda Bou Dagher-Kharrat, firstname.lastname@example.org