HYBRIDS AND PHYLOGENETIC SYSTEMATICS II. THE IMPACT OF HYBRIDS ON CLADISTIC ANALYSIS

Untying Gordian knots: unraveling reticulate polyploid plant evolution by genomic data using the large Ranunculus auricomus species complex

Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids. We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses. Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution. The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

Detection of F1 hybrids from single-genome data reveals frequent hybridization in Hymenoptera and particularly ants

Hybridization occupies a central role in many fundamental evolutionary processes, such as speciation or adaptation. Yet, despite its pivotal importance in evolution, little is known about the actual prevalence and distribution of current hybridization across the tree of life. Here we develop and implement a new statistical method enabling the detection of F1 hybrids from single-individual genome sequencing data. Using simulations and sequencing data from known hybrid systems, we first demonstrate the specificity of the method, and identify its statistical limits. Next, we showcase the method by applying it to available sequencing data from more than 1,500 species of Arthropods, including Hymenoptera, Hemiptera, Coleoptera, Diptera, and Archnida. Among these taxa, we find Hymenoptera, and especially ants, to display the highest number of candidate F1 hybrids, suggesting higher rates of recent hybridization between previously isolated gene pools in these groups. The prevalence of F1 hybrids was heterogeneously distributed across ants, with taxa including many candidates tending to harbor specific ecological and life-history traits. This work shows how large-scale genomic comparative studies of recent hybridization can be implemented, uncovering the determinants of first-generation hybridization across whole taxa.

Novel Approaches for Species Concepts and Delimitation in Polyploids and Hybrids

Hybridization and polyploidization are important processes for plant evolution. However, classification of hybrid or polyploid species has been notoriously difficult because of the complexity of processes and different evolutionary scenarios that do not fit with classical species concepts. Polyploid complexes are formed via combinations of allopolyploidy, autopolyploidy and homoploid hybridization with persisting sexual reproduction, resulting in many discrete lineages that have been classified as species. Polyploid complexes with facultative apomixis result in complicated net-work like clusters, or rarely in agamospecies. Various case studies illustrate the problems that apply to traditional species concepts to hybrids and polyploids. Conceptual progress can be made if lineage formation is accepted as an inevitable consequence of meiotic sex, which is established already in the first eukaryotes as a DNA restoration tool. The turnaround of the viewpoint that sex forms species as lineages helps to overcome traditional thinking of species as "units". Lineage formation and self-sustainability is the prerequisite for speciation and can also be applied to hybrids and polyploids. Species delimitation is aided by the improved recognition of lineages via various novel -omics methods, by understanding meiosis functions, and by recognizing functional phenotypes by considering morphological-physiological-ecological adaptations.

HybPhaser: A workflow for the detection and phasing of hybrids in target capture data sets

PREMISE

Hybrids contain divergent alleles that can confound phylogenetic analyses but can provide insights into reticulated evolution when identified and phased. We developed a workflow to detect hybrids in target capture data sets and phase reads into parental lineages using a similarity and phylogenetic framework.

METHODS

We used Angiosperms353 target capture data for Nepenthes, including known hybrids to test the novel workflow. Reference mapping was used to assess heterozygous sites across the data set and to detect hybrid accessions and paralogous genes. Hybrid samples were phased by mapping reads to multiple references and sorting reads according to similarity. Phased accessions were included in the phylogenetic framework.

RESULTS

All known Nepenthes hybrids and nine additional samples had high levels of heterozygous sites, had reads associated with multiple divergent clades, and were phased into accessions resembling divergent haplotypes. Phylogenetic analysis including phased accessions increased clade support and confirmed parental lineages of hybrids.

DISCUSSION

HybPhaser provides a novel approach to detect and phase hybrids in target capture data sets, which can provide insights into reticulations by revealing origins of hybrids and reduce conflicting signal, leading to more robust phylogenetic analyses.

Unravelling hybridization in Phytophthora using phylogenomics and genome size estimation

The genus Phytophthora comprises many economically and ecologically important plant pathogens. Hybrid species have previously been identified in at least six of the 12 phylogenetic clades. These hybrids can potentially infect a wider host range and display enhanced vigour compared to their progenitors. Phytophthora hybrids therefore pose a serious threat to agriculture as well as to natural ecosystems. Early and correct identification of hybrids is therefore essential for adequate plant protection but this is hampered by the limitations of morphological and traditional molecular methods. Identification of hybrids is also important in evolutionary studies as the positioning of hybrids in a phylogenetic tree can lead to suboptimal topologies. To improve the identification of hybrids we have combined genotyping-by-sequencing (GBS) and genome size estimation on a genus-wide collection of 614 Phytophthora isolates. Analyses based on locus- and allele counts and especially on the combination of species-specific loci and genome size estimations allowed us to confirm and characterize 27 previously described hybrid species and discover 16 new hybrid species. Our method was also valuable for species identification at an unprecedented resolution and further allowed correct naming of misidentified isolates. We used both a concatenation- and a coalescent-based phylogenomic method to construct a reliable phylogeny using the GBS data of 140 non-hybrid Phytophthora isolates. Hybrid species were subsequently connected to their progenitors in this phylogenetic tree. In this study we demonstrate the application of two validated techniques (GBS and flow cytometry) for relatively low cost but high resolution identification of hybrids and their phylogenetic relations.

Genome-wide genotyping of a novel Mexican Chile Pepper collection illuminates the history of landrace differentiation after Capsicum annuum L. domestication

Studies of genetic diversity among phenotypically distinct crop landraces improve our understanding of fruit evolution and genome structure under domestication. Chile peppers (Capsicum spp. L.) are economically valuable and culturally important species, and extensive phenotypic variation among landraces exists in southern Mexico, a center of C. annuum diversity. We collected 103 chile pepper seed accessions from 22 named landraces across 27 locations in southern Mexico. We genotyped these accessions with genotyping by sequencing (GBS), yielding 32,623 filtered single-nucleotide polymorphisms. Afterward, we genotyped 32 additional C. annuum accessions from a global collection for comparison to the Mexican collection. Within the Mexican collection, genetic assignment analyses showed clear genetic differentiation between landraces and clarified the unique nature of the Tusta landrace. Further clustering analyses indicated that the largest fresh-use Chile de Agua and dry-use Costeño landraces were part of separate clades, indicating that these two landraces likely represent distinct populations. The global accessions showed considerable admixture and limited clustering, which may be due to the collapse of use-type divisions outside of Central America. The separation of the Mexican landraces in part by fruit morphology related to use highlights the relevance of this use-type morphological diversity for plant breeders and the utility of fruit development variation for evolutionary biologists.

A genetic legacy of introgression confounds phylogeny and biogeography in oaks

Oaks (Quercus L.) have long been suspected to hybridize in nature, and widespread genetic exchange between morphologically defined species is well documented in two- to six-species systems, but the phylogenetic consequences of hybridization in oaks have never been demonstrated in a phylogenetically diverse sample. Here, we present phylogenomic analyses of a ca 30 Myr clade that strongly support morphologically defined species and the resolution of novel clades of white oaks; however, historical hybridization across clade boundaries is detectable and, undiagnosed, would obscure the imprint of biogeographic history in the phylogeny. Phylogenetic estimation from restriction-site-associated DNA sequencing data for 156 individuals representing 81 species supports two successive intercontinental disjunctions of white oaks: an early vicariance between the Eurasian and American white oaks, and a second, independent radiation represented by two relictual species. A suite of subsampled and partitioned analyses, however, supports a more recent divergence of the Eurasian white oaks from within the American white oaks and suggests that historic introgression between the Eurasian white oaks and a now-relictual lineage biases concatenated phylogenetic estimates. We demonstrate how divergence and reticulation both influence our understanding of the timing and nature of diversification and global colonization in these ecologically and economically important taxa.

A reassessment of the phylogeny and circumscription of Zaluzianskya (Scrophulariaceae)

The genus Zaluzianskya (Scrophulariaceae s.s.) encompasses a diversity of floral and ecological traits. However, this diversity, as described by the current taxonomic circumscription of Zaluzianskya, is an underestimate. We present molecular data suggesting that this genus requires expansion via incorporation of species from other genera and recognition of unnamed cryptic species. This study advances prior molecular phylogenies of the southern African genus through the addition of DNA regions and 51 populations that had not previously been sampled in a published phylogeny. A total of 82 species of Zaluzianskya and related genera are included, adding 48 to those previously sampled. Results are presented from analyses of five DNA regions, including nuclear ITS and four rapidly evolving chloroplast regions (trnL-trnF, rpl16, rps16, and trnS-trnfM). Our primary finding is that the genus Phyllopodium is polyphyletic as currently circumscribed, with some species placed within Zaluzianskya and others grouping with Polycarena, indicating the need for further phylogenetic work on these genera. Preliminary support for the incorporation of Reyemia into Zaluzianskya is reinforced here by the first molecular analysis to include both species of Reyemia and a strong sampling of species across Zaluzianskya and major clades of tribe Limoselleae. The two disjunct, tropical African species of Zaluzianskya are also confirmed as members of this genus. Finally, a broad sampling of 21 populations of Z. microsiphon establishes their phylogenetic division into two to five separate lineages. Hybridization, coevolution, and cryptic speciation may each play a role in the evolution of Z. microsiphon. Further resolution within a clade comprising sections Nycterinia and Macrocalyx is needed to better understand their relationships.

Multilocus phylogeny and phylogenomics of Eriochrysis P. Beauv. (Poaceae-Andropogoneae): Taxonomic implications and evidence of interspecific hybridization

Species delimitation is a vital issue concerning evolutionary biology and conservation of biodiversity. However, it is a challenging task for several reasons, including the low interspecies variability of markers currently used in phylogenetic reconstructions and the occurrence of reticulate evolution and polyploidy in many lineages of flowering plants. The first phylogeny of the grass genus Eriochrysis is presented here, focusing on the New World species, in order to examine its relationships to other genera of the subtribe Saccharinae/tribe Andropogoneae and to define the circumscriptions of its taxonomically complicated species. Molecular cloning and sequencing of five regions of four low-copy nuclear genes (apo1, d8, ep2-ex7 and ep2-ex8, kn1) were performed, as well as complete plastome sequencing. Trees were reconstructed using maximum parsimony, maximum likelihood, and Bayesian inference analyses. The present phylogenetic analyses indicate that Eriochrysis is monophyletic and the Old World E. pallida is sister to the New World species. Subtribe Saccharinae is polyphyletic, as is the genus Eulalia. Based on nuclear and plastome sequences plus morphology, we define the circumscriptions of the New World species of Eriochrysis: E. laxa is distinct from E. warmingiana, and E. villosa is distinct from E. cayennensis. Natural hybrids occur between E. laxa and E. villosa. The hybrids are probably tetraploids, based on the number of paralogues in the nuclear gene trees. This is the first record of a polyploid taxon in the genus Eriochrysis. Some incongruities between nuclear genes and plastome analyses were detected and are potentially caused by incomplete lineage sorting and/or ancient hybridization. The set of low-copy nuclear genes used in this study seems to be sufficient to resolve phylogenetic relationships and define the circumscriptions of other species complexes in the grass family and relatives, even in the presence of polyploidy and reticulate evolution. Complete plastome sequencing is also a promising tool for phylogenetic inference.

Phylogenetic analysis of Saccharum s.l. (Poaceae; Andropogoneae), with emphasis on the circumscription of the South American species

PREMISE OF THE STUDY

Polyploidy and reticulate evolution are often a complication for discovering phylogenetic relationships between genera and species. Despite the huge economic importance of sugarcane (Saccharum officinarum-Poaceae, Andropogoneae), the limits of the genus Saccharum and its species are complex and largely unresolved, involving both polyploidy and reticulate evolution. This study aimed to assess the phylogenetic relationships of Saccharum s.l., including Erianthus and Tripidium, as well as investigate the taxonomic circumscription of the South American species of the genus.

METHODS

Molecular cloning and sequencing of five regions of four low-copy nuclear loci were performed, including Aberrant panicle organization1 (apo1), Dwarf8 (d8), two exons of Erect panicle2 (ep2-ex7 and ep2-ex8), and Retarded palea1 (rep1). Concatenated trees were reconstructed using Maximum Parsimony, Maximum Likelihood, and Bayesian Inference analyses.

KEY RESULTS

The allopolyploid origin of Saccharum was demonstrated using evidence from nuclear genes. The samples of Saccharum s.l. grouped in two distinct clades, with S. arundinaceum and S. ravennae (= Tripidium, or Erianthus sect. Ripidium) apart from all other species analyzed of the genus. Saccharum angustifolium, S. asperum, and S. villosum correspond to distinct clades (different species). The plants with intermediate morphology between S. angustifolium and S. villosum presented a pattern of paralogues consistent with a hybrid origin.

CONCLUSIONS

Saccharum s.l. is polyphyletic and Tripidium should be recognized as a distinct genus. However, no strong evidence was found to support the segregation of Erianthus. The taxonomic circumscription of the South American species of the genus was resolved and the occurrence of natural hybrids was documented. Better understanding of the phylogenetic relationships of Saccharum and relatives may be useful for sugarcane breeders to identify potential taxa for interspecific and intergeneric crosses in the genetic improvement of sugarcane.

Nuclear ribosomal ITS functional paralogs resolve the phylogenetic relationships of a late-Miocene radiation cycad Cycas (Cycadaceae)

Cycas is the most widespread and diverse genus among the ancient cycads, but the extant species could be the product of late Miocene rapid radiations. Taxonomic treatments to date for this genus are quite controversial, which makes it difficult to elucidate its evolutionary history. We cloned 161 genomic ITS sequences from 31 species representing all sections of Cycas. The divergent ITS paralogs were examined within each species and identified as putative pseudogenes, recombinants and functional paralogs. Functional paralogs were used to reconstruct phylogenetic relationships with pseudogene sequences as molecular outgroups, since an unambiguous ITS sequence alignment with their closest relatives, the Zamiaceae, is unachievable. A fully resolved and highly supported tree topology was obtained at the section level, with two major clades including six minor clades. The results fully supported the classification scheme proposed by Hill (2004) at the section level, with the minor clades representing his six sections. The two major clades could be recognised as two subgenera. The obtained pattern of phylogenetic relationships, combined with the different seed dispersal capabilities and paleogeography, allowed us to propose a late Miocene rapid radiation of Cycas that might have been promoted by vicariant events associated with the complex topography and orogeny of South China and adjacent regions. In contrast, transoceanic dispersals might have played an important role in the rapid diversification of sect. Cycas, whose members have evolved a spongy layer in their seeds aiding water dispersals.

Promiscuous speciation with gene flow in silverside fish genus Odontesthes (Atheriniformes, Atherinopsidae) from south western Atlantic Ocean basins

The present paper integrates phylogenetic and population genetics analyses based on mitochondrial and nuclear molecular markers in silversides, genus Odontesthes, from a non-sampled area in the SW Atlantic Ocean to address species discrimination and to define Managements Units for sustainable conservation. All phylogenetic analyses based on the COI mitochondrial gene were consistent to support the monophyly of the genus Odontesthes and to include O. argentinensis, O. perugiae-humensis and some O. bonariensis haplotypes in a basal polytomy conforming a major derivative clade. Microsatellites data revealed somewhat higher genetic variability values in the O. argentinensis-perugia populations than in O. bonariensis and O. perugia-humensis taxa. Contrasting population genetics structuring emerged from mitochondrial and microsatellites analyses in these taxa. Whereas mitochondrial data supported two major groups (O. argentinensis-perugia-humensis vs. O. bonariensis-perugiae-humensis populations), microsatellite data detected three major genetic entities represented by O. bonariensis, O. perugiae-humensis and an admixture of populations belonging to O. argentinensis-perugiae respectively. Therefore, the star COI polytomy in the tree topology involving these taxa could be interpreted by several hypothetic scenarios such as the existence of shared ancestral polymorphisms, incomplete lineage sorting in a radiating speciation process and/or reticulation events. Present findings support that promiscuous and recent contact between incipient species sharing asymmetric gene flow exchanges, blurs taxa boundaries yielding complicated taxonomy and Management Units delimitation in silverside genus Odontesthes from SW Atlantic Ocean basins.

Phylogeny of a genomically diverse group of elymus (poaceae) allopolyploids reveals multiple levels of reticulation

The grass tribe Triticeae (=Hordeeae) comprises only about 300 species, but it is well known for the economically important crop plants wheat, barley, and rye. The group is also recognized as a fascinating example of evolutionary complexity, with a history shaped by numerous events of auto- and allopolyploidy and apparent introgression involving diploids and polyploids. The genus Elymus comprises a heterogeneous collection of allopolyploid genome combinations, all of which include at least one set of homoeologs, designated St, derived from Pseudoroegneria. The current analysis includes a geographically and genomically diverse collection of 21 tetraploid Elymus species, and a single hexaploid species. Diploid and polyploid relationships were estimated using four molecular data sets, including one that combines two regions of the chloroplast genome, and three from unlinked nuclear genes: phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I. Four gene trees were generated using maximum likelihood, and the phylogenetic placement of the polyploid sequences reveals extensive reticulation beyond allopolyploidy alone. The trees were interpreted with reference to numerous phenomena known to complicate allopolyploid phylogenies, and introgression was identified as a major factor in their history. The work illustrates the interpretation of complicated phylogenetic results through the sequential consideration of numerous possible explanations, and the results highlight the value of careful inspection of multiple independent molecular phylogenetic estimates, with particular focus on the differences among them.

Phylogeny of the Asparagales based on three plastid and two mitochondrial genes

PREMISE OF THE STUDY

The Asparagales, with ca. 40% of all monocotyledons, include a host of commercially important ornamentals in families such as Orchidaceae, Alliaceae, and Iridaceae, and several important crop species in genera such as Allium, Aloe, Asparagus, Crocus, and Vanilla. Though the order is well defined, the number of recognized families, their circumscription, and relationships are somewhat controversial.

METHODS

Phylogenetic analyses of Asparagales were based on parsimony and maximum likelihood using nucleotide sequence variation in three plastid genes (matK, ndhF, and rbcL) and two mitochondrial genes (atp1 and cob). Branch support was assessed using both jackknife analysis implementing strict-consensus (SC) and bootstrap analysis implementing frequency-within-replicates (FWR). The contribution of edited sites in the mitochondrial genes to topology and branch support was investigated.

KEY RESULTS

The topologies recovered largely agree with previous results, though some clades remain poorly resolved (e.g., Ruscaceae). When the edited sites were included in the analysis, the plastid and mitochondrial genes were highly incongruent. However, when the edited sites were removed, the two partitions became congruent.

CONCLUSIONS

Some deeper nodes in the Asparagales tree remain poorly resolved or unresolved as do the relationships of certain monogeneric families (e.g., Aphyllanthaceae, Ixioliriaceae, Doryanthaceae), whereas support for many families increases. However, the increased support is dominated by plastid data, and the potential influence of mitochondrial and biparentially inherited single or low-copy nuclear genes should be investigated.

Phylogenetic relationships and natural hybridization among the North American woody bamboos (Poaceae: Bambusoideae: Arundinaria)

In spite of the ecological and economic importance of temperate bamboos, relatively little is known about their population biology or evolutionary history. Recently, hybridization has emerged as a potential source of diversity in this group, as well as an underlying cause of taxonomic problems. As part of a broader phylogenetic study of the temperate bamboos, we report the results of an analysis of the North American Arundinaria gigantea species complex, including estimates of genetic variation and molecular evidence of natural hybridization among A. gigantea, A. tecta, and A. appalachiana. The study involved a comparative analysis of amplified fragment length polymorphisms (AFLPs) and chloroplast DNA sequences representing diversity within and among all three species plus individuals with intermediate or unusual morphological characteristics (putative hybrids). Molecular results support the recognition of three species previously defined on the basis of morphology, anatomy, and ecology, with most of the molecular variance accounted for by among-species variation. Molecular evidence also demonstrates that A. tecta and A. appalachiana are sister species, forming a clade that is significantly divergent from A. gigantea. The role of hybridization in the phylogenetic history of Arundinaria is discussed along with implications for the evolution and taxonomy of the temperate woody bamboos.

Species delimitation: a case study in a problematic ant taxon

Species delimitation has been invigorated as a discipline in systematics by an influx of new character sets, analytical methods, and conceptual advances. We use genetic data from 68 markers, combined with distributional, bioclimatic, and coloration information, to hypothesize boundaries of evolutionarily independent lineages (species) within the widespread and highly variable nominal fire ant species Solenopsis saevissima, a member of a species group containing invasive pests as well as species that are models for ecological and evolutionary research. Our integrated approach uses diverse methods of analysis to sequentially test whether populations meet specific operational criteria (contingent properties) for candidacy as morphologically cryptic species, including genetic clustering, monophyly, reproductive isolation, and occupation of distinctive niche space. We hypothesize that nominal S. saevissima comprises at least 4-6 previously unrecognized species, including several pairs whose parapatric distributions implicate the development of intrinsic premating or postmating barriers to gene flow. Our genetic data further suggest that regional genetic differentiation in S. saevissima has been influenced by hybridization with other nominal species occurring in sympatry or parapatry, including the quite distantly related Solenopsis geminata. The results of this study illustrate the importance of employing different classes of genetic data (coding and noncoding regions and nuclear and mitochondrial DNA [mtDNA] markers), different methods of genetic data analysis (tree-based and non-tree based methods), and different sources of data (genetic, morphological, and ecological data) to explicitly test various operational criteria for species boundaries in clades of recently diverged lineages, while warning against over reliance on any single data type (e.g., mtDNA sequence variation) when drawing inferences.

Distinguishing Terminal Monophyletic Groups from Reticulate Taxa: Performance of Phenetic, Tree-Based, and Network Procedures

Hybridization is a well-documented, natural phenomenon that is common at low taxonomic levels in the higher plants and other groups. In spite of the obvious potential for gene flow via hybridization to cause reticulation in an evolutionary tree, analytical methods based on a strictly bifurcating model of evolution have frequently been applied to data sets containing taxa known to hybridize in nature. Using simulated data, we evaluated the relative performance of phenetic, tree-based, and network approaches for distinguishing between taxa with known reticulate history and taxa that were true terminal monophyletic groups. In all methods examined, type I error (the erroneous rejection of the null hypothesis that a taxon of interest is not monophyletic) was likely during the early stages of introgressive hybridization. We used the gradual erosion of type I error with continued gene flow as a metric for assessing relative performance. Bifurcating tree-based methods performed poorly, with highly supported, incorrect topologies appearing during some phases of the simulation. Based on our model, we estimate that many thousands of gene flow events may be required in natural systems before reticulate taxa will be reliably detected using tree-based methods of phylogeny reconstruction. We conclude that the use of standard bifurcating tree-based methods to identify terminal monophyletic groups for the purposes of defining or delimiting phylogenetic species, or for prioritizing populations for conservation purposes, is difficult to justify when gene flow between sampled taxa is possible. As an alternative, we explored the use of two network methods. Minimum spanning networks performed worse than most tree-based methods and did not yield topologies that were easily interpretable as phylogenies. The performance of NeighborNet was comparable to parsimony bootstrap analysis. NeighborNet and reverse successive weighting were capable of identifying an ephemeral signature of reticulate evolution during the early stages of introgression by revealing conflicting phylogenetic signal. However, when gene flow was topologically complex, the conflicting phylogenetic signal revealed by these methods resulted in a high probability of type II error (inferring that a monophyletic taxon has a reticulate history). Lastly, we present a novel application of an existing nonparametric clustering procedure that, when used against a density landscape derived from principal coordinate data, showed superior performance to the tree-based and network procedures tested.

Amplified fragment length polymorphisms and sequence data in the phylogenetic analysis of polyploids: multiple origins of Veronica cymbalaria (Plantaginaceae)

The origin of polyploid Veronica cymbalaria (Plantaginaceae) was investigated using DNA sequence data and amplified fragment length polymorphism (AFLP) fingerprints to reveal the parentage of this taxon. The use of AFLP fingerprints in phylogenetic analysis is problematic and various methods have therefore been compared. DNA sequence data (for the internal transcribed spacer (ITS) region and the plastid trnL-F region (trnL intron, 3'exon, and trnL-F spacer)) and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis of the ITS region suggested a reliable hypothesis for the evolution of the V. cymbalaria complex. This hypothesis allowed evaluation of the effect of different distance measures (Jaccard and Nei-Li) in phenetic, character-state weighted parsimony, and Bayesian analyses of AFLP markers. The study establishes that tetraploid V. cymbalaria originated at least twice in the eastern Mediterranean, with one parent differing in the two separate origins. Hexaploid V. cymbalaria originated even more often. The results illustrate that even subtle differences in the analyses of AFLP markers can lead to drastically different conclusions. The study reveals multiple origins of a Mediterranean polyploid species. Furthermore, it demonstrates that the analysis of a complex marker system such as AFLP fingerprints using only one type of analysis can easily be misleading.

Phylogenetic relationships and possible hybrid origin of Lycoris species (Amaryllidaceae) revealed by its sequences

To examine interspecific relationships and test the hypothesis of hybrid origin within Lycoris species, this study used data from parsimony analyses with nuclear ITS sequences for 19 taxa representing 14 species of Lycoris and two outgroup taxa. The ITS sequences resolved three infrageneric clades. One clade included L. chinensis, L. longituba, L. longituba var. flava, L. anhuiensis, and L. aurea; the second one consisted of L. sprengeri, L. radiata, L. radiata var. radiata, L. radiata var. pumila, L. haywardii, L. rosea, L. sanguinea var. sanguinea, and L. sanguinea var. koreana; and the third included L. caldwellii, L. straminea, L. albiflora, L. flavescens, and two hybrids. The results strongly support the hypothesis that L. straminea originated from hybridization between L. chinensis and L. radiata var. pumila, and the allotriploid L. caldwellii and L. albiflora derived from hybridization between L. chinensis and L. sprengeri. As nucleotide additivity was observed in the artificial hybrids and several presumed hybrids, the likelihood of hybrid origin of Lycoris species is supported.

Genetic differentiation of geographically separated populations of the southern green stink bug Nezara viridula (Hemiptera: Pentatomidae)

Genetic variation in the southern green stink bug Nezara viridula (Linnaeus) from 11 geographically separated sampling locations (Slovenia, France, Greece, Italy, Madeira, Japan, Guadeloupe, Galapagos, California, Brazil and Botswana) was studied by sequencing 16S and 28S rDNA, cytochrome b and cytochrome c oxidase subunit I gene fragments and random amplified polymorphic DNA (RAPD) analysis. Sequencing revealed 11 distinct haplotypes clustering into lineages A, B and C. Lineage C was characteristic for a single analysed specimen from Botswana. Lineage B was detected in Japan, and it probably arose in Asia. Haplotypes of European and American specimens belonged to lineage A; specimens from France, Slovenia, Madeira and Brazil shared highly similar haplotypes (>99%) from subgroup A1, while all the specimens from Greece, California, Galapagos and Guadeloupe shared a haplotype from subgroup A2. RAPD data were more variable but consistent with mtDNA sequences, revealing the same clustering. They separated the Botswanian specimen from Japanese specimens and from a group of more closely related specimens from Europe and America. Sequence and RAPD results both support the African origin of N. viridula, followed by dispersal to Asia (lineage B) and, more recently, by expansion to Europe and America (lineage A). RAPD analysis revealed two highly supported subgroups in Japan, congruent with mtDNA lineages A2 and B, suggesting multiple colonization of Japan. Invariant sequences at the 28S rDNA combined with other results do not support the hypothesis that cryptic (sibling) species exist within the populations investigated in this study.

Testing Introgressive Hybridization Hypotheses Using Statistical Network Analysis of Nuclear and Cytoplasmic Haplotypes in the Leaf Beetle Timarcha goettingensis Species Complex

Previous studies of leaf beetles (Chrysomelidae) in the Timarcha goettingensis species complex using mitochondrial (cox2) and nuclear (ITS-2 rRNA) markers revealed two main clades confined to the Iberian Peninsula and the rest of Europe but showing incongruent distributions indicative of gene exchange between both groups. Because of the anastomosing nature of hybridization, which disrupts the cladistic structure of character variation, phylogenetic trees might be inappropriate to represent and study this process. Here we test for evidence of hybridization in the T. goettingensis complex by analyzing the extra homoplasy arising in hybrid genomes from the simultaneous analysis of genetically independent markers. Haplotype networks obtained by Templeton's statistical parsimony analysis were generated for combined (concatenated) cox2 and ITS-2 sequences from 167 individuals of the T. goettingensis complex. Networks were used to detect runs of homoplasious characters physically clustered along a nucleotide sequence, as evidence for recombination between both gene partitions. A hypergeometric tail probability for the chance occurrence of physically clustered character changes on the connections linking networks of genotypes was applied. The test recognized two instances of statistically significant clustering, indicating the presence of cox2-ITS-2 mosaic genotypes and reticulation of both main T. goettingensis clades, supporting the reticulate origin of samples of T. maritima in southwestern France and T. sinuatocollis/T. monserratensis in the eastern Pyrenees. Although the assessment of reticulation in DNA sequences does not provide direct proof for hybridization, the geographical distribution of mosaic genotypes in the vicinity of "pure" genotypes supports the effect of gene flow between the two divergent lineages. The study demonstrates the utility of statistical parsimony networks for the detection of hybrids in the growing number of phylogeographic studies based on multiple gene markers.

Hybridization in the recent past

The question we address in this article is how hybridization in the recent past can be detected in recently evolved species. Such species may not have evolved genetic incompatibilities and may hybridize with little or no fitness loss. Hybridization can be recognized by relatively small genetic differences between sympatric populations because sympatric populations have the opportunity to interbreed whereas allopatric populations do not. Using microsatellite DNA data from Darwin's finches in the Galapagos archipelago, we compare sympatric and allopatric genetic distances in pairs of Geospiza and Camarhynchus species. In agreement with the hybridization hypothesis, we found a statistically strong tendency for a species to be more similar genetically to a sympatric relative than to allopatric populations of that relative. Hybridization has been studied directly on two islands, but it is evidently more widespread in the archipelago. We argue that introgressive hybridization may have been a persistent feature of the adaptive radiation through most of its history, facilitating evolutionary diversification and occasionally affecting both the speed and direction of evolution.

Reconstructing patterns of reticulate evolution in plants

Until recently, rigorously reconstructing the many hybrid speciation events in plants has not been practical because of the limited number of molecular markers available for plant phylogenetic reconstruction and the lack of good, biologically based methods for inferring reticulation (network) events. This situation should change rapidly with the development of multiple nuclear markers for phylogenetic reconstruction and new methods for reconstructing reticulate evolution. These developments will necessitate a much greater incorporation of population genetics into phylogenetic reconstruction than has been common. Population genetic events such as gene duplication coupled with lineage sorting and meiotic and sexual recombination have always had the potential to affect phylogenetic inference. For tree reconstruction, these problems are usually minimized by using uniparental markers and nuclear markers that undergo rapid concerted evolution. Because reconstruction of reticulate speciation events will require nuclear markers that lack these characteristics, effects of population genetics on phylogenetic inference will need to be addressed directly. Current models and methods that allow hybrid speciation to be detected and reconstructed are discussed, with a focus on how lineage sorting and meiotic and sexual recombination affect network reconstruction. Approaches that would allow inference of phylogenetic networks in their presence are suggested.

Hybridization and gene flow between a day- and night-flowering species of Zaluzianskya (Scrophulariaceae s.s., tribe Manuleeae)

Despite apparent ethological isolation based on specialized pollination systems, hybridization between day-flowering Zaluzianskya microsiphon and night-flowering Z. natalensis has been proposed due to intermediate individuals found in sympatric populations of these species. The extent of this putative hybridization was investigated using inter-simple sequence repeat (ISSR) markers and principal components analysis (PCA) of morphological traits. The species are genetically similar, but show some intra- and interspecific variation in band frequencies. Neighbor-joining analyses of the ISSR data demonstrated that although the species largely formed distinct groups, several individuals from the sympatric populations of each species and the "hybrids" clustered together rather than with members of their own species. These results are consistent with hybridization, although they could also indicate historical similarity. Nine loci were present only in individuals of Z. microsiphon, the "hybrids," and sometimes the sympatric individuals of Z. natalensis. In contrast, only one locus showed the reverse pattern. This suggests unidirectional gene flow from Z. microsiphon to Z. natalensis, which is also supported by population-level examinations of four loci. Ordination revealed separate phenotype clusters for each species, with hybrid individuals located in between but often closer to the Z. natalensis cluster. One hypothesis is that hybrids are backcrossing with Z. natalensis, leading to introgression of Z. microsiphon genetic material.

Radiation of the Australian Salicornioideae (Chenopodiaceae)--based on evidence from nuclear and chloroplast DNA sequences

In phylogenetic analyses of nuclear ITS and chloroplast trnL DNA sequences, the mostly endemic Australian genera; Halosarcia, Pachycornia, Sclerostegia, Tecticornia, and Tegicornia of the subfamily Salicornioideae (Chenopodiaceae) together form a monophyletic group, congruent with the hypothesis that they evolved from a common ancestor. However, limited genetic differentiation evident in both nrDNA and cpDNA sequences among these taxa suggests a possible rapid radiation. Based on fossil pollen records and climatic models of other authors, it is hypothesized that the expansion of the Australian endemic Salicornioideae most likely occurred during the Late Miocene to Pliocene, when increasing aridity caused the formation of extensive salt lakes along endorheic paleodrainage channels. Moreover, Australian Sarcocornia representatives were supported as monophyletic, nested within a paraphyletic Sarcocornia clade that also comprised European Salicornia in the ITS analysis. This suggests that Sarcocornia arrived in Australia subsequent to the ancestor of the Australian endemic genera most likely via long-distance dispersal.

Molecular and cytological examination of Calopogon (Orchidaceae, Epidendroideae): circumscription, phylogeny, polyploidy, and possible hybrid speciation

The orchid genus Calopogon R.Br. (Orchidaceae), native to eastern North America and the northern Caribbean, currently contains five species and up to three varieties. Using nuclear internal transcribed spacer (ITS) ribosomal DNA sequences, amplified fragment length polymorphisms (AFLPs), chloroplast DNA restriction fragments, and chromosome counts, we present a phylogenetic and taxonomic study of the genus. Calopogon multiflorus and C. pallidus are consistently sister species, but the relationships of C. barbatus, C. oklahomensis, and C. tuberosus are not as clear. In the ITS analysis C. oklahomensis is sister to C. barbatus, whereas it is sister to C. tuberosus in the plastid restriction fragment analysis. Furthermore, all species were found to have chromosome numbers of 2n = 38 and 40, with the exception of the putatively hybrid-derived C. oklahomensis with 2n = 114 and 120. The hexaploidy of the latter, plus the discrepancy in its position between the ITS and plastid restriction fragment trees, could suggest that it is of hybrid origin. However, the presence of unique morphological and molecular characters might indicate that it is either an ancient hybrid or not of hybrid derivation at all. Finally, using these molecular methods all taxa appear to generally be discrete groups, with the exception of C. tuberosus vars. latifolius and tuberosus, the former of which is best combined with the latter.

Molecular evidence for bicontinental hybridogenous genomic constitution in Lepidium sensu stricto (Brassicaceae) species from Australia and New Zealand

Lepidium sensu stricto (s.s.) (Brassicaceae) (ca. 150 species) is distributed worldwide with endemic species on every continent. It is represented in Australia and New Zealand by 19 and seven native species, respectively. In the present study we used a nuclear ribosomal internal transcribed spacer (ITS) phylogeny in comparison with a cpDNA phylogeny to unravel the origin of Australian/New Zealand species. Although phylogenetic relationships within Lepidium s.s. were not fully resolved, the cpDNA data were in agreement with a Californian origin of Lepidium species from Australia/New Zealand. Strongly conflicting signals between the cp- and nuclear DNA phylogenetic analysis clearly indicated hybridogenous genomic constitution of Australian Lepidium s.s. species: All 18 studied Australian/New Zealand Lepidium s.s. species examined shared a Californian cpDNA type. While eleven Australian/New Zealand species appeared to harbor a Californian ITS type, a group of seven species shared a South African ITS type. This pattern is most likely explained by two trans-oceanic dispersals of Lepidium from California and Africa to Australia/New Zealand and subsequent hybridization followed by homogenization of the ribosomal DNA either to the Californian or South African ITS type in the two different lineages. Calibration of our molecular trees indicates a Pliocene/Pleistocene origin of Lepidium in Australia/New Zealand. Low levels of cpDNA and ITS sequence divergence and unresolved topologies within Australian/New Zealand species suggest a rapid and recent radiation of Lepidium after the hybridization event. This coincides with dramatic climatic changes in that geological epoch shaping the composition of the vegetation.

The phylogeny of Gaura (Onagraceae) based on ITS, ETS, and trnL-F sequence data

Gaura (Onagraceae: Onagreae) is a small North American genus of 21 species consisting mostly of night-blooming, moth-pollinated annuals and perennials. The current infrageneric classification based on differences in habit, floral symmetry, and fruit morphology recognizes eight sections within the genus. We examine the phylogenetic relationships of all 21 species of Gaura using DNA sequence data from the internal transcribed spacer region (ITS), the external transcribed spacer region (ETS), and the plastid trnL-F region. Combined analysis of these regions indicate Gaura is monophyletic only if it includes Stenosiphon, a monotypic genus comprised of S. linifolius. Within Gaura, our studies indicate that sections Gauridium, Schizocarya, Campogaura, Stipogaura, Xenogaura, and Gaura are monophyletic, but sections Xerogaura and Pterogaura are not and should be reevaluated. In addition, molecular data provide support for the hypothesis that G. sinuata and G. drummondii arose via interspecific hybridization followed by genome doubling; their influence on phylogenetic reconstruction is discussed.

Ribosomal ITS sequences and plant phylogenetic inference

One of the most popular sequences for phylogenetic inference at the generic and infrageneric levels in plants is the internal transcribed spacer (ITS) region of the 18S-5.8S-26S nuclear ribosomal cistron. The prominence of this source of nuclear DNA sequence data is underscored by a survey of phylogenetic publications involving comparisons at the genus level or below, which reveals that of 244 papers published over the last five years, 66% included ITS sequence data. Perhaps even more striking is the fact that 34% of all published phylogenetic hypothesis have been based exclusively on ITS sequences. Notwithstanding the many important contributions of ITS sequence data to phylogenetic understanding and knowledge of genome relationships, a number of molecular genetic processes impact ITS sequences in ways that may mislead phylogenetic inference. These molecular genetic processes are reviewed here, drawing attention to both underlying mechanism and phylogenetic implications. Among the most prevalent complications for phylogenetic inference is the existence in many plant genomes of extensive sequence variation, arising from ancient or recent array duplication events, genomic harboring of pseudogenes in various states of decay, and/or incomplete intra- or inter-array homogenization. These phenomena separately and collectively create a network of paralogous sequence relationships potentially confounding accurate phylogenetic reconstruction. Homoplasy is shown to be higher in ITS than in other DNA sequence data sets, most likely because of orthology/paralogy conflation, compensatory base changes, problems in alignment due to indel accumulation, sequencing errors, or some combination of these phenomena. Despite the near-universal usage of ITS sequence data in plant phylogenetic studies, its complex and unpredictable evolutionary behavior reduce its utility for phylogenetic analysis. It is suggested that more robust insights are likely to emerge from the use of single-copy or low-copy nuclear genes.

Origin(s) of the diploid hybrid species Helianthus deserticola (Asteraceae)

Homoploid hybrid speciation has traditionally been considered a rare event, dependent on the establishment of both a novel, balanced genotype and reproductive isolating barriers between the new species and its progenitors. However, more recent studies have shown that synthetic hybrids converge toward the chromosomal structure of natural hybrids after only a few generations, suggesting that this phenomenon may be more frequent than previously assumed. Here, the possibility that the diploid hybrid species Helianthus deserticola arose from more than one hybrid speciation event was investigated using patterns of variation from cpDNA, 18 nuclear microsatellite loci, and population interfertility. Helianthus deserticola contains cpDNA haplotypes characteristic of both parental species, is polyphyletic with one parental species based on nine microsatellite loci, and has a high degree of interfertility among populations. The data are consistent with either a single origin followed by introgression with the parental species or multiple origins. Analysis of microsatellite variation places the origin of H. deserticola between 170 000 and 63 000 years before present, making it unlikely that anthropogenic disturbances influenced its origin. Finally, the hybrid species generally has lower levels of genetic diversity but higher levels of differentiation among populations than either parental species.

The origin and evolution of Eragrostis tef (Poaceae) and related polyploids: evidence from nuclear waxy and plastid rps16

Tef (Eragrostis tef; Poaceae) is an allotetraploid (2n = 4x = 40) cereal crop whose origin within the large genus Eragrostis is unknown. Previous studies have suggested a total of 14 wild Eragrostis species as potential progenitors. Phylogenetic analysis of sequence data from the nuclear gene waxy and the plastid locus rps16 strongly supports the widely held hypothesis of a close relationship between tef and E. pilosa, a wild allotetraploid. Eragrostis heteromera, another previously proposed progenitor, is shown by the waxy data to be a close relative of one of the tef genomes. Other putative progenitors included in the taxon sample are not supported as closely related to tef. Plastid sequences from five varieties of tef and four E. pilosa accessions are identical and therefore are uninformative with respect to the question of multiple origins of these polyploids. The waxy phylogeny also resolves the relationships among other allopolyploids, supporting a close relationship between the morphologically similar allotetraploids E. macilenta, E. minor, and E. mexicana. Eragrostis cilianensis, another morphologically similar allopolyploid, appears to have shared one diploid progenitor with these species but derived its other genome from an unrelated diploid.

Intraspecific gene genealogies: trees grafting into networks

Intraspecific gene evolution cannot always be represented by a bifurcating tree. Rather, population genealogies are often multifurcated, descendant genes coexist with persistent ancestors and recombination events produce reticulate relationships. Whereas traditional phylogenetic methods assume bifurcating trees, several networking approaches have recently been developed to estimate intraspecific genealogies that take into account these population-level phenomena.

Gene genealogies reveal global phylogeographic structure and reproductive isolation among lineages of Fusarium graminearum, the fungus causing wheat scab

During the past decade, the plant disease called scab or Fusarium head blight of wheat and barley has reached epidemic proportions in North America and elsewhere in the world. Scab is an economically devastating plant disease, not only because it causes significant reduction in seed yields and quality, but also because infested seeds are often contaminated with trichothecene and estrogenic mycotoxins that pose a serious threat to animal health and food safety. To test whether the primary etiological agent of scab, the fungus Fusarium graminearum, is panmictic throughout its range, allelic genealogies were constructed from six single-copy nuclear genes from strains selected to represent the global genetic diversity of this pathogen. Excluding one hybrid strain, all six genealogies recovered the same seven biogeographically structured lineages, suggesting that they represent phylogenetically distinct species among which gene flow has been very limited during their evolutionary history. Parsimony analysis of the combined data set comprising 7,120 aligned nucleotide characters resolved most relationships among the seven lineages of the F. graminearum clade and related fusaria included in the study. Phylogenetic evidence is also presented for introgressive hybridization and intragenic recombination among lineages of the F. graminearum clade in nature.

Nuclear ribosomal DNA (nrDNA) concerted evolution in natural and artificial hybrids of Armeria (Plumbaginaceae)

Nuclear ribosomal DNA (nrDNA) internal transcribed spacer (ITS) sequences from artificial hybrids and backcrosses between Armeria villosa ssp. longiaristata and A. colorata were studied to assess the possible effects of concerted evolution in natural hybrids. F1 artificial hybrids show the expected pattern of additive polymorphisms for five of the six variable sites as estimated from direct sequences. However, homogenization of polymorphism is already observed in the F2, and is biased towards A. colorata except for one site. In backcrosses, an expected tendency towards homogenization of polymorphic sites in the direction of the recurrent parent is observed for five sites, although this does not necessarily imply concerted evolution. Conversely, the sixth site appears to elude such a mechanism and thus provides additional support for the occurrence of biased concerted evolution. Our findings are relevant to interpreting phylogeographic patterns involving gene flow and are also consistent with the hypothesis of a hybrid origin of A. villosa ssp. carratracensis.

Biogeographic patterns in Mediterranean and Macaronesian species of Saxifraga (Saxifragaceae) inferred from phylogenetic analyses of ITS sequences

A biogeographic study of Saxifraga section Saxifraga was performed based on phylogenetic analyses of ITS (internal transcribed spacer) sequences of nuclear ribosomal DNA. ITS sequences from 21 species and 31 populations were examined to identify colonization patterns for the two species of Saxifraga occurring in Macaronesia and for S. globulifera in the west Mediterranean basin. Phylogenetic analysis of the sequence data yield a single most parsimonious tree with many of the major clades well supported by bootstrap and decay values. The ITS tree provided resolution at specific and populational levels that points to two biogeographic patterns within the genus. In contrast to the molecular evidence provided by other authors for a Mediterranean origin of several Macaronesian genera of angiosperms, our results indicate that the Madeiran archipelago was colonized a single time by a species of Saxifraga originating from the Eurosiberian region. On the other hand, the molecular evidence also suggests that populations of S. globulifera from North Africa have been isolated for a long time from populations occurring in the Iberian Peninsula, and that the endemic S. reuteriana has evolved from the Iberian populations of S. globulifera. The Mediterranean Sea has probably been an effective isolating barrier for some plant groups that occur in Europe and North Africa.

Phylogeny of Rubus (rosaceae) based on nuclear ribosomal DNA internal transcribed spacer region sequences

We used nuclear ribosomal DNA internal transcribed spacer region (ITS 1 - 5.8S - ITS 2; ITS) sequences to generate the first phylogeny of Rubus based on a large, molecular data set. We sampled 57 taxa including 20 species of subgenus Rubus (blackberries), one to seven species from each of the remaining 11 subgenera, and the monotypic and closely related Dalibarda. In Rubus, ITS sequences are most informative among subgenera, and variability is low between closely related species. Parsimony analysis indicates that Rubus plus Dalibarda form a strongly supported clade, and D. repens may nest within Rubus. Of the subgenera with more than one species sampled, only subgenus Orobatus appears monophyletic. Three large clades are strongly supported: one contains all sampled species of nine of the 12 subgenera; another includes extreme Southern Hemisphere species of subgenera Comaropsis, Dalibarda, and Lampobatus; and a third clade consists of subgenus Rubus plus R. alpinus of subgenus Lampobatus. Rubus ursinus appears to be a hybrid between a close relative of R. macraei (subgenus Idaeobatus, raspberries) and an unidentified subgenus Rubus species. ITS sequences are generally consistent with biogeography and ploidy, but traditionally important morphological characters, such as stem armature and leaf type, appear to have limited phylogenetic value in Rubus.

Multiple independent losses of the plastid rpoC1 intron in Medicago (Fabaceae) as inferred from phylogenetic analyses of nuclear ribosomal DNA internal transcribed spacer sequences

A previous polymerase chain reaction based survey for the occurrence of the intron in chloroplast gene rpoC1 revealed its absence in one of the eight species of Medicago (Fabaceae; Trifolieae) examined. We extend the survey of Medicago to include 65 of the 86 species, representing all 12 sections and seven of the eight subsections recognized in the most recent comprehensive treatment of the genus. Our results indicate that 17 species from five sections lack the intron and that three of these sections are heterogeneous with regard to intron content. DNA sequencing across the rpoC1 intron-exon boundary in three of these species reveals the precise excision of the intron from the gene. Phylogenies derived from nuclear ribosomal DNA internal transcribed spacer sequences, estimated using maximum parsimony and maximum likelihood methods, suggest that most of the previously recognized sections in Medicago are not monophyletic as currently circumscribed. Furthermore, these results suggest that the rpoC1 intron has been lost independently a minimum of three times during the evolution of the group. The occurrence of multiple independent intron losses severely reduces the utility of this character as a phylogenetic marker in Medicago.Key words: Medicago, Fabaceae, chloroplast DNA, rpoC1 intron, nuclear ribosomal DNA internal transcribed spacer (ITS) sequences.

Molecular phylogenetics of Thlaspi s.l. (Brassicaceae) based on chloroplast DNA restriction site variation and sequences of the internal transcribed spacers of nuclear ribosomal DNA

Systematics of the genus Thlaspi s.l. is difficult and controversial. Previous hypotheses have been based on morphological and anatomical data. We have analyzed sequence variation of the internal transcribed spacer (ITS) regions of nuclear ribosomal DNA (nrDNA) among 13 Thlaspi s.l. taxa, representing all sections of the genus. Phylogenetic relationships among ITS sequences of the Thlaspi s.l. taxa studied are in general concordance with a previously published chloroplast DNA based phylogeny of this group. Most-parsimonious trees from ITS and chloroplast DNA data support three groups that are congruent with lineages (Thlaspi s.str., Noccaea–Raparia, Microthlaspi) previously described by Meyer on the basis of seed anatomy. The ITS data grouped Microthlaspi granatense outside the Microthlaspi clade and, therefore, Microthlaspi appeared paraphyletic on the ITS tree, in contrast with the chloroplast DNA phylogeny. We speculate that concerted evolutionary forces have acted among different nrDNA arrays (brought together in M. granatense by hybridization with a related taxon), resulting in the fixation of the alien species nrDNA type in M. granatense, which, however, maintains a Microthlaspi chloroplast genome type. Both molecular data sets detected intraspecific variation among Microthlaspi perfoliatum accessions of different geographic origin and different ploidy levels. Our molecular evidence would suggest the hybrid origin of polyploid M. perfoliatum from diploid M. perfoliatum and M. natolicum. Key words: chloroplast DNA, restriction site variation, sequence analysis of nuclear ribosomal DNA ITS regions, Thlaspi s.l. (Brassicaceae), molecular phylogeny, congruence.

Polymorphism and concerted evolution in a tandemly repeated gene family: 5S ribosomal DNA in diploid and allopolyploid cottons

5S RNA genes and their nontranscribed spacers are tandemly repeated in plant genomes at one or more chromosomal loci. To facilitate an understanding of the forces that govern 5S rDNA evolution, copy-number estimation and DNA sequencing were conducted for a phylogenetically well-characterized set of 16 diploid species of cotton (Gossypium) and 4 species representing allopolyploid derivatives of the diploids. Copy number varies over twentyfold in the genus, from approximately 1,000 to 20,000 copies/2C genome. When superimposed on the organismal phylogeny, these data reveal examples of both array expansion and contraction. Across species, a mean of 12% of nucleotide positions are polymorphic within individual arrays, for both gene and spacer sequences. This shows, in conjunction with phylogenetic evidence for ancestral polymorphisms that survive speciation events, that intralocus concerted evolutionary forces are relatively weak and that the rate of interrepeat homogenization is approximately equal to the rate of speciation. Evidence presented also shows that duplicated 5S rDNA arrays in allopolyploids have retained their subgenomic identity since polyploid formation, thereby indicating that interlocus concerted evolution has not been an important factor in the evolution of these arrays. A descriptive model, one which incorporates the opposing forces of mutation and homogenization within a selective framework, is outlined to account for the empirical data presented. Weak homogenizing forces allow equivalent levels of sequence polymorphism to accumulate in the 5S gene and spacer sequences, but fixation of mutations is nearly prohibited in the 5S gene. As a consequence, fixed interspecific differences are statistically underrepresented for 5S genes. This result explains the apparent paradox that despite similar levels of gene and spacer diversity, phylogenetic analysis of spacer sequences yields highly resolved trees, whereas analyses based on 5S gene sequences do not.