Ecological niche comparison and molecular phylogeny segregate the invasive moss species Campylopus introflexus (Leucobryaceae,Bryophyta) from its closest relatives

The delimitation of the invasive moss species Campylopus introflexus from its closest relative, Campylopus pilifer, has been long debated based on morphology. Previous molecular phylogenetic reconstructions based on the nuclear ribosomal internal transcribed spacers (ITS) 1 and 2 showed that C. pilifer is split into an Old World and a New World lineage, but remained partly inconclusive concerning the relationships between these two clades and C. introflexus. Analyses of an extended ITS dataset displayed statistically supported incongruence between ITS1 and ITS2. ITS1 separates the New World clade of C. pilifer from a clade comprising C. introflexus and the Old World C. pilifer. Ancestral state reconstruction showed that this topology is morphologically supported by differences in the height of the dorsal costal lamellae in leaf cross‐section (despite some overlap). ITS2, in contrast, supports the current morphological species concept, i.e., separating C. introflexus from C. pilifer, which is morphologically supported by the orientation of the hyaline hair point at leaf apex as well as costal lamellae height. Re‐analysis of published and newly generated plastid atpB‐rbcL spacer sequences supported the three ITS lineages. Ecological niche modeling proved a useful approach and showed that all three molecular lineages occupy distinct environmental spaces that are similar, but undoubtedly not equivalent. In line with the ITS1 topology, the C. pilifer lineage from the New World occupies the most distinct environmental niche, whereas the niches of Old World C. pilifer and C. introflexus are very similar. Taking the inferences from ecological niche comparisons, phylogenetics, and morphology together, we conclude that all three molecular lineages represent different taxa that should be recognized as independent species, viz. C. introflexus, C. pilifer (Old World clade), and the reinstated C. lamellatus Mont. (New World clade).     

Phylogeny of Impatiens (Balsaminaceae): integrating molecular and morphological evidence into a new classification

Impatiens L. is one of the largest angiosperm genera, containing over 1000 species, and is notorious for its taxonomic difficulty. Here, we present, to our knowledge, the most comprehensive phylogenetic analysis of the genus to date based on a total evidence approach. Forty‐six morphological characters, mainly obtained from our own investigations, are combined with sequence data from three genetic regions, including nuclear ribosomal ITS and plastid atpB‐rbcL and trnL‐F. We include 150 Impatiens species representing all clades recovered by previous phylogenetic analyses as well as three outgroups. Maximum‐parsimony and Bayesian inference methods were used to infer phylogenetic relationships. Our analyses concur with previous studies, but in most cases provide stronger support. Impatiens splits into two major clades. For the first time, we report that species with three‐colpate pollen and four carpels form a monophyletic group (clade I). Within clade II, seven well‐supported subclades are recognized. Within this phylogenetic framework, character evolution is reconstructed, and diagnostic morphological characters for different clades and subclades are identified and discussed. Based on both morphological and molecular evidence, a new classification outline is presented, in which Impatiens is divided into two subgenera, subgen. Clavicarpa and subgen. Impatiens; the latter is further subdivided into seven sections.     

Phylogeny,biogeography, and morphological ancestral character reconstruction in the Mediterranean genus Fumana (Cistaceae)

Fumana is a diverse genus of the Cistaceae family, consisting of 21 currently accepted species. In this study, nuclear (ITS) and plastid (matK, trnT‐L) molecular markers were used to reconstruct the phylogeny and to estimate divergence times, including 19 species of Fumana. Phylogenetic analyses (Bayesian Inference, Maximum Parsimony and Maximum Likelihood) confirmed the monophyly of Fumana and did not support the infrageneric divisions previously established. The results support four main clades that group species that differ in vegetative and reproductive characters. Given the impossibility to define morphological characters common to all species within the clades, our proposal is to reject infrageneric divisions. Molecular dating and ancestral area analyses provide evidence for a Miocene diversification of the genus in the north‐western Mediterranean. Ancestral state reconstructions revealed ancestral character states for some traits related to xeric and arid habitats, suggesting a preadaptation to the Mediterranean climate.     

Comparative phylogeographies of six species of hinged terrapins (Pelusios spp.) reveal discordant patterns and unexpected differentiation in the P. castaneus/P. chapini complex and P. rhodesianus

Using up to 2117 bp of mitochondrial DNA and up to 2012 bp of nuclear DNA, we analysed phylogeographic differentiation of six widely distributed species of African hinged terrapins (Pelusios spp.) representing different habitat types. Two taxa each live in savannahs or in forests and mesic savannahs, respectively, and the remaining two species occur in intermediate habitats. The species living in forests and mesic savannahs do not enter dry savannahs, whereas the savannah species may occur in dry and wet savannahs and even in semi‐arid steppe regions. We found no obvious correlation between habitat type and phylogeographic pattern: one savannah species (P. rhodesianus) shows phylogeographic structure, i.e. pronounced genetic differences among geographically distinct populations, and the other (P. nanus) not. One species inhabiting forests and mesic savannahs (P. carinatus) has phylogeographic structure, the other (P. gabonensis) not. The same pattern is true for the two ecologically intermediate species, with phylogeographic structure present in P. castaneus and absent in P. chapini. Nuclear evidence suggests that the latter two taxa with abutting and partially overlapping ranges are distinct, while mtDNA is only weakly differentiated. Pelusios castaneus shows pronounced phylogeographic structure, which could reflect Pleistocene range interruptions correlated with the fluctuating forest cover in West and Central Africa. Our results do not support the recognition of an extinct subspecies of P. castaneus for the Seychelles. Pelusios carinatus contains two well supported clades, which are separated by the Congo River. This species is closely related to P. rhodesianus, a taxon consisting of two deeply divergent mitochondrial clades. One of these clades is paraphyletic with respect to P. carinatus, but the two clades of P. rhodesianus are not differentiated in the studied nuclear markers and, again, paraphyletic with respect to P. carinatus. Using mtDNA sequences from the type material of P. rhodesianus, we were able to allocate this name to one of the two clades. However, owing to the confusing relationships of P. rhodesianus and P. carinatus, we refrain from taxonomic decisions.     

Molecular phylogenetics of slit‐faced bats (Chiroptera: Nycteridae) reveal deeply divergent African lineages

The bat family Nycteridae contains only the genus Nycteris, which comprises 13 currently recognized species from Africa and the Arabian Peninsula, one species from Madagascar, and two species restricted to Malaysia and Indonesia in South‐East Asia. We investigated genetic variation, clade membership, and phylogenetic relationships in Nycteridae with broad sampling across Africa for most clades. We sequenced mitochondrial cytochrome b (cytb) and four independent nuclear introns (2,166 bp) from 253 individuals. Although our samples did not include all recognized species, we recovered at least 16 deeply divergent monophyletic lineages using independent mitochondrial and multilocus nuclear datasets in both gene tree and species tree analyses. Mean pairwise uncorrected genetic distances among species‐ranked Nycteris clades (17% for cytb and 4% for concatenated introns) suggest high levels of phylogenetic diversity in Nycteridae. We found a large number of designated clades whose members are distributed wholly or partly in East Africa (10 of 16 clades), indicating that Nycteris diversity has been historically underestimated and raising the possibility that additional unsampled and/or undescribed Nycteris species occur in more poorly sampled Central and West Africa. Well‐resolved mitochondrial, concatenated nuclear, and species trees strongly supported African ancestry for SE Asian species. Species tree analyses strongly support two deeply diverged subclades that have not previously been recognized, and these clades may warrant recognition as subgenera. Our analyses also strongly support four traditionally recognized species groups of Nycteris. Mitonuclear discordance regarding geographic population structure in Nycteris thebaica appears to result from male‐biased dispersal in this species. Our analyses, almost wholly based on museum voucher specimens, serve to identify species‐rank clades that can be tested with independent datasets, such as morphology, vocalizations, distributions, and ectoparasites. Our analyses highlight the need for a comprehensive revision of Nycteridae.     

Internal Transcribed Spacer Sequence Analysis of Puccinia helianthi Schw. and Its Application in Detection of Sunflower Rust

Two basidiomycete‐specific primers ITS1‐F and ITS4‐B were used in identification of the genus Puccinia. The primers showed good specificity for the genus with an 816‐bp product that was amplified exclusively. Twenty sequences of internal transcribed spacer (ITS) regions of Puccinia helianthi isolates from China remain unchanged. The whole ITS length (including ITS1 sequence 194 bp, 5.8S rRNA gene 156 bp, ITS2 sequence 206 bp) was 556 bp. By comparing the aligned ITS sequences of several Puccinia isolates from China, Spain and the United States, ITS homogeneity among these sunflower rust isolates was >99%. Genetic homology and phylogeny of P. helianthi with other Puccinia spp. was investigated. Nineteen sequences of rDNA ITS1 and ITS2 were determined and used as phylogenetic markers. Phylogenetic analysis of ITS regions showed that Puccinia spp. of sunflower was clustered in one clade with P. komarovii and P. violae, divergent from Puccinia spp. of Chrysanthemum, P. tenaceti of tansy (Tanacetum vulgare) and Puccina spp. of big sagebrush (Artemisia tridentate) indicating sunflower rust had distant phylogenetic relationships with other Compositae rusts. With the specified primers SR‐1 and SR‐2, either from purified urediniospores or symptomless (but infected) sunflower leaves could be examined specifically. Therefore, results of this study help in detection and polygenetic study of rust fungi occurring on sunflower.     

Morphological and genetic diversity of Beaufort Sea diatoms with high contributions from the Chaetoceros neogracilis species complex

Seventy‐five diatom strains isolated from the Beaufort Sea (Canadian Arctic) in the summer of 2009 were characterized by light and electron microscopy (SEM and TEM), as well as 18S and 28S rRNA gene sequencing. These strains group into 20 genotypes and 17 morphotypes and are affiliated with the genera Arcocellulus, Attheya, Chaetoceros, Cylindrotheca, Eucampia, Nitzschia, Porosira, Pseudo‐nitzschia, Shionodiscus, Thalassiosira, and Synedropsis. Most of the species have a distribution confined to the northern/polar area. Chaetoceros neogracilis and Chaetoceros gelidus were the most represented taxa. Strains of C. neogracilis were morphologically similar and shared identical 18S rRNA gene sequences, but belonged to four distinct genetic clades based on 28S rRNA, ITS‐1 and ITS‐2 phylogenies. Secondary structure prediction revealed that these four clades differ in hemi‐compensatory base changes (HCBCs) in paired positions of the ITS‐2, suggesting their inability to interbreed. Reproductively isolated C. neogracilis genotypes can thus co‐occur in summer phytoplankton communities in the Beaufort Sea. C. neogracilis generally occurred as single cells but also formed short colonies. It is phylogenetically distinct from an Antarctic species, erroneously identified in some previous studies as C. neogracilis, but named here as Chaetoceros sp. This work provides taxonomically validated sequences for 20 Arctic diatom taxa, which will facilitate future metabarcoding studies on phytoplankton in this region.     

Identification of species and materia medica within Angelica L. (Umbelliferae) based on phylogeny inferred from DNA barcodes

DNA barcodes have been increasingly used in authentication of medicinal plants, while their wide application in materia medica is limited in their accuracy due to incomplete sampling of species and absence of identification for materia medica. In this study, 95 leaf accessions of 23 species (including one variety) and materia medica of three Pharmacopoeia‐recorded species of Angelica in China were collected to evaluate the effectiveness of four DNA barcodes (rbcL, matK, trnH‐psbA and ITS). Our results showed that ITS provided the best discriminatory power by resolving 17 species as monophyletic lineages without shared alleles and exhibited the largest barcoding gap among the four single barcodes. The phylogenetic analysis of ITS showed that Levisticum officinale and Angelica sinensis were sister taxa, which indicates that L. officinale should be considered as a species of Angelica. The combination of ITS + rbcL + matK + trnH‐psbA performed slight better discriminatory power than ITS, recovering 23 species without shared alleles and 19 species as monophyletic clades in ML tree. Authentication of materia medica using ITS revealed that the decoction pieces of A. sinensis and A. biserrata were partially adulterated with those of L. officinale, and the temperature around 80 °C processing A. dahurica decoction pieces obviously reduced the efficiency of PCR and sequencing. The examination of two cultivated varieties of A. dahurica from different localities indicated that the four DNA barcodes are inefficient for discriminating geographical authenticity of conspecific materia medica. This study provides an empirical paradigm in identification of medicinal plants and their materia medica using DNA barcodes.     

Examination of prezygotic and postzygotic isolating barriers in tropical Ulva (Ulvophyceae,Chlorophyta): evidence for ongoing speciation

Phylogenetic clades based on DNA sequences such as the chloroplast rbcL gene and the nuclear ITS region are frequently used to delimit algal species. However, these molecular markers cannot accurately delimit boundaries among some Ulva species. Although Ulva reticulata and Ulva ohnoi occasionally bloom in tropical to warm‐temperate regions and are clearly distinguishable by their reticulate or plain blade morphology, they have few or no sequence divergences in these molecular markers and form a monophyletic clade. In this study, to clarify the speciation and species delimitation in the U. reticulata‐ohnoi complex clade, reproductive relationships among several sexual strains from the Philippines and Japan including offspring that originated from the type specimen of U. ohnoi were examined by culturing and hybridization in addition to the ITS‐based analysis. As a result, both prezygotic and postzygotic reproductive isolation were revealed to occur between genetically perforated U. reticulata and imperforate U. ohnoi. They were also separated on the basis of sequence analysis of the ITS region. That strongly supports that the two taxa are independent biological species. Although no prezygotic barrier among the Philippine and Japanese strains of U. reticulata was observed, unexpectedly zoospores produced by hybrid sporophytes in some of their combinations mostly failed to develop, indicating partial formation of a postzygotic barrier despite a 0.2% divergence in the ITS sequence. These findings suggest speciation is still ongoing in U. reticulata.     

Species concept and nomenclatural changes within the genera Elliptochloris and Pseudochlorella (Trebouxiophyceae) based on an integrative approach

The genera Elliptochloris and Pseudochlorella were erected for Chlorella‐like green algae producing two types of autospores and cell packages, respectively. Both genera are widely distributed in different soil habitats, either as free living or as photobionts of lichens. The species of these genera are often difficult to identify because of the high phenotypic plasticity and occasional lack of characteristic features. The taxonomic and nomenclatural status of these species, therefore, remains unclear. In this study, 34 strains were investigated using an integrative approach. Phylogenetic analyses demonstrated that the isolates belong to two independent lineages of the Trebouxiophyceae (Elliptochloris and Prasiola clades) and confirmed that the genera are not closely related. The comparison of morphology, molecular phylogeny, and analyses of secondary structures of SSU and ITS rDNA sequences revealed that all of the strains belong to three genera: Elliptochloris, Pseudochlorella, and Edaphochlorella. As a consequence of the taxonomic revisions, we propose two new combinations (Elliptochloris antarctica and Pseudochlorella signiensis) and validate Elliptochloris reniformis, which is invalidly described according to the International Code for Nomenclature (ICN), by designating a holotype. To reflect the high phenotypic plasticity of P. signiensis, two new varieties were described: P. signiensis var. magna and P. signiensis var. communis. Chlorella mirabilis was not closely related to any of these genera and was, therefore, transferred to the new genus Edaphochlorella. All of the taxonomic changes were highly supported by all phylogenetic analyses and were confirmed by the ITS‐2 Barcodes using the ITS‐2/CBC approach.     

Phylogeny and species delineation in the marine diatom Pseudo‐nitzschia (Bacillariophyta) using cox1, LSU,and ITS2 rRNA genes: A perspective in character evolution

Analyses of the mitochondrial cox1, the nuclear‐encoded large subunit (LSU), and the internal transcribed spacer 2 (ITS2) RNA coding region of Pseudo‐nitzschia revealed that the P. pseudodelicatissima complex can be phylogenetically grouped into three distinct clades (Groups I–III), while the P. delicatissima complex forms another distinct clade (Group IV) in both the LSU and ITS2 phylogenetic trees. It was elucidated that comprehensive taxon sampling (sampling of sequences), selection of appropriate target genes and outgroup, and alignment strategies influenced the phylogenetic accuracy. Based on the genetic divergence, ITS2 resulted in the most resolved trees, followed by cox1 and LSU. The morphological characters available for Pseudo‐nitzschia, although limited in number, were overall in agreement with the phylogenies when mapped onto the ITS2 tree. Information on the presence/absence of a central nodule, number of rows of poroids in each stria, and of sectors dividing the poroids mapped onto the ITS2 tree revealed the evolution of the recently diverged species. The morphologically based species complexes showed evolutionary relevance in agreement with molecular phylogeny inferred from ITS2 sequence–structure data. The data set of the hypervariable region of ITS2 improved the phylogenetic inference compared to the cox1 and LSU data sets. The taxonomic status of P. cuspidata and P. pseudodelicatissima requires further elucidation.     

Differential Organ Distribution,Pathogenicity and Benomyl Sensitivity of Colletotrichum spp. from Blackberry Plants in Northern Colombia

Blackberry anthracnose, caused by Colletotrichum spp., is an important disease of cultivated blackberry in the world. In Colombia, it is the number one limiting factor for commercial production. This study was conducted to determine the species of Colletotrichum infecting blackberry plants as well as the organ distribution, pathogenicity and response to benomyl of the isolated strains. Sixty isolates from stems (n = 20), thorns (n = 20) and inflorescences (n = 20) were identified as Colletotrichum acutatum and Colletotrichum gloeosporioides by a species‐specific polymerase chain reaction (PCR). Both Colletotrichum species were found in the same plant but on different organs. Colletotrichum gloeosporioides species predominated in thorn lesions (n = 16) and C. acutatum in stems (n = 15) and inflorescence (n = 15). Pathogenicity assays on detached blackberry organs demonstrated differences between the two species with an average period of lesion development of 8.7 days for C. gloeosporioides and 10.3 days for C. acutatum. Wound inoculated organs had 90% disease development compared to 17.5% in non‐wounded. All C. acutatum isolates (n = 34) were benomyl tolerant, whereas C. gloeosporioides isolates (n = 26) were 30.7% sensitive and 69.2% moderately tolerant. Phylogenetic analysis with ITS sequences of a subset of 18 strains showed that strains classified as C. gloeosporioides had 100% identity to Colletotrichum kahawae, which belongs to the C. gloeosporioides species complex, whereas C. acutatum strains clustered into two different groups, with high similarity to the A2 and the A4 molecular groups. These data demonstrate for the first time the differential distribution of both species complexes in blackberry plant organs and further clarifies the taxonomy of the strains.     

Deep sympatric mtDNA divergence in the autumnal moth (Epirrita autumnata)

Deep sympatric intraspecific divergence in mtDNA may reflect cryptic species or formerly distinct lineages in the process of remerging. Preliminary results from DNA barcoding of Scandinavian butterflies and moths showed high intraspecific sequence variation in the autumnal moth, Epirrita autumnata. In this study, specimens from different localities in Norway and some samples from Finland and Scotland, with two congeneric species as outgroups, were sequenced with mitochondrial and nuclear markers to resolve the discrepancy found between mtDNA divergence and present species‐level taxonomy. We found five COI sub‐clades within the E. autumnata complex, most of which were sympatric and with little geographic structure. Nuclear markers (ITS2 and Wingless) showed little variation and gave no indications that E. autumnata comprises more than one species. The samples were screened with primers for Wolbachia outer surface gene (wsp) and 12% of the samples tested positive. Two Wolbachia strains were associated with different mtDNA sub‐clades within E. autumnata, which may indicate indirect selection/selective sweeps on haplotypes. Our results demonstrate that deep mtDNA divergences are not synonymous with cryptic speciation and this has important implications for the use of mtDNA in species delimitation, like in DNA barcoding.     

Phylogenetic analysis of eastern Asian and eastern North American disjunct Lespedeza (Fabaceae) inferred from nuclear ribosomal ITS and plastid region sequences

Lespedeza (tribe Desmodieae, Fabaceae) follows a disjunct distribution in eastern Asia and eastern North America. Phylogenetic relationships among its species and related taxa were inferred from nuclear ribosomal internal transcribed spacer (ITS) and plastid sequences (trnH‐psbA, psbK‐psbI, trnK‐matK and rpoC1). We examined 35 species of Lespedeza, two of Kummerowia and one of Campylotropis, the sole constituents of the Lespedeza group. An analysis of these data revealed that the genus Campylotropis is sister to the other two genera. However, we were unable to resolve the relationships between Kummerowia and Lespedeza in the strict consensus trees of parsimony analyses based on plastid and combined DNA data. In the genus Lespedeza, the Old World subgenus Macrolespedeza is monophyletic, whereas the transcontinental subgenus Lespedeza is paraphyletic. Monophyly of eastern Asian species and of North American species is strongly supported. Although inconsistent with the traditional classification, this phylogenetic finding is consistent with seedling morphology. Three subgroups recognized in subgenus Macrolespedeza were unresolved in our phylogenetic trees. An incongruence length difference (ILD) test indicated that the two partitions (nuclear ITS and plastid sequences) were significantly incongruent, perhaps because of hybridization between species in Lespedeza. Most of the primary clades of tribe Desmodieae are Asian, implying that the relatively few New World ones, such as those in Lespedeza, are more recently derived from Asia. © 2010 The Linnean Society of London, Botanical Journal of the Linnean Society, 2010, 164 , 221–235.     

Comparative analyses of plastid and AFLP data suggest different colonization history and asymmetric hybridization between Betula pubescens and B. nana

Birches (Betula spp.) hybridize readily, confounding genetic signatures of refugial isolation and postglacial migration. We aimed to distinguish hybridization from range‐shift processes in the two widespread and cold‐adapted species B. nana and B. pubescens, previously shown to share a similarly east–west‐structured variation in plastid DNA (pDNA). We sampled the two species throughout their ranges and included reference samples of five other Betula species and putative hybrids. We analysed 901 individual plants using mainly nuclear high‐resolution markers (amplified fragment length polymorphisms; AFLPs); a subset of 64 plants was also sequenced for two pDNA regions. Whereas the pDNA variation as expected was largely shared between B. nana and B. pubescens, the two species were distinctly differentiated at AFLP loci. In B. nana, both the AFLP and pDNA results corroborated the former pDNA‐based hypothesis that it expanded from at least two major refugia in Eurasia, one south of and one east of the North European ice sheets. In contrast, B. pubescens showed a striking lack of geographic structuring of its AFLP variation. We identified a weak but significant increase in nuclear (AFLP) gene flow from B. nana into B. pubescens with increasing latitude, suggesting hybridization has been most frequent at the postglacial expansion front of B. pubescens and that hybrids mainly backcrossed to B. pubescens. Incongruence between pDNA and AFLP variation in B. pubescens can be explained by efficient expansion from a single large refugium combined with leading‐edge hybridization and plastid capture from B. nana during colonization of new territory already occupied by this more cold‐tolerant species.     

Plastid DNA sequence data help to clarify phylogenetic relationships and reticulate evolution in Lycoris (Amaryllidaceae)

The temperate East Asian genus Lycoris is a well known lineage of ornamental geophytes consisting of at least 20 species, some of which are thought to be of natural hybrid origin. Previous genetic studies have supported this hypothesis, but these have relied solely on the use of karyology and/or nuclear ribosomal ITS sequences. No plastid DNA data have been available to address interspecific relationships within Lycoris until now. In this study, 500 individuals from 29 populations representing 16 of the 20 published Lycoris spp. were sampled, and DNA sequences were generated for two plastid markers (trnS‐trnfM and trnC‐ycf6). From these data we inferred phylogenetic relationships among the sampled taxa at the species and population levels using concatenated phylogenetic methods. A well resolved and strongly supported phylogenetic reconstruction for Lycoris was obtained. Although the plastid DNA topology differs from that derived previously using ITS, both genomes produce trees that cluster Lycoris spp. into three clades. One of these, containing polyploid taxa such as L. albiflora, L. caldwellii, L. straminea and L. houdyshelii, shows strong evidence of reticulation, and we discuss the identity of potential parents of these allopolyploids. In contrast, we offer evidence that challenges the hypothesis that triploid individuals of L. radiata are the result of hybridization. Instead, they appear to be autotriploids that have arisen in more than one location. By comparing the phylogenetic results obtained using nuclear genomic data to those from the plastid genome, a much clearer picture of the role that hybridization and reticulation have played in the evolution of Lycoris is emerging. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 176 , 115–126.     

Allopolyploidy in the Wintergreen Group of tribe Gaultherieae (Ericaceae) inferred from low‐copy nuclear genes

DNA sequence data from the low‐copy nuclear genes waxy (GBSSI) and leafy were compared with plastid and ITS sequence data from prior studies to reconstruct phylogenetic relationships in the Wintergreen Group of tribe Gaultherieae (Ericaceae). We conducted phylogenetic analysis with 102 species that includes representatives of all 15 major clades previously identified within the Wintergreen Group and that together span its circum‐Pacific distribution. Results yielded two distinct homeologous copies of waxy for two of the clades, each in widely separated parts of the tree. It also yielded two copies of leafy for one of the clades; only one copy of leafy was found for the other clade, but it was placed in the same major clade as its waxy counterpart and well away from its placement in a prior plastid analysis. A combined four‐locus (waxy, leafy, ITS and plastid data) phylogenetic analysis of all available relevant data placed the copies of each of the clades in two distinct positions in the phylogeny with strong overall statistical support. In combination with evidence from morphology, reproductive biology and cytology, the results suggest that these clades arose through allopolyploid hybridization between lineages deep in the phylogeny but relatively close geographically. This finding confirms previous assumptions that hybridization has played an important role in the evolution of the Gaultherieae.