Are both sympatric species Ilex perado and Ilex canariensis secretly hybridizing? Indication from nuclear markers collected in Tenerife
© Manen; licensee BioMed Central Ltd. 2004
Received: 17 August 2004
Accepted: 18 November 2004
Published: 18 November 2004
Intra-specific and intra-individual polymorphism is frequently observed in nuclear markers of Ilex (Aquifoliaceae) and discrepancy between plastid and nuclear phylogenies is the rule in this genus. These observations suggest that inter-specific plastid or/and nuclear introgression played an important role in the process of evolution of Ilex. With the aim of a precise understanding of the evolution of this genus, two distantly related sympatric species collected in Tenerife (Canary Islands), I. perado and I. canariensis, were studied in detail. Introgression between these two species was previously never reported. One plastid marker (the atpB-rbcL spacer) and two nuclear markers, the ribosomal internal transcribed spacer (ITS) and the nuclear encoded plastid glutamine synthetase (nepGS) were analyzed for 13 and 27 individuals of I. perado and I. canariensis, respectively.
The plastid marker is intra-specifically constant and correlated with species identity. On the other hand, whereas the nuclear markers are conserved in I. perado, they are highly polymorphic in I. canariensis. The presence of pseudogenes and recombination in ITS sequences of I. canariensis explain this polymorphism. Ancestral sequence polymorphism with incomplete lineage sorting, or past or recent hybridization with an unknown species could explain this polymorphism, not resolved by concerted evolution. However, as already reported for many other plants, past or recent introgression of an alien genotype seem the most probable explanation for such a tremendous polymorphism.
Data do not allow the determination with certitude of the putative species introgressing I. canariensis, but I. perado is suspected. The introgression would be unilateral, with I. perado as the male donor, and the paternal sequences would be rapidly converted in highly divergent and consequently unidentifiable pseudogenes. At least, this study allows the establishment of precautionary measures when nuclear markers are used in phylogenetic studies of genera having experienced introgression such as the genus Ilex.
Aquifoliaceae comprise one genus, Ilex  and approximately 400 species. The fossil record indicates that the genus was cosmopolitan during the Eocene. It is now largely extinct in Australia, Europe and Africa where only few species persist. Most diversity is currently found in South-America and in Southeast-Asia. They are evergreen or deciduous trees or bushes living in warm-moist-temperate, sub-tropical, tropical or montane-tropical areas.
The molecular phylogeny of the genus Ilex [2, 3] shows that systematic relationships are still not well understood. The plastid phylogeny (inferred from the atpB-rbcL spacer, rbcL and trnL-F) is highly correlated with the geographic distribution of extant species. Four chloroplast clades are found: one exclusively Eurasian clade, one exclusively American clade and two different North-American/Asian clades (one of them comprising most of the deciduous species among other evergreen species). On the other hand, the nuclear phylogeny (inferred from ribosomal ITS and the 5S RNA spacer) is incongruent with the plastid phylogeny, suggesting frequent interlineage hybridizations. The nuclear phylogeny is not correlated with the geographic distribution of extant species.
Any of the plastid or the nuclear phylogeny corroborates previous morphological or biosystematic studies . Using chloroplast RFLPs, trnL-trnF sequencing and nuclear ITS sequencing, a study of Asian Ilex of the Bonin Island and of the Ryukyu Island  confirmed that hybridization played a role in this region, leading to interspecific introgressions independently observed on both Islands. RAPD data indicate that the Japanese species Ilex leucoclada M. is highly polymorphic . During its history, the genus Ilex probably experienced frequent incomplete lineage sorting and nuclear and/or cytoplasmic introgression, making the study of its history very complex.
Few data are reported on the chromosome number of Ilex . The basic haploid number is 20, with deviation to 17, 18 and 19. From the 27 chromosome numbers available for the genus Ilex, three species are tetraploid (I. anomala, I. verticillata and I. argentina) and one species is hexaploid (I. pedunculosa), indicating probable hybridizations between species having divergent genomic background (alloploidy).
A previous study  showed that individuals of many species of Ilex contain polymorphic nuclear sequences (ITS and 5S rDNA spacer). Except for I. purpurea and I. guianensis, only one individual was studied per species. The sampling being too low for a correct evaluation of this intraspecific polymorphism, an exhaustive study of one plastid marker (the atpB-rbcL spacer) and two nuclear markers (the ribosomal internal transcribed spacers, ITS, and the nuclear encoded plastid glutamine synthetase, nepGS) was undertaken on several individuals of I. perado and I. canariensis collected in Tenerife (Canary Islands). These species were chosen because, based on DNA data, they are not closely related [2, 3] and are growing sympatrically in Canary Islands. Both species are morphologically variable but few characters allow species identification. I. canariensis is endemic of Canary Islands, whereas I. perado has a wider distribution in Spain, Portugal, North-Africa and Canary Islands. Natural or artificial hybridization between both species was never reported. The data show that, contrarily to I. perado, I. canariensis has highly polymorphic ITS and nepGS sequences. The aim of this study was to (1) explain the polymorphism observed in ITS of I. canariensis by an investigation of its pattern of substitution and its functionality, (2) determine the evolutionary mechanisms responsible of this polymorphism and (3) focus on ITS evolution and consequences for phylogenetic reconstruction of the genus Ilex.
As the ITS sequences found in I. canariensis are extremely polymorphic, it was interesting to observe their relationships with available ITS sequences previously investigated species by Manen et al. . Figure 1 shows one of the most parsimonious ITS tree of I. perado and I. canariensis sequences found in Tenerife altogether with 43 ITS sequences of other Ilex species. Thick bars indicate internal branches conserved in the consensus tree. The closest possible outgroups for the genus Ilex, are Helwingia and Phyllonoma. However these genera are so isolated systematically that their use to root Ilex should be taken with care. On 13 individuals of I. perado, only one substitution (a transition) is observed in ITS. On the other hand, ITS sequences of most individuals of I. canariensis are polymorphic and few sequences are identical. The divergence between all ITS sequences observed in I. canariensis is much higher than between available ITS sequences of all other species investigated. The GC content is 57% for I. perado and from 45 to 62% for I. canariensis. Regarding their GC content, three groups of ITS are found in I. canariensis: a clade with 45% GC, a clade with 61% GC in average and several clades with 53–54% GC in average. The GC content of other investigated species range from 55 to 61% (Figure 1).
The ITS sequences of I. canariensis are distributed in two groups in the phylogram represented in Figure 1. One group forms a large clade conserved in the consensus tree but not sustained by bootstrap statistics. Another group forms a small clade (GC 45%) which branches variably: as indicated in Figure 1, or at the position indicated by the asterisk.
Secondary structure of ITS 2
Pattern of substitution
Substitution patterns of I. perado and I. canariensis ITS sequences from a reconstructed maximum likelihood ancestral ITS sequence of Ilex.
As expected for pseudogenes [10, 11], the observed rate of deamination-like substitutions at methylated cytosine sites (CpG and CpNpG sites) is higher than the expected rate of C -> T and G -> A substitutions at non-methylated sites for can 28B_45, can 39_3_53, can 44_4_54 and can 36_1_62 (Table 1). A chi-square homogeneity test  indicates that this is highly significant for can 28B_45 and can 36_1_53, which certainly represent pseudogenes.
Substitutions at conserved sites of the 5.8S rDNA
The alignment of fifty 5.8S sequences (modified from Muir et al. ) shows that 59 sites are conserved in vertebrates, invertebrates, fungi and plants and are expected to be functionally constrained. Substitutions observed at these sites would suggest non-functional pseudogenes. Contrarily to I. perado, many ITS sequences found in I. canariensis have substitutions at some of these conserved sites. Only the GC 61% clade of I. canariensis comprises non-substituted conserved 5.8S sites (sequences indicated in bold in Figure 1). Sequences of the GC 45% clade have 10–11 substitutions. Sequences of the GC 53–54% clades have 2 to 7 substitutions. Three sequences of the GC 61% clade have only one mutation, which may be PCR artifacts  and two sequences (can_20_1_58 and can_27_2_59) with a lower GC content (58 and 59%, respectively) have 2 mutations. Thus ITS sequences of I. canariensis having a GC content higher than 60% are expected to be functional genes, all other sequences with lower GC content are suspected to be pseudogenes.
Most ITS sequences of I. canariensis experienced frequent recombinations: in the entire ITS matrix of I perado and I. canariensis, the DnaSP program  detects 19 minimum possible recombination events (RM). From 0 to 8 minimum possible recombination events are calculated in the different ITS clusters (Figure 4A). No recombination was detected in I. perado and in clade 2 of I. canariensis. Clades 5 and 6 are highly recombined. An example of recombined ITS sequences of I. canariensis from clade 5 is shown in Figure 4B, where only informative sites are shown.
The LARD maximum-likelihood method  was applied to find the breakpoints in the alignment, which gave the highest likelihood under an evolutionary model incorporating recombination. Only 3 sequences can be analyzed with this program. Three ITS sequences shown in Figure 4 were submitted to LARD: can 90_6_54, can 25_2_54 and can 90_4_53. Two recombination points were located by the program (at the left of positions 242 and at the left of position 582) in accordance with the delimitation indicated in Figure 4. There is no recombination point between positions 392 and 455 for these particular 3 sequences.
Nuclear encoded plastid glutamine synthetase (nepGS) data
The high polymorphism of ITS sequences observed in I. canariensis is frequently reported for other plant groups [15, 16]. It might have several origins: an incomplete lineage sorting from ancestral polymorphism or an horizontal transfer (introgression) through inter-specific hybridization (alloploidy), both of them not resolved by concerted evolution. Before the discussion on the origin of this polymorphism, the characterization and the fate of these different ITS sequences will be first examined.
The genome of Ilex canariensis contains ITS pseudogenes
High polymorphism of ITS has been explained by the presence of divergent pseudogenes in Gossypium, Nicotiana, Tripsacum, Exospermum, Zygogonum, Zea , Quercus , Leucaena , Adinauclea, Haldina, Mitragyna  and others. Thus, this could also be the case for I. canariensis. Individual criteria are not sufficient to identify pseudogenes unambiguously  and different criteria were chosen: GC content, secondary structure of ITS 2, rate of substitution, pattern of substitution at methylated cytosine sites and substitutions at highly conserved sites of the 5.8S rDNA. ITS sequences with a GC content of 45% are unambiguously pseudogenes and satisfy to all other criteria. Moreover they have a large deletion in the ITS 1 region, which make these sequences certainly non-functional. The deletion allows an easy detection of this pseudogene on agarose gels and it is observed in many individuals of I. canariensis (Figure 2).
Other classes of ITS sequences with a GC content of 53–54% are also suspected to be pseudogenes by one or the other criteria but not by all of them, as expected regarding their relatively high GC content. For instance, some ITS 2 sequences of the GC 53–54% class still have a typical angiosperm secondary structure (for instance can_39_3_53, Figure 3), but have (1) an increased rate of nucleotide substitution, (2) deamination-like substitutions or (3) mutations at normally highly conserved 5.8S rDNA sites. Only the GC 61 % clade contains functional ITS sequences. Thus, it can be considered that the functional ITS GC content is 57 % for I. perado and above 60 %. for I. canariensis.
It is interesting to note that most I. canariensis individuals of the GC 61 % clade never have ITS pseudogenes in the GC 45 % or GC 53–54 % clades. This is probably because these individuals do not contain pseudogenes. For other individuals, a PCR selection for pseudogenes occurred, as reported for Nicotiana , in which ITS sequences with a weak secondary structure (pseudogenes) are preferentially used as templates. The inclusion of dimethylsulfoxide (DMSO) in PCR reactions [10, 19], but see , would allow amplification of functional ITS sequences in these individuals of I. canariensis.
In conclusion, the high divergence found in ITS sequences of I. canariensis with a GC content lower than 60% (clades 1, 2, 3, 4, 5 and 6) could be explained by a release of evolutionary constraint and a subsequent high rate of substitution. Indeed, ITS sequences have functional constraints in relation with the processing of the rRNA precursor producing the functional 18S, 26S and 5.8S subunits.
ITS sequences of Ilex canariensis are recombined
Evidence for recombination in divergent sequences is not obvious. It is difficult to recognize homoplasy generated by recombination from actual homoplasy (parallel history). Statistical methods (based on linkage desequilibrium, neutrality tests and substitution distribution along the locus) are still too rudimentary to precisely describe the recombination events in the set of ITS sequences found in I. canariensis. Moreover, recombinants could result from "jumping" PCR reaction [20–23], where prematurely terminated extension products can act as primer on paralogous templates. This has been shown on nepGS for Oxalis  and on four low-copy genes for Gossypium .
The minimum number of recombination events (RM) calculated with DnaSP  underestimates the total number of recombination events . Thus, there is no doubt that I. canariensis ITS sequences experienced intra-molecular recombinations (Figure 4). The factor RM has been also calculated for PCR products of each individual in order to detect possible jumping PCR artifacts. In few of them (specimens 20, 22, 24, 27, 28 and 38) recombinants have been detected in ITS sequences resulting from a unique PCR reaction (data not shown). This could be the result of jumping PCR. However most of them are multiple recombinants and not simple recombinants as it is expected in jumping PCR . As an example, the alignment represented in Figure 4 shows that specimen 90 comprises two different recombined ITS sequences resulting from the same PCR reaction, that could be the result of jumping PCR. DnaSP did not detect recombination between the four cloned ITS sequences of individual 90 because recombined fragments are paralogous sequences fragments found in other individuals. Moreover, the recombinants result from at least three crossover events and are suggested not PCR artifact. Thus, they represent true organismal intra-molecular recombinations.
The distribution of informative characters shown in Figure 4, as well as the use of programs PIST and LARD based on maximum-likelihood analyses, demonstrate unambiguously that sequences of clade 5 (Figure 4) experienced recombination events. This can not be generalized for other clades. Although DnaSP suggests recombination, an alignment demonstrating recombination, as for clade 5, was not possible for other clades, even with the help of PIST and LARD. This could be explained by the recent origin of the recombination events observed in clade 5 and by the fact that mutations did not yet obscured the recombined orthologous fragments. In this respect it is to be noticed that clade 5 shows much longer branches than other clades. This may indicate that, in clades with relatively shorter branches, mutations (or concerted evolution) did homogenize the recombined fragments, mimicking clonal divergence. Thus it can be considered that most clades also comprise recombined ITS sequences, as DnaSP suggests, but of more ancient origin than those of clade 5, and homogenized by mutation or concerted evolution.
Recombination in highly polymorphic ITS sequences seems a rule in plants. This is not surprising because the mechanisms of concerted evolution in rDNA arrays are based on crossing-over and gene conversion. It has been reported in Begonia , Microseris , Quercus , Amelanchier , Paeonia , Buddia, Gossypium, Nicotiana, Tripsacum , Armeria  and others.
In addition to the high rate of substitution of pseudogenes, at least some ITS sequences experienced recombination. This explains why the divergence between ITS sequences of I. canariensis is much higher than between ITS sequences of all other species investigated, knowing that, according to their GC content (see Figure 1), they all are potentially functional. This also explains the absence of a bootstrap support for a monophyletic clade of I. canariensis ITS sequences because of long branch problems due to accelerated rate of substitution and more certainly to recombination.
The origin of the ITS polymorphism in I. canariensis
Two evolutionary mechanisms could produce the observed ITS polymorphism: an ancestral polymorphism escaping lineage sorting or a past or recent introgression of an alien genotype escaping concerted evolution. Because of the influence of concerted evolution, ancestral polymorphism is not the most likely explanation of ITS polymorphism . On the other hand, a growing number of reports shows that ITS polymorphism is attributable to interspecific hybridization, although the parents are not always identifiable [15, 16].
Assuming that multiple ITS sequences found in I. canariensis are the result of experienced hybridization with another species, or an ancient polymorphism with incomplete sorting, the determination of the identity of the putative hybridizing species or the finding of genetic relationships of the putative polymorphism is not obvious. This is because ITS sequences enclosed in non-functional clusters have dramatically diverged from the putative functional sequences and are recombined. All available ITS sequences of 43 other species of Ilex, representing a good sampling of the genus [2, 3] were incorporated in the phylogenetic analysis, altogether with all ITS sequences found in I. perado and I. canariensis of Tenerife (Figure 1). Most functional (above 60% GC) and non-functional (53–54% GC) ITS clades of I. canariensis group together but with no bootsrap support. They group with an American lineage (I. brevicuspis, I. anomala, I. microdunta, I. integerrima, I. theezans, I. guianensis, I. brasiliensis and I. cassine). Only the GC 45% clade does not group with the bulk of I. canariensis ITS sequences. Its position is not defined and varies in the vicinity of a Eurasian lineage (I. latifolia, I. leucoclada, I. maximocziana, I. rugosa and I. perado). Thus data do not support a particular relationship of most I. canariensis ITS pseudogenes with another Ilex species, except for the pseudogenes with a GC content of 45%, that are frequently observed in I. canariensis.
In the case of hybridization involving the island species I. canariensis, the most probable candidate would be the sympatric species I. perado. It can not be ruled out however that the distribution of I. canariensis was much wider in the past [32, 33] and that this hybridization may have occurred with another unknown or extinct species of the Eurasian lineage represented here by I. latifolia, I. leucoclada, I. maximocziana, I. rugosa and I. perado. Pseudogene sequences (particularly the ITS sequences of clade GC 45%) being too divergent and of different nucleotide composition, the observed relationship of clade GC 45% with the group of species comprising I. perado is questionable because of possible spurious long branch attraction. However, the data of the nuclear encoded plastid glutamine synthetase (a nuclear single copy locus) are not conflicting with an introgression of I. perado in I. canariensis. All the eight polymorphic sites observed in I. canariensis always comprise one allele shared with I. perado. Another possibility is that these ITS pseudogenes represent a relictual ancestral polymorphism in the course of elimination by lineage sorting or concerted evolution. In fact ancestral polymorphism could also be the result of ancient introgressions. The data accumulated here do not allow a definitive conclusion.
If a putative cryptic hybridization between I. perado and I. canariensis is confirmed, the introgression would be unidirectional because ITS sequences of I. perado do not show any polymorphism. This situation is reminiscent of the unilateral hybridization observed between Begonia formosana and B. aptera, where on 60 ITS sequences analysed in natural or artificial hybrids, 58 sequences are clustering with the ovule donor B. formosana, and only 2 are found clustering with the pollen donor B. aptera . Unidirectional interspecific hybridization linked to unilateral incompatibility is frequently described in plants. However, this is not the only mechanism that can explain unidirectional hybridization. The flowering time of I. perado precedes the one of I. canariensis, thus the loading of still living I. perado pollen grains on young effective I. canariensis stigmates is more favored than the contrary. Moreover, there are much more male than female I. perado plants in Tenerife [34, 35]. These evidences could explain the proposed unidirectional introgression.
This study was undertaken with the aim to study and overcome the problem of ITS polymorphism found in many species of Ilex . Introgression [3, 4] and high polymorphism  have already been shown in several species of Ilex. Thus, precautionary measures should be taken when studying nuclear ITS sequences in the genus, particularly the search for recombinant and pseudogenes. Particular PCR conditions should be used [10, 19, 23]. Razafimandimbison and al.  were however unable to find PCR conditions to amplify a functional ITS sequence in 2 species of Rubiaceae. Amplified ITS sequences should be checked for function. A measure of the GC content (above 55% for a functional ITS sequence in Ilex is recommended. This study will probably make phylogenetic interpretations easier and will certainly help to the understanding of the complex evolutionary history of Ilex .
Thirteen individuals of I. perado ssp. platyphylla Webb. & Berth. and 27 individuals of I. canariensis Poir. were collected in Tenerife (Las Mercedez, Aqua Garcia and Aqua Mansa) and genomic DNA was extracted from dry leaves as previously reported .
The plastid atpB-rbcL spacer was sequenced for the 40 specimens of Tenerife according to Cuénoud et al. . In a first experiment all ITS sequences (ITS 1, 5.8S and ITS 2) were directly sequenced from the PCR fragment according to Manen et al. . All individuals of I. perado and five individuals of I. canariensis produced a perfectly readable ITS sequence with no polymorphism. On the other hand, 22 individuals of I. canariensis produced an unreadable highly polymorphic ITS sequence. Cloning in E. coli was necessary and four clones per individual were sequenced. For all specimens, the nuclear encoded plastid glutamine synthetase (nepGS) was amplified and sequenced according to Emshwiller and Doyle . Polymorphisms were also observed in most individuals of I. canariensis, but as indels are not involved in the polymorphism, sequences were readable and polymorphic sites were coded according to the international code for nucleotide polymorphism (see Figure 5). Sequences are deposited at GenBank (atpB-rbcL spacer: AJ786512-AJ786551; ITS: AJ786413-AJ786504; nepGS: AJ809595-AJ809628).
The ITS sequences of I. canariensis and I. perado were aligned with ITS sequences of 43 species of Ilex previously studied in Manen et al.  and Phyllonoma and Helwingia, the closest relatives of the genus Ilex (for the ITS matrix see additional file 1).
Plastid atpB-rbcL spacer and nuclear nepGS matrices only comprise the sequences found in Tenerife for I. perado and I. canariensis.
Maximum parsimony trees were calculated from the atpB-rbcL spacer matrix, the ITS matrix and the nepGS matrix, using PAUP 4.0b10  (heuristic search, TBR branch swapping with 10 random additions of sequences, only keeping the first 100 most parsimonious trees). Bootstrap statistics of the ITS tree were calculated from 1000 replications with the same method, except that the first 10 most parsimonious trees were kept.
Secondary structure of ITS 2
The secondary structure of ITS 2 was investigated using the minimum free-energy program Mfold , which has the advantage to provide sub-optimal folding. Sequences were constrained to force the pairing of the 5'-end of the 26S and the 3'-end of the 5.8S regions according to the results of Hershkovitz and Zimmer .
Pattern of substitution of ITS
In order to examine the pattern of substitution of ITS sequences found in I. perado and I. canariensis, these sequences were compared with a reconstructed Ilex ancestral sequence. This ancestral sequence was determined by maximum likelihood [PAUP 4.0b10 with base frequencies, ti/tv ratio, proportion of invariable sites and gamma shape parameter estimated under the HKY model (Hasegawa et al. 1985) allowing for different rate of transitions and transversions as well as unequal base frequencies] using the unique maximum parsimony tree obtained from the ITS matrix of Manen et al.  based on 45 species of Ilex and Helwingia and Phyllonoma as outgroup. For comparisons, the number of substitutions (Kimura 2-parameter) was calculated from this ancestral sequence for I. perado and I. canariensis ITS sequences using PAUP.
The high frequency of deamination-like substitutions (C -> T and G -> A at CpG and CpNpG sites) is typical for pseudogenes and was also calculated from the reconstructed Ilex ancestral sequence for I. perado and I. canariensis ITS sequences, according to Buckler et al.  and Muir et al. . This frequency was compared with the frequency of C -> T and G -> A substitutions at non-methylated sites.
Substitutions at conserved sites of the ribosomal 5.8 S subunit
Based on the alignment of fifty 5.8S sequences including vertebrates, invertebrates, fungi and plants (modified from Muir et al. ), 59 totally conserved sites were determined. The number of substitutions observed at these invariant sites in all I. perado and I. canariensis 5.8S sequences was calculated.
Detection of recombinations
The minimum number of recombination events RM  was calculated using the DnaSP program  in the entire ITS matrix of I. perado and I. canariensis and in different ITS clusters observed in I. canariensis. This program is based on linkage desequilibrium, neutrality tests and substitution distribution along the locus. ITS sequences showing strong evidence of recombination, as detected by DnaSP, were submitted to PIST  to calculate the probability of recombination by maximum likelihood: the tree score of these sequences is compared with the tree scores of 1000 simulated clonal sequences along the specified tree under the specified model of evolution. The score of recombined sequences will tend to have larger score than the simulated clonal sequences because richer in conflicting phylogenetic information. A maximum likelihood tree of selected I. canariensis ITS sequences (see results) was constructed with base frequencies, ti/tv ratio, proportion of invariable sites and gamma shape parameter estimated under the HKY model  and these parameters were used to calculate the tree scores of simulated clonal sequences.
The recombination points of three selected sequences (see results) showing evidence of recombination were calculated by maximum likelihood using LARD  with the HKY model . The program calculates a maximum likelihood unrooted tree of 3 sequences and searches for a tree with a better score assuming a recombination point in the input sequences. After the calculation of one recombination point, the sequence alignment was truncated at this point to search for other potential recombination points.
I wish to thank Yamama Naciri-Graven (CJB) for statistical evaluation of the pattern of substitutions, Gabrielle Barriera and Pierre-André Loizeau (CJB) for the confirmation of species determination and Gisèle Vuille-dit-Bille for technical assistance. This work was supported by the "Conservatoire et Jardin Botaniques" and the University of Geneva.
- Powell M, Savolainen V, Cuénoud P, Manen JF, Andrews S: The montain holly (Nemopanthus mucronatus : Aquifoliaceae) revisited with molecular data. Kew Bull. 2000, 55: 341-347.View ArticleGoogle Scholar
- Cuénoud P, Del Pero Martinez MA, Loizeau PA, Spichiger R, Andrews S, Manen JF: Molecular phylogeny and biogeography of the genus Ilex L. (Aquifoliaceae). Ann Bot. 2000, 85: 111-122. 10.1006/anbo.1999.1003.View ArticleGoogle Scholar
- Manen JF, Boulter MC, Naciri-Graven Y: The complex history of the genus Ilex L. (Aquifoliaceae): evidence from the comparison of plastid and nuclear DNA sequences and from fossil data. Plant Syst Evol. 2002, 235: 79-98. 10.1007/s00606-002-0225-x.View ArticleGoogle Scholar
- Setogushi H, Watanabe I: Intersectional gene flow between insular endemics of Ilex (Aquifoliaceae) on the Bonin Islands and the Ryukyu Islands. Am J Bot. 2000, 87: 793-810.View ArticleGoogle Scholar
- Torimaru T, Tomaru N, Nishimura N, Yamamoto S: Clonal diversity and genetic differentiation in Ilex leucoclada M. patches in an old-growth beech forest. Mol Ecol. 2003, 12: 809-818.View ArticlePubMedGoogle Scholar
- Index to Plant Chromosome Numbers. [http://mobot.mobot.org/W3T/Search/ipcn.html]
- Hershkovitz MA, Zimmer EA: Conservation patterns in angiosperm rDNA ITS 2 sequences. Nucl Acids Res. 1996, 24: 2857-2867. 10.1093/nar/24.15.2857.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucl Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Sokal RR, Rohlf FJ: Biometry. The principles and practice of statistics in Biological research. 1995, WH Freeman and Company. New York, 3Google Scholar
- Buckler ES, Ippolito A, Holtsford TP: The evolution of ribosomal DNA: Divergent paralogues and phylogenetic implications. Genetics. 1997, 145: 821-832.PubMedGoogle Scholar
- Muir G, Fleming CC, Schlötterer C: Three divergent rDNA clusters predate the species divergence in Quercus petraea (Matt.) Leibl. and Quercus robur L. Mol Biol Evol. 2001, 18: 112-119.View ArticlePubMedGoogle Scholar
- Rozas J, Rozas R: DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics. 1999, 15: 174-175. 10.1093/bioinformatics/15.2.174.View ArticlePubMedGoogle Scholar
- Worobey M: A novel approach to detecting and measuring recombination: new insights into evolution in viruses, bacteria, and mitochodria. Mol Biol Evol. 2001, 18: 1425-1434.View ArticlePubMedGoogle Scholar
- Holmes EC, Worobey M, Rambaut A: Phylogenetic evidence for recombination in dengue virus. Mol Biol Evol. 1999, 16: 405-409.View ArticlePubMedGoogle Scholar
- Bailey CD, Carr TG, Harris SA, Hughes CE: Characterization of angiosperm nrDNA polymorphism, paralogy, and pseudogenes. Mol Phylogenet Evol. 2003, 29: 435-455. 10.1016/j.ympev.2003.08.021.View ArticlePubMedGoogle Scholar
- Alvarez I, Wendel JF: Ribosomal ITS sequences and plant phylogenetic inference. Mol Phyl Evol. 2003, 29: 417-434. 10.1016/S1055-7903(03)00208-2.View ArticleGoogle Scholar
- Hughes CE, Bailey CD, Harris SA: Divergent and reticulate species relationships in Leucaena (Fabaceae) inferred from multiple data sources: insights into polyploid origins and nrDNA polymorphism. Am J Bot. 2002, 89: 1057-1073.View ArticlePubMedGoogle Scholar
- Razafimandimbison SG, Kellogg EA, Bremer B: Recent origin and phylogenetic utility of divergent ITS putative pseudogenes: A case study from Naucleeae (Rubiaceae). Syst Biol. 2004, 53: 177-192. 10.1080/10635150490423278.View ArticlePubMedGoogle Scholar
- Varadaraj K, Skinner DM: Denaturants or cosolvents improve the specificity of PCR amplification of a G+C-rich DNA using genetically-engineered DNA-polymerases. Gene. 1994, 140: 1-5. 10.1016/0378-1119(94)90723-4.View ArticlePubMedGoogle Scholar
- Pääbo S, Irwin DM, Wilson AC: DNA damage promotes jumping between templates during enzymatic amplification. J Biol Chem. 1990, 265: 4718-4721.PubMedGoogle Scholar
- Jansen R, Ledley FD: Disruption of phase during PCR amplification and cloning of heterozygous target sequences. Nucl Acids Res. 1990, 18: 5153-5156.PubMed CentralView ArticlePubMedGoogle Scholar
- Bradley RD, Hillis DM: Recombinant DNA sequences generated by PCR amplification. Mol Biol Evol. 1997, 14: 592-593.View ArticlePubMedGoogle Scholar
- Thompson JR, Marcelino LA, Polz MF: Heteroduplexes in mixed-template amplifications: formation, consequence and elimination by "reconditioning" PCR. Nucl Acids Res. 2002, 30: 2083-2088. 10.1093/nar/30.9.2083.PubMed CentralView ArticlePubMedGoogle Scholar
- Emshwiller E, Doyle JJ: Origins of domestication and polyploidy in oca (Oxalis tuberosa: Oxalidaceae). 2. Chloroplast-expressed glutamine synthetase data. Amer J Bot. 2002, 89: 1042-1056.View ArticleGoogle Scholar
- Cronn R, Cedroni M, Haselkorn T, Grover C, Wendel JF: PCR-mediated recombination in amplification products derived from polyploid cotton. Theor Appl Genet. 2002, 104: 482-489. 10.1007/s001220100741.View ArticlePubMedGoogle Scholar
- Hudson RR, Kaplan NL: Statistical properties of the number of recombination events in the history of a sample of DNA sequences. Genetics. 1985, 111: 147-164.PubMed CentralPubMedGoogle Scholar
- Chiang TY, Hong KH, Peng CI: Experimental hybridization reveals biased inheritance of the internal transcribed spacer of the nuclear ribosomal DNA in Begonia x taipeiensis. J Plant Res. 2001, 114: 343-351.View ArticleGoogle Scholar
- Roelofs D, Van Velzen J, Kuperus P, Bachmann K: Molecular evidence for an extinct parent of the tetraploid species Microseris acuminata and M. campestris (Asteraceae, Lactuceae). Mol Ecol. 1997, 6: 641-649. 10.1046/j.1365-294X.1997.00225.x.View ArticlePubMedGoogle Scholar
- Campbell CS, Wojciechowski MF, Baldwin BG, Alice LA, Donoghue MJ: Persistent nuclear ribosomal DNA sequence polymorphism in the Amelanchier agamic complex (Rosaceae). Mol Biol Evol. 1997, 14: 81-90.View ArticlePubMedGoogle Scholar
- Sang T, Crawford DJ, Stuessy TF: Documentation of reticulate evolution in peonies (Paeonia) using internal transcribed spacer sequences of nuclear ribosomal DNA: implications for biogeography and concerted evolution. Proc Natl Acad Sci USA. 1995, 92: 6813-6817.PubMed CentralView ArticlePubMedGoogle Scholar
- Nieto Feliner G, Gutiérrez Larenza B, Fuertes Aguilar J: Fine-scale geographical structure, intra-individual polymorphism and recombination in nuclear ribosomal internal transcribed spacers in Armeria (Plumbaginaceae). Ann Bot. 2004, 93: 189-200. 10.1093/aob/mch027.PubMed CentralView ArticlePubMedGoogle Scholar
- Sunding P: Origins of the macaronesian flora. In: Plants and islands. Edited by: Bramwell D. 1979, Academic Press Inc. London, 13-40.Google Scholar
- Bramwell D, Bramwell ZI: Wild Flowers of the Canary Islands. 1974, Stanley Thornes (Publishers) Ltd. London and BurfordGoogle Scholar
- Kunkel G: The Ilex complex in the Canary Islands and Madeira. Cuad Bot Canar. 1977, 28: 17-29.Google Scholar
- Andrews S: Hollies of the Canary Islands. Bull Holly Soc Amer. 1984, 21: 1-10.Google Scholar
- Emshwiller E, Doyles JJ: Chloroplast-expressed glutamine synthetase (ncpGS): potential utility for phylogenetic studies with an example from Oxalis (Oxalidaceae). Mol Phyl Evol. 1999, 12: 310-319. 10.1006/mpev.1999.0613.View ArticleGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. Sinauer Associates, Sunderland, Massachusetts. 2003Google Scholar
- Hasegawa M, Kishino M, Yano T: A hidden markov model approach to variation among sites in rate of evolution. Mol Biol Evol. 1996, 22: 160-174.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.