MBE Advance Access originally published online on July 23, 2008
Molecular Biology and Evolution 2008 25(10):2099-2108; doi:10.1093/molbev/msn157
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
Chalcone Synthase Gene Lineage Diversification Confirms Allopolyploid Evolutionary Relationships of European Rostrate Violets
* Nationaal Herbarium Nederland, Leiden University, Leiden, The Netherlands
Leerstoelgroep Biosystematiek, Wageningen University, Wageningen, The Netherlands
E-mail: vandenhof{at}nhn.leidenuniv.nl.
 |
Abstract |
|---|
 TOP Abstract IntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Phylogenetic relationships among and within the subsections of the genus Viola are still far from resolved. We present the first organismal phylogeny of predominantly western European species of subsection Rostratae based on the plastid trnS–trnG intron and intergenic spacer and the nuclear low-copy gene chalcone synthase (CHS) sequences. CHS is a key enzyme in the synthesis of flavonoids, which are important for flower pigmentation. Genes encoding for CHS are members of a multigene family. In Viola, 3 different CHS copies are present. CHS gene lineages obtained confirmed earlier hypotheses about reticulate relationships between species of Viola subsection Rostratae based on karyotype data. Comparison of the CHS gene lineage tree and the plastid species phylogeny of Viola reconstructed in this study indicates that the different CHS copies present in Viola are the products of both recent and more ancient duplications.
Key Words: chalcone synthase • gene lineage diversification • phylogeny • Viola subsection Rostratae • allopolyploidy • trnS–trnG
|
Introduction |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Speciation through hybridization is considered a common process in higher plants. Although hybrids between distantly related taxa are usually sterile, they can become fertile again by doubling their chromosome numbers. The resulting chimeric species can have 2 or more sets of chromosomes derived from different parental species; this is called allopolyploidy (Stebbins 1971
; Song et al. 1995; Bennett 2004; Hegarty and Hiscock 2005). In contrast with allopolyploidy, which may occur in connection with hybridization between taxa that are not very closely related, hybrid speciation without a change in chromosome number may occur in cases where the parental species are closely related and their primary hybrid is somewhat fertile; this process is called homoploid hybrid stabilization (Rieseberg 1997; Rieseberg et al. 2003; Abbott et al. 2005).
Polyploid evolution has been an important factor in the evolutionary history of land plants and continues to be so also in extant lineages such as the plant genus Viola (Violaceae). In fact, the first report of an infrageneric series of polyploid levels was from Viola (Miyaji 1913). The base chromosome number of Viola is believed to be x = 6 or x = 7, but the vast majority of north-temperate taxa have been shown to be paleo-allotetraploid with secondary base numbers of x = 10 or x = 12 (Nordal and Jonsell 1998; Marcussen and Nordal 1998; Karlsson et al. 2008). These are hereafter referred to as secondary diploids. Further, polyploidy based on these secondary diploid chromosome numbers has been demonstrated especially within the species-rich subsections of section Viola (Karlsson et al. 2008).
Within Section Viola subsection Rostratae Kupffer (sometimes treated as the separate section Trigonocarpea [Godr.] Vl. V. Nikitin), most species have retained the secondarily diploid chromosome number of 2n = 20. However, subsequent polyploidization events have led to the formation of higher ploids with chromosome numbers of 2n = 40 (octoploid), 60 (dodecaploid), or even 58 (subdodecaploid); these are hereafter referred to as secondary tetraploids and (sub)hexaploids, respectively. Nearly, all these secondary polyploids, a total of 10 species, are native to western Eurasia, and their relatively recent polyploid parentages have been investigated in a series of cytological studies in the late 1950s and early 1960s (fig. 1) (Valentine 1950, 1958; Moore and Harvey 1961; Harvey 1966). The subsection consists of about 50 species with a northern temperate distribution in North America and Eurasia. Most species have white to dark lilac flowers and grow in woodlands. Subsection Rostratae is characterized primarily by symplesiomorphic characters. Phylogenetic analyses based on nuclear ribosomal Internal Transcribed Spacer (nrITS) sequences have recovered that the subsection is paraphyletic with respect to a number of other north-temperate groups (Ballard et al. 1999; Yoo et al. 2005). In Europe, where subsection Rostratae is morphologically most diverse, the subsection has traditionally been subdivided in a variable number of morphologically defined groups, usually at the series level. Series Rosulantes is characterized by having a basal rosette and flowers produced only from the lateral aerial shoots; this growth form is found also in other sections and may be considered as primitive within the genus. Series Mirabiles differs from the Rosulantes in producing flowers also from the basal leaf rosette, series Arosulatae in lacking the basal rosette altogether, and series Repentes in being stoloniferous and producing flowers from the rosettes. However, the recognition of series is problematic for 2 main reasons. First, the series typically define small groups of species by a very limited number of autapomorphies, thereby rendering the remaining groups paraphyletic and defined by synplesiomorphies only. Second, several of these series cannot be considered monophyletic because of the alloploid relationships between taxa of different series.
|
Series Arosulatae defines a small group of 5 western Eurasian species. Species of the series are specialists of temporarily flooded habitats, rather than woodland, and are easily characterized by lacking leaf rosettes and by their leaf and stipule characters. Traditionally, 7 species have been recognized, but 2 (the East Asian Viola acuminata and the submediterranean Viola jordanii) must be omitted on the basis of having a leaf rosette and the quite different choices of habitat. The remaining 5 species are all Central European. Viola stagnina is the only secondary diploid with 2n = 20, whereas 3 species are secondary tetraploids with 2n = 40 (Viola canina, Viola elatior, and Viola pumila) next to a secondary subhexaploid with 2n = 58 (Viola lactea) (Moore and Harvey 1961
).
Especially, the study by Ballard et al. (1999) indicated that the taxonomy of the genus Viola needs revision and that more molecular phylogenetic studies are called for. Although the nrITS region used by Ballard et al. (1999) was useful for recognizing infrageneric groups of the genus Viola, nrITS is generally not useful for examining evolutionary relationships among polyploid lineages. This is because recombination and concerted evolution between orthologous nrITS copies often lead to retention of only one copy type and erasion of the other parental copy (Wendel et al. 1995; Álvarez and Wendel 2003). This is usually also the case in Viola as nearly all investigated species have retained only one nrITS copy type regardless of ploidal level (Ballard et al.1999; Malécot et al. 2007). The nrITS as a phylogenetic marker is therefore not suitable for recovering the reticulate relationships within the genus. Plastid markers generally also have the problem of retention of a single parental copy as these are usually uniparentally inherited in plants; furthermore, sequence variation in plastid markers is usually low (Corriveau and Coleman 1988; Taberlet et al. 2007). Again, reticulate relationships therefore remain obscure.
Álvarez and Wendel (2003) suggested using single- or low-copy nuclear markers to circumvent the problem of concerted evolution causing misleading phylogenetic reconstructions of polyploid species. Phylogenetic analysis of paralogous and orthologous copies of single- or low-copy genes in alloploid species is also a good method to reveal the parental contributors to alloploid genomes. This method has been successfully applied in numerous studies (e.g., Popp and Oxelman 2001; Smedmark et al. 2005).
We utilized the low-copy nuclear chalcone synthase (CHS) gene as a phylogenetic marker in Viola subsection Rostratae. As an independent data set we chose the trnS–trnG intergenic spacer and intron as plastid phylogenetic marker as this region proved to be sufficiently informative in Viola to assess interspecific relationships.
CHS is the first enzyme in the flavonoid synthesis pathway and is encoded by a small gene family (Durbin et al. 1995). Flavonoids are important secondary metabolites responsible for a multitude of tasks in plants, ranging from flower and fruit coloration and protection against UV radiation to pathogen defense and pollen development (Harborne 1994). In Viola cornuta, 3 different CHS gene copies were found to be expressed from early stages of flower coloration onward (Farzad et al. 2003). In general, genes of the CHS family consist of 1 intron flanked by 2 exons. There is high variation in the number of CHS copies among angiosperms. In asterids, the number of CHS copies ranges from a single copy in Antirrhinum (Sommer and Saedler 1986) to 6 copies in Ipomoea (Clegg and Durbin 2003) and 8 in Petunia (Koes et al. 1987). Similarly for the rosids, both Arabidopsis (Wang et al. 2007) and Populus (Tuskan et al. 2006) have 2 CHS copies, whereas both Vitis (Sparvoli et al. 1994; Jaillon et al. 2007) and V. cornuta cultivars (Farzad et al. 2003, 2005) have 3 CHS copies.
We collected different CHS paralogues in species of Viola subsection Rostratae and analyzed these phylogenetically to 1) test earlier hypotheses about reticulate relationships of several allopolyploid taxa based on karyotype data in subsection Rostratae (e.g., between V. stagnina, a possible Dutch endemic, and its closest relatives), 2) make a comparison with a species phylogeny of Viola subsection Rostratae based on sequences of the plastid trnS–trnG intron and intergenic spacer to infer how many duplications of CHS took place during the evolution of Viola.
|
Materials and Methods |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Taxon Sampling
In total, 30 Viola taxa with a predominantly western European origin were sampled, of which 21 taxa belong to Viola subsection Rostratae. The 9 taxa outside subsection Rostratae represent sections Andinium, Boreali-Americanae, Chamaemelanium, Erpetion, and Melanium and subsection Viola of section Viola. These species appeared to be either closely or more distantly related to the species of subsection Rostratae in a previous molecular phylogenetic study of Viola (Ballard et al. 1999
). DNA was obtained from freshly collected material from the field and from herbarium collections. For reconstruction of the CHS gene lineage tree, 2 different parts of the gene were sampled, the intron and exon 2. Exon 2 lineages available in GenBank from representatives of major Angiosperm clades were included in the analysis to find out whether the different CHS copies present in Viola are the products of recent or more ancient duplications. The following lineages were sampled: Gymnosperms: GbCHS (Ginkgo biloba, AY647263 [GenBank] ) and PsCHS (Pinus sylvestris, X60754 [GenBank] ); Monocots: IhCHS (Iris x hollandica, AB232914 [GenBank] ), HvCHS (Hordeum vulgare, X58339 [GenBank] ), and ZmCHS (Zea mays, AY728478 [GenBank] , X60204 [GenBank] ); Core eudicots: VvCHS (Vitis vinifera, AB015872 [GenBank] , AB066275 [GenBank] , EF192464 [GenBank] , AM 454341, X75969 [GenBank] ); Rosids: GmCHS (Glycine max, AY262686 [GenBank] ), PsCHS (Pisum sativum, D88263 [GenBank] , D88262 [GenBank] , D88261 [GenBank] , D88260 [GenBank] , X63333 [GenBank] ), PtCHS (Populus spp., DQ371804 [GenBank] , EF147137 [GenBank] , EF147091 [GenBank] , DQ371802 [GenBank] ), and VcCHS (V. cornuta cultivar, AY497407 [GenBank] , AY497414 [GenBank] ); Asterids: AmCHS (Antirrhinum majus, X03710 [GenBank] ), DcCHS (Daucus carota, D16255 [GenBank] ), and PhCHS (Petunia hybrida, X14597 [GenBank] ).
The phylogenetic analyses performed were all rooted differently. There were several reasons for this. First of all, plastid sequences of non-Violaceae were not used for phylogenetic analyses. The Angiosperm Phylogeny Group topology was used to constrain the analyses instead (see below). For the plastid phylogeny, closely related genera of Viola were used as outgroups. Second, CHS intron sequences outside Viola could not be aligned with CHS intron Viola data because of too high sequence divergence. For the CHS intron analyses, we therefore tentatively used Viola CHS3 as outgroup. Third, CHS gene duplication events could only be assessed with a broad taxonomic sampling. For the CHS exon 2 analyses, we therefore used gymnosperm lineages for rooting.
DNA Extraction, Polymerase Chain Reaction Amplification, Cloning, and Sequencing
Total genomic DNA was extracted using the Dneasy Plant Mini Kit (Qiagen, Hilden, Germany) and the cetyltrimethylammonium bromide (CTAB) method of Doyle JJ and Doyle JL (1987) with some modifications. Leaf material was ground using a Ratch Mill. In total, 750 µl CTAB buffer (Doyle JJ and Doyle JL 1987) was added to the ground material together with proteinase K and RNase. After incubating for 30 min at 60 °C, 750 µl of chloroform–isoamyl alcohol (24:1) was added. The samples were briefly vortexed and then centrifuged for 10 min at 12,000 rpm. The upper aqueous layer was transported to a clean 2-ml tube. A total of 500 µl chloroform–isoamyl alcohol (24:1) was added, and the samples were again centrifuged at 12,000 rpm. After 5 min spinning, the upper phase was transferred to a new 2-ml tube. The DNA was then precipitated by adding cold 500 µl isopropanol. The samples were shaken 5–10 min and subsequently centrifuged for 15 min at 12,000 rpm. The supernatant was removed, and 70% ethanol was added. The samples were subsequently shaken vigorously for 2 min, after which the ethanol was poured off. The remaining ethanol was removed by evaporation. The resulting DNA pellet was dissolved in 200 µl 0.1x Tris-EDTA buffer.
In total, one plastid region (trnS–trnG spacer and intron) and one nuclear region (CHS intron and exon 2) were amplified and sequenced. Polymerase chain reaction (PCR) amplification of the plastid spacer and intron was performed with primers designed by Shaw et al. (2005). For the trnS spacer, the primers trnSGCU (5'-AGA TAG GGA TTC GAA CCC TCG GT-3') and trnG2S (5'-TTT TAC CAC TAA ACT ATA CCC GC-3') were used. The trnG intron was amplified with the primers trnGUUC (5'-GTA GCG GGA ATC GAA CCC GCA TC-3') and trnG2G (GCG GGT ATA GTT TAG TGG TAA AA).
The primers CHSX1F (5'-AGG AAA AAT TCA AGC GCA TG-3') and CHSX2RN (5'-TTC AGT CAA GTG CAT GTA ACG-3') designed by Strand et al. (1997) were used for amplifying the CHS intron. The primers CHS forward (5'-TAY CAR CAR GGN TGY TTY GC-3') and CHS reverse (5'-GGR TGD GCD ATC CAR AAV A-3') from Farzad et al. (2003) were used to amplify exon 2 of the CHS gene (fig. 2a). The generated intron and exon sequences did not overlap as the "intermediate" part turned out to be too large and too heterogeneous for this. PCR fragments for several Viola species of the CHS intron are shown in figure 2b. Per individual, 12 clones were analyzed, and consensus sequences were compiled from 3 to 7 individual clones.
|
PCR amplification conditions for the trnS spacer consisted of denaturing for 50 s at 95 °C, annealing for 1 min at 53 °C, and extension for 2 min at 72 °C. This cycle was repeated 35 times. The trnG intron was amplified with the same conditions except for the annealing temperature which was 56 °C. PCR amplification conditions for the CHS intron consisted of denaturing for 1 min at 95 °C, annealing for 90 s at 53 °C, and extension for 2 min at 72 °C. This cycle was repeated 35 times. The conditions for amplifying the CHS exon 2 consisted of denaturing for 45 s at 95 °C, annealing for 1 min at 55 °C, and extension for 1 min at 72 °C. This cycle was repeated 40 times.
PCR products were purified using the Promega Wizard Purification System; cloned using the pGEM-T Easy Vector System; and sequenced using the M13 primers with 30 s denaturing at 95 °C, annealing for 30 s at 50 °C, and extension for 1 min at 72 °C. This cycle was repeated 35 times. PCR products were purified and analyzed on an ABI 377 (Applied Biosystems Inc., Foster City, CA) or a MegaBACE Sequence Analyzer 4.0 (Amersham Biosciences, Uppsala, Sweden) automated sequencer using the manufacturers' protocols.
Phylogenetic Analyses
DNA sequences were aligned using McClade 4.06 (Maddison DR and Maddison WP 2003) with the pairwise alignment option and manual adjustment where necessary. Individual insertion and deletion events were manually added as additional binary characters.
MrModeltest version 2.2 (Nylander 2004) was used to find the best model of sequence evolution (Posada and Crandall 1998). The models used for Bayesian analyses were the symmetrical model with separate gamma distributions and a separate proportion of invariant sites for CHS exon 2 (SYMIG model), the General Time Reversible model with gamma distribution for CHS intron (GTRG model), and the General Time Reversible model with gamma distribution and a separate proportion of invariant sites for trnS-trnG (GTRIG model). Maximum parsimony (MP) analyses were carried out with PAUP* 4.0b10 (Swofford 2001). Phylogenies were obtained using the heuristic search option, with 20 random sequence additions and Tree Bisection-Reconnection branch swapping. After each sequence addition, a maximum of 10,000 trees were saved.
For MP, bootstrap support (Felsenstein 1985) was calculated with 2,000 bootstrap replicates, using only 10 random sequence additions each bootstrap replicate. After every random sequence addition replicate, a maximum of 2,500 trees were saved. Bayesian inference analyses were performed using MrBayes 3.1 (Huelsenbeck and Ronquist 2001). Markov chain Monte Carlo (MCMC) analyses were run for 8 million generations with 5 simultaneous MCMCs, saving one tree per 100 generations. The burn-in values were identified using the program Tracer 1.3 (Rambaut and Drummond 2004).
To convert the CHS gene tree composed of multiple paralogous lineages from allopolyploid taxa into a species tree to assess gene duplications, GeneTree version 1.3 (Page and Charleston 1997) was used. The analyses were run with default settings. GeneTree requires fully resolved organismal and gene trees as input. For the organismal tree, 1 of the 200,000 fully resolved most parsimonious trees (MPTs) of trnS–trnG data was chosen randomly. This analysis was constrained for all nonviolets to the latest angiosperm phylogeny topology as depicted on the Angiosperm Phylogeny Web site (version 8, June 2007) (http://www.mobot.org/MOBOT/research/Apweb/). For the CHS exon 2 gene tree, the 95% most probable Bayesian tree was used.
Homology Assessment of CHS Copies
The CHS lineages found were assigned to different copies based on size, sequence divergence, and phylogenetic position (Helariutta et al. 1996; Doyle and Davis 1998; Smedmark et al. 2005). CHS fragments within one species with only minor divergence and gaps were interpreted as alleles. The different CHS copies of V. cornuta published by Farzad et al. (2003) always ended up in a single clade in all analyses performed here. We therefore used a single representative sequence only. When size difference and sequence divergence were more apparent, for example, by the presence of large indel events, the CHS fragment was treated as a paralogous copy. Our classification of alleles and paralogous copies was further confirmed by topological positions in the phylogenies obtained.
|
Results |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
trnS–trnG
MP analyses of the trnS–trnG alignment produced a total of 200,000 MPTs with 537 steps (consistency index [CI] = 0.8239; retention index [RI] = 0.7312). The majority rule consensus tree (data not shown) has a similar topology as the Bayesian tree (fig. 3). We plotted both the bootstrap support (BS) values and posterior probability index (PPI) values on the latter. All species sampled of Viola subsection Rostratae ended up in 5 different, poorly to well-supported subclades (<50–98% BS; 0.56–1.00 PPI). The largest subclade consists of V. stagnina, V. elatior, V. lactea, V. canina, Viola sieheana, V. jordanii, Viola oligyrtia, Viola rupestris, and V. pumila.
|
CHS Intron and Exon 2
Three copies of the CHS intron were found in the sampled species of Viola subsection Rostratae (figs. 2b and 4
). Copies 1 (CHS1, 600 base pairs [bp]) and 2 (CHS2, 735–1,100 bp) have a relatively similar sequence identity and were found in all Viola species sampled with the exception of Viola biflora, Viola banksii, and Viola pubescens. The third copy was not found in all species and was the largest sized (CHS3, 775–1,160 bp). The CHS3 copy is probably present in more taxa, but amplification failures probably led to an undersampling of this particular copy. It differed quite substantially from the other 2 copies in size and sequence similarity. In contrast with CHS1, multiple paralogs/orthologs were found in CHS2 and CHS3.
|
The complete CHS intron alignment consisted of 2,980 bp after exclusion of a 64-bp segment that was too variable for proper alignment. A total of 302 characters were phylogenetically informative, of which 12 were indel characters (indels varying in size between 5 and 422 bp). Two indels, found in CHS1 and CHS2, seemed to be the result of slip-strand mispairing as many repeats were found in these regions. The first (TGATTT) and second repeats (TGTT) were repeated up to 4 times. The other 10 indels lacked a repetitive structure. Most of the indels occurred in CHS2. In CHS3, 2 large indels were found of 182 and 422 bp, respectively.
MP analyzed of CHS intron sequences produced 200,000 MPTs (921 steps, CI = 0.7828; RI = 0.8824). The majority rule consensus tree (fig. 4) had a similar topology as the Bayesian tree (data not shown).
At least 2 different copies of CHS exon 2 were found in Viola. Both copies had a similar sizes (984 bp) and were retrieved from almost all Viola species analyzed. They differed substantially in sequence similarity. Of the 984 bp retrieved, 498 were phylogenetically informative. One autapomorphic gap was found. Bayesian analyses of CHS exon 2 sequences produced a topology similar to MP (data not shown), in which 2 main clades were present comprising the different copies of CHS exon 2.
Reconciliation of Gene Tree and Species Tree
A reconciled tree (fig. 5) reconstructed with the program GeneTree (Page and Charleston 1997) was used to visualize CHS exon 2 duplications during the evolution of Viola. It seems that assuming 6 CHS gene duplication events is sufficient to make the gene and species tree congruent.
|
|
Discussion |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Polyploidy in Viola Subsection Rostratae
The internal topology of the CHS2 intron clades (fig. 4) is in general agreement with previously inferred relationships between parental species and their allopolyploid hybrids in Viola subsection Rostratae (Moore and Harvey 1961
). Moore and Harvey (1961) could recognize the parental karyotypes in the genomes of artificially constructed allopolyploid Viola hybrids by the unique size and shape of the chromosomes. Subsequently, they used observations on chromosome pairing to formulate hypotheses regarding the origin of allopolyploids (fig. 1). In their study, 5 different types of genomes (A–E) could be recognized, each referring to the secondary diploid level (2n = 20). Only the A and C genomes occurred in extant secondary diploids, namely in Viola reichenbachiana (A) and V. stagnina (C). The other 4 genome types were found only in combination with other genome types in the secondary tetraploid (2n = 40) or subhexaploid taxa (2n = 58). From this, they concluded that V. stagnina (C) contributed a C genome to V. canina (BC) and its close relative V. lactea (BCE) and possibly also to V. pumila (CD). Similarly, V. reichenbachiana (A) would have contributed an A genome to Viola riviniana f. riviniana (AB). Thus, the species possessing the B genome was involved in the origin of both V. riviniana f. riviniana (AB) and V. canina (BC). The authors attributed the 3 "missing" genomes (B, D, and E) to secondary diploids species that might have become subsequently extinct, at least in Europe. Viola elatior was not included in this study. In the CHS2 intron tree, V. canina was found to have one orthologue in common with V. stagnina (corresponding to genome C) and a second in common with V. riviniana (corresponding to genome B). The second orthologue of V. riviniana was found to be closely related to V. reichenbachiana (corresponding to genome A). Like V. canina, V. pumila had one orthologue in common with V. stagnina (genome C), whereas its second copy was found to be closely related to V. canina and V. riviniana suggesting that genome D could have been derived from genome B.
The phylogenetic position of the gene lineages retrieved from V. elatior suggests that this particular species probably contains the C genome because its CHS2 intron copy ended up close to V. stagnina. The chromosome number of V. elatior (2n = 40) suggests that it is a secondary tetraploid, but we were not successful in detecting more than one orthologue in this species. Unpublished isozyme studies also reveal a lower number of allozymic bands than usual in species of this ploidal level. These findings, together with earlier observations of quadrivalents in the meiosis of the species (Clausen 1927), indicate that V. elatior may be an autopolyploid derivative of some stagnina-like ancestral species possessing the C genome.
In contrast with the Arosulatae series, which was found to be monophyletic for CHS1 and part of CHS2, CHS lineages retrieved from species assigned to the Rosulantes and Mirabiles series did hardly ever end up in the same clades. This is probably caused by the fact that the morphological characters used to delimit these series are phylogenetically uninformative. Unknown hybridization and polyploidization events within and between these series probably also cause paraphyly.
The small series Mirabiles consists of one secondarily diploid (2n = 20) species, Viola mirabilis, and 2 secondary tetraploids (2n = 40). The local endemic Viola pseudo–mirabilis of Les Grands Causses in southern France has been variously interpreted on morphological grounds as an intermediate between V. mirabilis and V. riviniana (Valentine et al. 1968) or as a polyploid derivative of V. mirabilis and V. reichenbachiana (Espeut 1999). The latter view has later been confirmed by unpublished isozyme data and a chromosome count of 2n = 40 (Verlaque and Espeut 2007). In the present study, V. pseudo–mirabilis ends up as sister to V. riviniana (61% BS; 0.99 PPI), which is not in contradiction to the previous findings because V. riviniana is itself a polyploid derivative of V. reichenbachiana. The other secondary tetraploid Mirabiles species, Viola willkommii, endemic to northern Spain, is morphologically rather similar to V. mirabilis but differs considerably in choice of habitat. This particular species is a secondary tetraploid and has been supposed to have originated from V. mirabilis and V. rupestris (Marcussen T, personal communication). Our CHS data confirm the parentage of V. mirabilis, but the second parent V. willkommii remains unclear. Two paralogues of the CHS2 intron were retrieved from V. willkommii, of which one ended up in a strongly supported clade (96% BP; 0.97 PPI) with V. mirabilis and the other in an unresolved polytomy. The CHS1 intron fragment retrieved from V. willkommii ended up in a basal dichotomy with the lineage of V. rupestris. Future genome type data should be collected of the Mirabiles series to confirm hypotheses about reticulate relationships suggested by the CHS gene lineages obtained here.
Closest Relatives of V. stagnina
Viola subsection Rostratae as mentioned before has been taxonomically subdivided in series Arosulatae, Mirabiles, Repentes, and Rosulantes (Valentine 1958). Viola stagnina is considered to be a member of the Arosulatae series together with V. canina, V. elatior, V. lactea, and V. pumila, based on the lack of a basal leaf rosette. This taxonomic placement is supported by the fact that in our study, 2 Arosulatae clades are present in the CHS intron gene lineage tree (fig. 4), one consisting of CHS1 intron lineages and one of CHS2 intron lineages. These lineages were all retrieved from species assigned to the Arosulatae series.
Based on previous morphological and karyological studies and our results, we can conclude that the closest relatives of V. stagnina are V. pumila, V. elatior, V. canina, and V. lactea. Fingerprinting techniques are currently being applied to assess gene flow between different European populations of V. stagnina and its close relatives to determine whether a Dutch variety of V. stagnina deserves a different taxonomic status.
CHS Lineage Diversification in Viola
In the reconciled tree (fig. 5), 2 main clades are present. Unfortunately, GeneTree does not provide statistical support for individual nodes. The congruence in topology with Huang et al. (2004) and Yamazaki et al. (2001) indicates that a general phylogenetic signal was recovered, though. The first clade consists of monocot, core eudicots, and Viola representatives. This clade indicates that at least one duplication event in the CHS gene family took place before the split between the monocots and the eudicots. The second main clade consists of rosid, asterid, and Viola representatives. This indicates that another duplication event in the CHS gene family took place during the split between the core eudicots and rosids/asterids. Similar results were also found by Huang et al. (2004).
Yang and Gu (2006) also describe multiple rounds of CHS gene duplications during the evolution of the angiosperms. According to these authors, the most ancestral gene lineages originated during the divergence of different plant families, such as Solanaceae, Convolvulaceae, and Asteraceae in a first round of duplications. Derived CHS genes further duplicated and diverged, which led to the occurrence of various CHS plant family–specific genes in subsequent rounds of duplication.
The CHS intron gene lineage tree suggests that 3 CHS copies are present in Viola, whereas the CHS exon 2 gene lineage tree (data not shown) only contains 2 copies. The oldest duplication event in clade I (fig. 5) suggests the possibility of the presence of a third copy of CHS in Viola. This might also explain why CHS intron copies 1 and 2 have a close resemblance. The fact that not all copies were retrieved does not mean they are not there, but could explain why taxa from different sections end up nested in section Viola in the reconciled tree. Our interpretation of the CHS paralogues in Viola is different from Farzad et al. (2003, 2005). We consider the CHS paralogues in V. cornuta to be different alleles, whereas the latter study identified them as different copies. The plant analyzed by Farzad et al. (2003, 2005) was a garden cultivar of hybrid origin. The nominal species is in itself a high polyploid (Marcussen et al. 2008). Furthermore, our analyses of a larger sample of different Viola lineages showed that the interpretation of Farzad et al. (2005) was incomplete. Farzad et al. (2005) showed that the 3 CHS paralogous in V. cornuta are all still expressed and fully functional. Expression patterns were found to be slightly different, which might indicate subfunctionalization. Subfunctionalization of duplicate CHS genes in angiosperms appears to have happened by differentiation of their regulatory elements (Yang and Gu 2006). It would be interesting to further investigate the mechanisms of subfunctionalization in a wider array of Viola species.
|
Appendix |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Voucher information of taxa sampled. Behind the accession numbers of the CHS intron sequences, the copy numbers are indicated, corresponding to those in figure 4.
|
Acknowledgements |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
We would like to acknowledge T. Marcussen for providing material and for his valuable feedback on earlier versions of this manuscript. Furthermore, we like to acknowledge P. van Beers, R. L. Eckstein, H. den Held, F. Hellberg, W. J. Holverda, the late R. van der Meijden, R. van Moorsel, I. Nordal, P. B. Pelser, N. Schidlo, E. J. Weeda, and B. Wijlens for providing material and assistance in the field; Kew Botanical Gardens for providing DNA samples; and M. C. M. Eurlings, B. J. van Heuven, L. McIvor, S. Bollendorff, and D. L. V. Co for technical support. We would like to acknowledge the following organizations and departments for providing access to plant material and collecting permits: Staatsbosbeheer, Natuurmonumenten, Overijssels Landschap, and Kalmar county environmental department.
|
Footnotes |
|---|
Neelima Sinha, Associate Editor
|
References |
|---|
|
TOPAbstractIntroductionMaterials and MethodsResultsDiscussionAppendixAcknowledgementsReferences |
|---|
Abbott RJ, Ireland HE, Joseph L, Davies MS, Rogers HJ. Recent plant speciation in Britain and Ireland: origins, establishment and evolution of four new hybrid species. Biol Environ Proc R Ir Acad (2005) 105:173–183.[CrossRef]
Álvarez I, Wendel JF. Ribosomal ITS sequences and plant phylogenetic inference. Mol Phylogenet Evol (2003) 29:417–437.[CrossRef][Web of Science][Medline]
Ballard HE, Sytsma KJ, Kowal RR. Shrinking the violets: phylogenetic relationships of infrageneric groups in Viola (Violaceae) based on internal transcribed spacer DNA sequences. Syst Bot (1999) 23:439–458.[CrossRef]
Bennett MD. Perspectives on polyploidy in plants—ancient and neo. Biol J Linn Soc (2004) 82:411–423.[CrossRef][Web of Science]
Clausen J. Chromosome number and the relationship of species in the genus Viola. Ann Bot (1927) 61:677–714.
Clegg MT, Durbin ML. Tracing floral adaptations from ecology to molecules. Nat Rev Genet (2003) 4:206–215.[CrossRef][Web of Science][Medline]
Corriveau JL, Coleman AW. Rapid screening method to detect potential biparental inheritance of plastid DNA and results from over 200 angiosperm species. Am J Bot (1988) 75:1443–1458.[CrossRef][Web of Science]
Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull Bot Soc Am (1987) 19:11–15.
Doyle JJ, Davis JI. Homology in molecular phylogenetics: a parsimony perspective. In: Molecular systematics of plants II—Soltis DE, Soltis PS, Doyle JJ, eds. (1998) Norwell (MA): Kluwer Academic Publishers. 101–131.
Durbin ML, Learn GH, Huttley GA, Clegg MT. Evolution of the chalcone synthase gene family in the genus Ipomoea. Proc Natl Acad Sci USA (1995) 92:3338–3342.
Espeut M. Errata et addenda de l'article: "Approche du genre Viola dans le Midi Mediterranéen français". Monde Pl (1999) 467:7–9.
Farzad M, Griesbach R, Hammond J, Weiss MR, Elmendorf HG. Differential expression of three key anthocyanin biosynthetic genes in a color-changing flower, Viola cornuta cv. Yesterday, today and tomorrow. Plant Sci (2003) 165:1333–1342.
Farzad M, Soria-Hernanz DF, Altura M, Hamilton MB, Weiss MR, Elmendorf HG. Molecular evolution of the chalcone synthase gene family and identification of the expressed copy in flower petal tissue of Viola cornuta. Plant Sci (2005) 168:1127–1134.
Felsenstein J. Confidence-limits on phylogenies—an approach using the bootstrap. Evolution (1985) 39:783–791.[CrossRef][Web of Science]
Harborne JB. The flavonoids: advances in research since 1986 (1994) London: Chapman and Hall.
Harvey MJ. Cytotaxonomic relationships between the European and North American rostrate violets. New Phytol (1966) 65:469–476.[CrossRef][Web of Science]
Hegarty MJ, Hiscock SJ. Hybrid speciation in plants: new insights from molecular studies. New Phytol (2005) 165:411–423.[CrossRef][Web of Science][Medline]
Helariutta Y, Kotilainen M, Elomaa P, Kalkkinen N, Bremer K, Teeri TH, Albert VA. Duplication and functional divergence in the chalcone synthase gene family of Asteraceae: evolution with substrate change and catalytic simplification. Proc Natl Acad Sci USA (1996) 93:9033–9038.
Huang JX, Qu LJ, Yang J, Gu HY. A preliminary study on the origin and evolution of chalcone synthase (CHS) gene in Angiosperms. Acta Bot Sin (2004) 46:10–19.
Huelsenbeck JP, Ronquist F. MrBayes: Bayesian inference of phylogeny. Bioinformatics (2001) 17:754–755.
Jaillon O, Aury J-M, Noel B, et al, (56 co-authors). The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature (2007) 449:463–468.[CrossRef][Web of Science][Medline]
Karlsson T, Marcussen T, Wind P, Jonsell B, Forthcoming. Violaceae—Karlsson T, Nordica Flora, eds. (2008) Vol. 6. Stockholm (Sweden): The Bergius Foundation.
Koes RE, Spelt CE, Mol J, Gerats A. The chalcone synthase multigene family of Petunia hybrida (V30): sequence homology, chromosomal localization and evolutionary aspects. Plant Mol Biol (1987) 10:159–169.[Medline]
Maddison DR, Maddison WP. Macclade 4.06 (2003) Sunderland (MA): Sinauer Associates.
Malécot V, Marcussen T, Munzinger J, Yockteng R, Henry M. On the origin of the sweet-smelling Parma Violet cultivars (Violaceae): wide intraspecific hybridization, sterility, and sexual reproduction. Am J Bot (2007) 94:29–41.
Marcussen T, Nordal I. Viola suavis, a new species in the Nordic flora, with analyses of the relation to other species in the subsection Viola (Violaceae). Nordic Journal of Botany (1998) 18:221–237.[CrossRef][Web of Science]
Miyaji Y. Untersuchungen über die chromosomezahlen bei einigen Viola-Arten. Bot Mag Tokyo (1913) 27:443–460.
Moore DM, Harvey MJ. Cytogenetic relationships of Viola lactea SM. and other West European Arosulate violets. New Phytol (1961) 60:85–95.[Medline]
Nordal I, Jonsell B. A phylogeographic analysis of Viola rupestris: three post-glacial migration routes into the Nordic area? Bot J Linn Soc (1998) 128:105–122.[Web of Science]
Nylander JAA. MrModeltest v2 (2004) Uppsala (Sweden): Evolutionary Biology Centre: Uppsala University.
Page RDM, Charleston MA. From gene to organismal phylogeny: reconciled trees and the gene tree/species tree problem. Mol Phylogenet Evol (1997) 7:231–240.[CrossRef][Web of Science][Medline]
Popp M, Oxelman B. Inferring the history of the polyploid Silene aegaea (Caryophyllaceae) using plastid and homoeologous nuclear DNA sequences. Mol Phylogenet Evol (2001) 20:474–481.[CrossRef][Web of Science][Medline]
Posada D, Crandall KA. Modeltest: testing the model of DNA substitution. Bioinformatics (1998) 14:817–818.
Rambaut A, Drummond AJ. Tracer. Version 1.3 (2004) Oxford: University of Oxford.
Rieseberg LH. Hybrid origins of plant species. Annu Rev Ecol Syst (1997) 28:359–389.[CrossRef][Web of Science]
Rieseberg LH, Raymond O, Rosenthal DM, Lai Z, Livingstone K, Nakazato T, Durphy JL, Schwarzback AE, Donovan LA, Lexer C. Major ecological transitions in wild sunflowers facilitated by hybridization. Science (2003) 301:1211–1216.
Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE, Small RL. The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am J Bot (2005) 92:142–166.
Smedmark JEE, Eriksson T, Bremer B. Allopolyploid evolution in Geinae (Colurieae: Rosaceae)—building reticulate species trees from bifurcating gene trees. Org Divers Evol (2005) 5:275–283.[CrossRef]
Sommer H, Saedler H. Structure of the chalcone synthase gene of Antirrhinum majus. Mol Gen Genet (1986) 202:429–434.[CrossRef][Web of Science]
Song K, Lu P, Tank K, Osborn TC. Rapid genome change in synthetic polyploids of Brassica and its implications for polyploid evolution. Proc Natl Acad Sci USA (1995) 92:7719–7723.
Sparvoli F, Martin C, Scienza A, Gavazzi G, Tonelli C. Cloning and molecular analysis of structural genes involved in flavonoid and stilbene biosynthesis in grape (Vitis vinifera L.). Plant Mol Biol (1994) 24:743–755.[CrossRef][Web of Science][Medline]
Stebbins GL. Chromosomal evolution in higher plants (1971) London: Edward Arnold.
Strand AE, Leebens-mack J, Milligan BG. Nuclear DNA-based markers for plant evolutionary biology. Mol Ecol (1997) 6:113–118.[CrossRef][Medline]
Swofford DL. PAUP*: phylogenetic analysis using parsimony (* and other methods). Version 4 (2001) Sunderland (MA): Sinauer Associates.
Taberlet P, Coissac E, Pompanon F, Gielly L, Miquel C, Valentini A, Vermat T, Corthier G, Brochmann C, Willerslev E. Power and limitations of chloroplast trnL (UUA) intron for plant DNA barcoding. Nucleic Acids Res (2007) 35:e1–e8.
Tuskan GA, DiFazio S, Jansson S, et al, (104 co-authors). The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science (2006) 313:1596–1604.
Valentine DH. The experimental taxonomy of two species of Viola. New Phytol (1950) 49:193–212.
Valentine DH. Cytotaxonomy of the rostrate violets. Proc Linn Soc Lond (1958) 169:132–134.
Valentine DH, Merxmuller H, Schmidt A. Viola L. In: Flora Europaea 2—Tutin TG, Heywood VH, Burges NA, Moore DM, Valentine DH, Walters SM, Webb DA, eds. (1968) Cambridge: Cambridge Press. 270–282.
Verlaque R, Espeut M. IAPT/IOPB chromosome data 3. Taxon (2007) 56:209.
Wang WK, Schaal BA, Chiou YM, Murakami N, Ge XJ, Huang CC, Chiang TY. Diverse selective modes among orthologs/paralogs of the chalcone synthase (Chs) gene family of Arabidopsis thaliana and its relative A. halleri ssp. gemmifera. Mol Phylogenet Evol (2007) 44:503–520.[CrossRef][Web of Science][Medline]
Wendel JF, Schnabel A, Seelanan T. Bi-directional interlocus concerted evolution following allopolyploid speciation in cotton (Gossypium). Proc Natl Acad Sci USA (1995) 92:280–284.
Yamazaki Y, Suh D-Y, Sitthithaworn W, Ishiguro K, Kobayashi Y, Shibuya M, Ebizuka Y, Sankawa U. Diverse chalcone synthase superfamily enzymes from the most primitive vascular plant, Psilotum nudum. Planta (2001) 214:75–84.[Web of Science][Medline]
Yang J, Gu H. Duplication and divergent evolution of the CHS and CHS like genes in the chalcone synthase (CHS) superfamily. Chin Sci Bull (2006) 51:505–509.[CrossRef]
Yoo KO, Jang SK, Lee WT. Phylogeny of Korean Viola based on ITS sequences. Korean J Plant Taxon (2005) 35:7–23.
CiteULike 
Connotea 
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||