Intra-individual polymorphism in diploid and apomictic polyploid hawkweeds (Hieracium, Lactuceae, Asteraceae): disentangling phylogenetic signal, reticulation, and noise

Intragenomic rDNA variation - the product of concerted evolution, mutation, or something in between?

The classical model of concerted evolution states that hundreds to thousands of ribosomal DNA (rDNA) units undergo homogenization, making the multiple copies of the individual units more uniform across the genome than would be expected given mutation frequencies and gene redundancy. While the universality of this over 50-year-old model has been confirmed in a range of organisms, advanced high throughput sequencing techniques have also revealed that rDNA homogenization in many organisms is partial and, in rare cases, even apparently failing. The potential underpinning processes leading to unexpected intragenomic variation have been discussed in a number of studies, but a comprehensive understanding remains to be determined. In this work, we summarize information on variation or polymorphisms in rDNAs across a wide range of taxa amongst animals, fungi, plants, and protists. We discuss the definition and description of concerted evolution and describe whether incomplete concerted evolution of rDNAs predominantly affects coding or non-coding regions of rDNA units and if it leads to the formation of pseudogenes or not. We also discuss the factors contributing to rDNA variation, such as interspecific hybridization, meiotic cycles, rDNA expression status, genome size, and the activity of effector genes involved in genetic recombination, epigenetic modifications, and DNA editing. Finally, we argue that a combination of approaches is needed to target genetic and epigenetic phenomena influencing incomplete concerted evolution, to give a comprehensive understanding of the evolution and functional consequences of intragenomic variation in rDNA.

Evolutionary history and patterns of geographical variation, fertility, and hybridization in Stuckenia (Potamogetonaceae)

Aquatic plant species are often widespread, even across continents. They pose a challenge to species delimitation and taxonomy due to their reduced morphology and high phenotypic plasticity. These difficulties are even more pronounced in the case of interspecific hybridization. We investigate the aquatic plant genus Stuckenia for the first time on a worldwide scale. Expert species determination is aided by sequencing of nuclear ribosomal ITS and 5S-NTS regions and the plastid intergenic spacers rpl20-5'rps12 and trnT-trnL. Nuclear markers are used to infer hybridization, and the maternal origin of hybrids is addressed with plastid markers. Pure species are subjected to phylogenetic analyses. Two main Stuckenia lineages are found: one consists of S. amblyphylla, S. filiformis, S. pamirica, and S. vaginata, the other includes S. pectinata and S. striata. The widespread species S. pectinata, S. filiformis, and S. vaginata show intraspecific genetic variation, which is structured geographically. Many intraspecific hybrids, which are usually fertile, occur between those genotypes. Interspecific hybrids, which are consistently sterile, are detected among all widespread species; some are reported for the first time in several countries and regions. They originated multiple times from reciprocal crosses and reflect the geographical origins of parental genotypes. Intraspecific genetic variation can be higher than interspecific differences between closely related species. Comparison of phenotypic variation in the field and in cultivation with genotypic variation shows that numerous conspicuous forms have been overestimated taxonomically. These are resolved as phenotypes responding to unusual environments, have recurrently evolved adaptations, or represent extreme forms of continuous variation of the recognized species. However, some specific regional lineages, which have evolved from variable species, may be interpreted as early steps of the speciation process. Hybridization has been underestimated in some regions as a source of Stuckenia diversity, and the respective hybrid plants have been misidentified as intraspecific taxa or even as separate species. Many erroneous entries in sequence databases are detected and summarized. This work provides a sound basis for species delimitation and hybrid recognition in this difficult genus.

A Multigene Phylogeny of Native American Hawkweeds (Hieracium Subgen. Chionoracium, Cichorieae, Asteraceae): Origin, Speciation Patterns, and Migration Routes

Native American hawkweeds are mainly mountainous species that are distributed all over the New World. They are severely understudied with respect to their origin, colonization of the vast distribution area, and species relationships. Here, we attempt to reconstruct the evolutionary history of the group by applying seven molecular markers (plastid, nuclear ribosomal and low-copy genes). Phylogenetic analyses revealed that Chionoracium is a subgenus of the mainly Eurasian genus Hieracium, which originated from eastern European hawkweeds about 1.58-2.24 million years ago. Plastid DNA suggested a single origin of all Chionoracium species. They colonized the New World via Beringia and formed several distinct lineages in North America. Via one Central American lineage, the group colonized South America and radiated into more than a hundred species within about 0.8 million years, long after the closure of the Isthmus of Panama and the most recent uplift of the Andes. Despite some incongruences shown by different markers, most of them revealed the same crown groups of closely related taxa, which were, however, largely in conflict with traditional sectional classifications. We provide a basic framework for further elucidation of speciation patterns. A thorough taxonomic revision of Hieracium subgen. Chionoracium is recommended.

Horizontally Acquired nrDNAs Persist in Low Amounts in Host Hordeum Genomes and Evolve Independently of Native nrDNA

Nuclear ribosomal DNA (nrDNA) has displayed extraordinary dynamics during the evolution of plant species. However, the patterns and evolutionary significance of nrDNA array expansion or contraction are still relatively unknown. Moreover, only little is known of the fate of minority nrDNA copies acquired between species via horizontal transfer. The barley genus Hordeum (Poaceae) represents a good model for such a study, as species of section Stenostachys acquired nrDNA via horizontal transfer from at least five different panicoid genera, causing long-term co-existence of native (Hordeum-like) and non-native (panicoid) nrDNAs. Using quantitative PCR, we investigated copy number variation (CNV) of nrDNA in the diploid representatives of the genus Hordeum. We estimated the copy number of the foreign, as well as of the native ITS types (ribotypes), and followed the pattern of their CNV in relation to the genus' phylogeny, species' genomes size and the number of nrDNA loci. For the native ribotype, we encountered an almost 19-fold variation in the mean copy number among the taxa analysed, ranging from 1689 copies (per 2C content) in H. patagonicum subsp. mustersii to 31342 copies in H. murinum subsp. glaucum. The copy numbers did not correlate with any of the genus' phylogeny, the species' genome size or the number of nrDNA loci. The CNV was high within the recognised groups (up to 13.2 × in the American I-genome species) as well as between accessions of the same species (up to 4×). Foreign ribotypes represent only a small fraction of the total number of nrDNA copies. Their copy numbers ranged from single units to tens and rarely hundreds of copies. They amounted, on average, to between 0.1% (Setaria ribotype) and 1.9% (Euclasta ribotype) of total nrDNA. None of the foreign ribotypes showed significant differences with respect to phylogenetic groups recognised within the sect. Stenostachys. Overall, no correlation was found between copy numbers of native and foreign nrDNAs suggesting the sequestration and independent evolution of native and non-native nrDNA arrays. Therefore, foreign nrDNA in Hordeum likely poses a dead-end by-product of horizontal gene transfer events.

Molecular Evolution and Organization of Ribosomal DNA in the Hawkweed Tribe Hieraciinae (Cichorieae, Asteraceae)

Molecular evolution of ribosomal DNA can be highly dynamic. Hundreds to thousands of copies in the genome are subject to concerted evolution, which homogenizes sequence variants to different degrees. If well homogenized, sequences are suitable for phylogeny reconstruction; if not, sequence polymorphism has to be handled appropriately. Here we investigate non-coding rDNA sequences (ITS/ETS, 5S-NTS) along with the chromosomal organization of their respective loci (45S and 5S rDNA) in diploids of the Hieraciinae. The subtribe consists of genera Hieracium, Pilosella, Andryala, and Hispidella and has a complex evolutionary history characterized by ancient intergeneric hybridization, allele sharing among species, and incomplete lineage sorting. Direct or cloned Sanger sequences and phased alleles derived from Illumina genome sequencing were subjected to phylogenetic analyses. Patterns of homogenization and tree topologies based on the three regions were compared. In contrast to most other plant groups, 5S-NTS sequences were generally better homogenized than ITS and ETS sequences. A novel case of ancient intergeneric hybridization between Hispidella and Hieracium was inferred, and some further incongruences between the trees were found, suggesting independent evolution of these regions. In some species, homogenization of ITS/ETS and 5S-NTS sequences proceeded in different directions although the 5S rDNA locus always occurred on the same chromosome with one 45S rDNA locus. The ancestral rDNA organization in the Hieraciinae comprised 4 loci of 45S rDNA in terminal positions and 2 loci of 5S rDNA in interstitial positions per diploid genome. In Hieracium, some deviations from this general pattern were found (3, 6, or 7 loci of 45S rDNA; three loci of 5S rDNA). Some of these deviations concerned intraspecific variation, and most of them occurred at the tips of the tree or independently in different lineages. This indicates that the organization of rDNA loci is more dynamic than the evolution of sequences contained in them and that locus number is therefore largely unsuitable to inform about species relationships in Hieracium. No consistent differences in the degree of sequence homogenization and the number of 45S rDNA loci were found, suggesting interlocus concerted evolution.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.

Evolutionary history and genetic diversity of apomictic allopolyploids in Hieracium s.str.: morphological versus genomic features

PREMISE

The origin of allopolyploids is believed to shape their evolutionary potential, ecology, and geographical ranges. Morphologically distinct apomictic types sharing the same parental species belong to the most challenging groups of polyploids. We evaluated the origins and variation of two triploid taxa (Hieracium pallidiflorum, H. picroides) presumably derived from the same diploid parental pair (H. intybaceum, H. prenanthoides).

METHODS

We used a suite of approaches ranging from morphological, phylogenetic (three unlinked molecular markers), and cytogenetic analyses (in situ hybridization) to genome size screening and genome skimming.

RESULTS

Genotyping proved the expected parentage of all analyzed accessions of H. pallidiflorum and H. picroides and revealed that nearly all of them originated independently. Genome sizes and genome dosage largely corresponded to morphology, whereas the maternal origin of the allopolyploids had no discernable effect. Polyploid accessions of both parental species usually contained genetic material from other species. Given the phylogenetic distance of the parents, their chromosomes appeared only weakly differentiated in genomic in situ hybridization (GISH), as well as in overall comparisons of the repetitive fraction of their genomes. Furthermore, the repeatome of a phylogenetically more closely related species (H. umbellatum) differed significantly more.

CONCLUSIONS

We proved (1) multiple origins of hybridogeneous apomicts from the same diploid parental taxa, and (2) allopolyploid origins of polyploid accessions of the parental species. We also showed that the evolutionary dynamics of very fast evolving markers such as satellite DNA or transposable elements does not necessarily follow patterns of speciation.

Evidence of Intra-individual SSU Polymorphisms in Dark-spored Myxomycetes (Amoebozoa)

The nuclear small subunit rRNA gene (SSU or 18S) is a marker frequently used in phylogenetic and barcoding studies in Amoebozoa, including Myxomycetes. Despite its common usage and the confirmed existence of divergent copies of ribosomal genes in other protists, the potential presence of intra-individual SSU variability in Myxomycetes has never been studied before. Here we investigated the pattern of nucleotide polymorphism in the 5' end fragment of SSU by cloning and sequencing a total of 238 variants from eight specimens, each representing a species of the dark-spored orders Stemonitidales and Physarales. After excluding singletons, a relatively low SSU intra-individual variability was found but our data indicate that this might be a widely distributed phenomenon in Myxomycetes as all samples analyzed possessed various ribotypes. To determine if the occurrence of multiple SSU variants within a single specimen has a negative effect on the circumscription of species boundaries, we conducted phylogenetic analyses that revealed that clone variation may be detrimental for inferring phylogenetic relationships among some of the specimens analyzed. Despite that intra-individual variability should be assessed in additional taxa, our results indicate that special care should be taken for species identification when working with closely related species.

Reproductive pathways in Hieracium s.s. (Asteraceae): strict sexuality in diploids and apomixis in polyploids

Background and Aims

Apomixis or asexual seed reproduction is a key evolutionary mechanism in certain angiosperms providing them with reproductive assurance and isolation. Nevertheless, the frequency of apomixis is largely unknown, especially in groups with autonomous apomixis such as the diploid-polyploid genus Hieracium.

Methods

Using flow cytometric analyses, we determined the ploidy level and reproductive pathways (sexual vs. apomictic) for 7616 seeds originating from 946 plants belonging to >50 taxa sampled at 130 sites across Europe.

Key Results

Diploid seeds produced by diploids were formed exclusively by the sexual pathway after double fertilization of reduced embryo sacs. An absolute majority of tri- and tetraploid seeds (99.6 %) produced by tri- and tetrapolyploid taxa were formed by autonomous apomixis. Only 20 polyploid seeds (0.4 %) were formed sexually. These seeds, which originated on seven polyploid accessions of four taxa, were formed after fertilization of either unreduced embryo sacs through a so-called triploid bridge or reduced embryo sacs, and frequently resulted in progeny with an increased ploidy. In addition, the formation of seedlings with increased ploidy (4x and 6x) was found in two triploid plants. This is the first firm evidence on functional facultative apomixis in polyploid members of Hieracium sensu stricto (s.s.).

Conclusions

The mode of reproduction in Hieracium s.s. is tightly associated with ploidy. While diploids produce seeds exclusively sexually, polyploids produce seeds by obligate or almost obligate apomixis. Strict apomixis can increase the reproductive assurance and this in turn can increase the colonization ability of apomicts. Nevertheless, our data clearly show that certain polyploid plants are still able to reproduce sexually and contribute to the formation of new cytotypes and genotypes. The finding of functional facultative apomicts is essential for future studies focused on evolution, inheritance and ecological significance of apomixis in this genus.

The endangered Florida pondweed (Potamogeton floridanus) is a hybrid: Why we need to understand biodiversity thoroughly

Thorough understanding of biodiversity is a fundamental prerequisite for biological research. A lack of taxonomic knowledge and species misidentifications are particularly critical for conservation. Here we present an example of Potamogeton floridanus, the Florida Pondweed, an endangered taxon endemic to a small area in the Florida panhandle, whose taxonomic status remained controversial for more than a century, and all previous attempts to elucidate its identity have failed. We applied molecular approaches to tackle the origin of the mysterious taxon and supplemented them with morphological and anatomical investigations of both historical herbarium collections and plants recently collected in the type area for a comprehensive taxonomic reassessment. Sequencing of two nuclear ribosomal markers and one chloroplast non-coding spacer resulted in the surprising discovery that P. floridanus is a hybrid of P. pulcher and P. oakesianus, with the former being the maternal parent. The hybrid colony is currently geographically isolated from the distribution range of P. oakesianus. We show that previous molecular analyses have failed to reveal its hybrid identity due to inadequate nuclear DNA sequence editing. This is an example how the uncritical use of automized sequence reads can hamper molecular species identifications and also affect phylogenetic tree construction and interpretation. This unique hybrid taxon, P. ×floridanus, adds another case study to the debate on hybrid protection; consequences for its conservation are discussed.

Methanol Extracts of 28 Hieracium Species from the Balkan Peninsula - Comparative LC-MS Analysis, Chemosystematic Evaluation of their Flavonoid and Phenolic Acid Profiles and Antioxidant Potentials

INTRODUCTION

Hieracium s. str. represents one of the largest and most complex genera of flowering plants. As molecular genetics seems unlikely to disentangle intricate relationships within this reticulate species complex, analysis of flavonoids and phenolic acids, known as good chemosystematic markers, promise to be more reliable. Data about pharmacological activity of Hieracium species are scarce.

OBJECTIVE

Evaluation of the chemosystematic significance of flavonoids and phenolic acids of methanol extracts of aerial flowering parts of 28 Hieracium species from the Balkans. Additionally, investigation of antioxidant potentials of the extracts.

METHODS

Comparative qualitative and quantitative analysis of flavonoids and phenolic acids was performed by LC-MS. Multivariate statistical data analysis included non-metric multidimensional scaling (nMDS), unweighted pair-group arithmetic averages (UPGMA) and principal component analysis (PCA). Antioxidant activity was evaluated using three colorimetric tests.

RESULTS

Dominant phenolics in almost all species were luteolin type flavonoids, followed by phenolic acids. Although the investigated Hieracium species share many compounds, the current classification of the genus was supported by nMDS and UPGMA analyses with a good resolution to the group level. Hieracium naegelianum was clearly separated from the other investigated species. Spatial and ecological distances of the samples were likely to influence unexpected differentiation of some groups within H. sect. Pannosa. The vast majority of dominant compounds significantly contributed to differences between taxa. The antioxidant potential of the extracts was satisfactory and in accordance with their phenolics composition.

CONCLUSIONS

Comparative LC-MS analysis demonstrated that flavonoids and phenolic acids are good indicators of chemosystematic relationships within Hieracium, particularly between non-hybrid species and groups from the same location. Copyright © 2017 John Wiley & Sons, Ltd.

Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb

In all eukaryotes, the highly repeated 35S ribosomal DNA (rDNA) sequences encoding 18S-5.8S-26S ribosomal RNA (rRNA) typically show high levels of intragenomic uniformity due to homogenisation processes, leading to concerted evolution of 35S rDNA repeats. Here, we compared 35S rDNA divergence in several seed plants using next generation sequencing and a range of molecular and cytogenetic approaches. Most species showed similar 35S rDNA homogeneity indicating concerted evolution. However, Cycas revoluta exhibits an extraordinary diversity of rDNA repeats (nucleotide sequence divergence of different copies averaging 12 %), influencing both the coding and non-coding rDNA regions nearly equally. In contrast, its rRNA transcriptome was highly homogeneous suggesting that only a minority of genes (<20 %) encode functional rRNA. The most common SNPs were C > T substitutions located in symmetrical CG and CHG contexts which were also highly methylated. Both functional genes and pseudogenes appear to cluster on chromosomes. The extraordinary high levels of 35S rDNA diversity in C. revoluta, and probably other species of cycads, indicate that the frequency of repeat homogenisation has been much lower in this lineage, compared with all other land plant lineages studied. This has led to the accumulation of methylation-driven mutations and pseudogenisation. Potentially, the reduced homology between paralogs prevented their elimination by homologous recombination, resulting in long-term retention of rDNA pseudogenes in the genome.

Understanding diploid diversity: A first step in unraveling polyploid, apomictic complexity in Amelanchier

PREMISE OF THE STUDY

Delimitation of Amelanchier species is difficult because of polyploidy and gametophytic apomixis. A first step in unraveling this species problem is understanding the diversity of the diploids that contributed genomes to polyploid apomicts. This research helps clarify challenging species-delimitation problems attending polyploid, apomictic complexity.

METHODS

We sampled 431 diploid accessions from 13 species, of which 10 are North American and three are Old World. Quantitative morphological analyses tested the null hypothesis of no discrete groups. Using three to nine diploid accessions per species, we constructed phylogenies with DNA sequences from ETS, ITS, the second intron of LEAFY, and chloroplast regions rpoB-trnC, rpl16, trnD-trnT, and ycf6-psbM.

KEY RESULTS

Most Amelanchier diploid taxa are morphologically and ecogeographically distinct and genetically exclusive lineages. They rarely hybridize with one another. Nuclear and chloroplast DNA sequences almost completely resolve the Amelanchier phylogeny. The backbone is the mostly western North American clade A, eastern North American clade B, and Old World clade O. DNA sequences and morphology support clades A and O as sister taxa. Despite extensive paralogy, our LEAFY data are phylogenetically informative and identify a clade (T) of three arborescent taxa within clade B.

CONCLUSIONS

Amelanchier diploids differ strikingly from polyploid apomicts, in that hybridization among them is rare, and they form taxa that would qualify as species by most species concepts. Knowledge of diploid morphology, phylogeny, and ecogeography provides a foundation for understanding the evolutionary history of polyploid apomicts, their patterns of diversification, and their species status.

Evolution of apomixis loci in Pilosella and Hieracium (Asteraceae) inferred from the conservation of apomixis-linked markers in natural and experimental populations

The Hieracium and Pilosella (Lactuceae, Asteraceae) genera of closely related hawkweeds contain species with two different modes of gametophytic apomixis (asexual seed formation). Both genera contain polyploid species, and in wild populations, sexual and apomictic species co-exist. Apomixis is known to co-exist with sexuality in apomictic Pilosella individuals, however, apomictic Hieracium have been regarded as obligate apomicts. Here, a developmental analysis of apomixis within 16 Hieracium species revealed meiosis and megaspore tetrad formation in 1 to 7% of ovules, for the first time indicating residual sexuality in this genus. Molecular markers linked to the two independent, dominant loci LOSS OF APOMEIOSIS (LOA) and LOSS OF PARTHENOGENESIS (LOP) controlling apomixis in Pilosella piloselloides subsp. praealta were screened across 20 phenotyped Hieracium individuals from natural populations, and 65 phenotyped Pilosella individuals from natural and experimental cross populations, to examine their conservation, inheritance and association with reproductive modes. All of the tested LOA and LOP-linked markers were absent in the 20 Hieracium samples irrespective of their reproductive mode. Within Pilosella, LOA and LOP-linked markers were essentially absent within the sexual plants, although they were not conserved in all apomictic individuals. Both loci appeared to be inherited independently, and evidence for additional genetic factors influencing quantitative expression of LOA and LOP was obtained. Collectively, these data suggest independent evolution of apomixis in Hieracium and Pilosella and are discussed with respect to current knowledge of the evolution of apomixis.

Reticulate evolution in North American black-fruited hawthorns (Crataegus section Douglasia; Rosaceae): evidence from nuclear ITS2 and plastid sequences

BACKGROUND AND AIMS

The taxonomic complexity of Crataegus (hawthorn; Rosaceae, Maleae), especially in North America, has been attributed by some to hybridization in combination with gametophytic apomixis and polyploidization, whereas others have considered the roles of hybridization and apomixis to be minimal. Study of the chemical composition and therapeutic value of hawthorn extracts requires reproducible differentiation of entities that may be difficult to distinguish by morphology alone. This study sought to address this by using the nuclear ribosomal spacer region ITS2 as a supplementary DNA barcode; however, a lack of success prompted an investigation to discover why this locus gave unsatisfactory results.

METHODS

ITS2 was extensively cloned so as to document inter- and intraindividual variation in this locus, using hawthorns of western North America where the genus Crataegus is represented by only two widely divergent groups, the red-fruited section Coccineae and the black-fruited section Douglasia. Additional sequence data from selected loci on the plastid genome were obtained to enhance further the interpretation of the ITS2 results.

KEY RESULTS

In the ITS2 gene tree, ribotypes from western North American hawthorns are found in two clades. Ribotypes from diploid members of section Douglasia occur in one clade (with representatives of the east-Asian section Sanguineae). The other clade comprises those from diploid and polyploid members of section Coccineae. Both clades contribute ribotypes to polyploid Douglasia. Data from four plastid-derived intergenic spacers demonstrate the maternal parentage of these allopolyploids.

CONCLUSIONS

Repeated hybridization between species of section Douglasia and western North American members of section Coccineae involving the fertilization of unreduced female gametes explains the observed distribution of ribotypes and accounts for the phenetic intermediacy of many members of section Douglasia.

High penetrance of a pan-canina type rDNA family in intersection Rosa hybrids suggests strong selection of bivalent chromosomes in the section Caninae

All dogroses (Rosa sect. Caninae) are characterized by the peculiar canina meiosis in which genetic material is unevenly distributed between female and male gametes. The pan-canina rDNA family (termed beta) appears to be conserved in all dogroses analyzed so far. Here, we have studied rDNAs in experimental hybrids obtained from open pollination of F1 plants derived from 2 independent intersectional crosses between the pentaploid dogrose species (2n = 5x = 35) Rosa rubiginosa as female parent (producing 4x egg cells due to the unique asymmetrical canina meiosis) and the tetraploid (2n = 4x = 28) garden rose R. hybrida 'André Brichet' as male parent (producing 2x pollen after normal meiosis). We analyzed the structure of rDNA units by molecular methods [CAPS and extensive sequencing of internal transcribed spacers (ITS)] and determined the number of loci on chromosomes by FISH. FISH showed that R. rubiginosa and 'André Brichet' harbored 5 and 4 highly heteromorphic rDNA loci, respectively. In the second generation of hybrid lines, we observed a reduced number of loci (4 and 5 instead of the expected 6). In R. rubiginosa and 'André Brichet', 2-3 major ITS types were found which is consistent with a weak homogenization pressure maintaining high diversity of ITS types in this genus. In contrast to expectation (the null hypothesis of Mendelian inheritance of ITS families), we observed reduced ITS diversity in some individuals of the second generation which might derive from self-fertilization or from a backcross to R. rubiginosa. In these individuals, the pan-canina beta family appeared to be markedly enriched, while the paternal families were lost or diminished in copies. Although the mechanism of biased meiotic transmission of certain rDNA types is currently unknown, we speculate that the bivalent-forming chromosomes carrying the beta rDNA family exhibit extraordinary pairing efficiency and/or are subjected to strong selection in Caninae polyploids.

Evolution of the polyploid north-west Iberian Leucanthemum pluriflorum clan (Compositae, Anthemideae) based on plastid DNA sequence variation and AFLP fingerprinting

BACKGROUND AND AIMS

The genus Leucanthemum is a species-rich polyploid complex from southern and central Europe, comprising 41 species with ploidy ranging from 2x to 22x. The present contribution aims at reconstructing the evolutionary history of a geographically isolated species group (the L. pluriflorum clan) from the north-west Iberian Peninsula comprising the diploid L. pluriflorum, the tetraploids L. ircutianum subsp. pseudosylvaticum and L. × corunnense (a putative hybrid taxon based on crossing between L. pluriflorum and L. merinoi), and the hexaploids L. sylvaticum and L. merinoi.

METHODS

Chromosome number variation (determined flow cytometrically) and sequence variation were analysed for two intergenic spacer regions on the plastid genome (psbA-trnH and trnC-petN) for individuals from 54 populations in combination with amplified fragment length polymorphism (AFLP) fingerprinting of 246 representative individuals from these populations.

KEY RESULTS

Plastid sequence data revealed that all surveyed members of the L. pluriflorum clan possess plastid haplotypes that are closely related to each other and distinctly separated from other Leucanthemum species. AFLP fingerprinting resulted in allopolyploid fragment patterns for most of the polyploid populations, except for the tetraploid L. × corunnense and a further tetraploid population in northern Galicia, which cluster with the diploids rather than with the other polyploids. In silico modelling of (auto)tetraploid AFLP genotypes further corroborates the allopolyploid nature of L. ircutianum subsp. pseudosylvaticum, L. sylvaticum and L. merinoi.

CONCLUSIONS

The present study provides evidence for recognizing one diploid (L. pluriflorum), one autotetraploid (L. corunnense), one allotetraploid (L. pseudosylvaticum) and one allohexaploid (L. sylvaticum with the two geographically and ecologically differentiated subspecies subsp. sylvaticum and subsp. merinoi) in the L. pluriflorum clan. It also has implications for the understanding of biogeographical patterns in the Iberian Peninsula.

Retreating or standing: responses of forest species and steppe species to climate change in arid Eastern Central Asia

Background

The temperature in arid Eastern Central Asia is projected to increase in the future, accompanied by increased variability of precipitation. To investigate the impacts of climate change on plant species in this area, we selected two widespread species as candidates, Clematis sibirica and C. songorica, from montane coniferous forest and arid steppe habitats respectively.

Methodology/Principal Findings

We employed a combined approach of molecular phylogeography and species distribution modelling (SDM) to predict the future responses of these two species to climate change, utilizing evidence of responses from the past. Genetic data for C. sibirica shows a significant phylogeographical signal (N_ST > F_ST, P<0.05) and demographic contraction during the glacial-interglacial cycles in the Pleistocene. This forest species would likely experience range reduction, though without genetic loss, in the face of future climate change. In contrast, SDMs predict that C. songorica, a steppe species, should maintain a consistently stable potential distribution under the Last Glacial Maximum (LGM) and the future climatic conditions referring to its existing potential distribution. Molecular results indicate that the presence of significant phylogeographical signal in this steppe species is rejected and this species contains a high level of genetic differentiation among populations in cpDNA, likely benefiting from stable habitats over a lengthy time period.

Conclusions/Significance

Evidence from the molecular phylogeography of these two species, the forest species is more sensitive to past climate changes than the steppe species. SDMs predict that the forest species will face the challenge of potential range contraction in the future more than the steppe species. This provides a perspective on ecological management in arid Eastern Central Asia, indicating that increased attention should be paid to montane forest species, due to their high sensitivity to disturbance.

Patterns of plastid and nuclear variation among apomictic polyploids of Hieracium: evolutionary processes and taxonomic implications

BACKGROUND AND AIMS

Apomictic species (with asexual seed production) make up for 20-50 % of all taxonomically recognized species in northern Europe, but the phylogenetic relationships of apomictic species and the mode of evolution and speciation remain largely unknown and their taxonomy is consequently disputed.

METHODS

In the present study, plastid psbD-trnT sequences (349 accessions) and 12 nuclear microsatellite loci (478 accessions) were used to create an overview of the molecular variation in (mainly) northern European members of the most species-rich of all plant genera, Hieracium s.s. The results are discussed and interpreted in the context of morphological and cytological data on the same species.

KEY RESULTS AND CONCLUSIONS

The complete psbD-trnT alignment was 1243 bp and 50 polymorphisms defined 40 haplotypes. All haplotypes found in the sections of the genus distributed in the northern European lowlands fell into one of two main groups, group H and group V, mutually separated by seven or eight polymorphisms. All accessions belonging to H. sects. Foliosa, Hieracioides (viz. H. umbellatum) and Tridentata and all but one accession of triploid species of H. sects. Oradea and Vulgata showed haplotypes of group V. Haplotypes of group H were found in all accessions of H. sects. Bifida and Hieracium and in all tetraploid representatives of H. sects. Oreadea and Vulgata. Additional haplotypes were found in accessions of the genus Pilosella and in southern European and Alpine sections of Hieracium. In contrast, the distribution of individual haplotypes in the two major groups appeared uncorrelated with morphology and current taxonomy, but polymorphisms within species were only rarely encountered. In total, 160 microsatellite alleles were identified. Levels of variation were generally high with only nine pairs of accessions being identical at all loci (in all cases representing accessions of the same species). In the neighbor-joining analysis based on the microsatellite data, accessions of the same species generally clustered together and some smaller groups of species congruent with morphology and/or current taxonomy were recovered but, except for H. sect. Oreadea, most larger groups were not correlated with morphology. Although the plastid DNA sequences show too little variation and the nuclear microsatellites are too variable to resolve relationships successfully among species or to fully understand processes of evolution, it is concluded that both species and sections as defined by morphology are largely congruent with the molecular data, that gene flow between the sections is rare or non-existent and that the tetraploid species may constitute the key to understanding evolution and speciation in this genus.

Reconstruction of phylogenetic relationships in a highly reticulate group with deep coalescence and recent speciation (Hieracium, Asteraceae)

Phylogeny reconstruction based on multiple unlinked markers is often hampered by incongruent gene trees, especially in closely related species complexes with high degrees of hybridization and polyploidy. To investigate the particular strengths and limitations of chloroplast DNA (cpDNA), low-copy nuclear and multicopy nuclear markers for elucidating the evolutionary history of such groups, we focus on Hieracium s.str., a predominantly apomictic genus combining the above-mentioned features. Sequences of the trnV-ndhC and trnT-trnL intergenic spacers were combined for phylogenetic analyses of cpDNA. Part of the highly variable gene for squalene synthase (sqs) was applied as a low-copy nuclear marker. Both gene trees were compared with previous results based on the multicopy external transcribed spacer (ETS) of the nuclear ribosomal DNA. The power of the different markers to detect hybridization varied, but they largely agreed on particular hybrid and allopolyploid origins. The same crown groups of species were recognizable in each dataset, but basal relationships were strongly incongruent among cpDNA, sqs and ETS trees. The ETS tree was considered as the best approximation of the species tree. Both cpDNA and sqs trees showed basal polytomies as well as merging or splitting of species groups of non-hybrid taxa. These patterns can be best explained by a rapid diversification of the genus with ancestral polymorphism and incomplete lineage sorting. A hypothetical scenario of Hieracium speciation based on all available (including non-molecular) evidence is depicted. Incorporation of seemingly contradictory information helped to better understand species origins and evolutionary patterns in this notoriously difficult agamic complex.

Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids

BACKGROUND

Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units.

METHODS

We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci.

RESULTS AND CONCLUSIONS

In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.

Intraradical dynamics of two coexisting isolates of the arbuscular mycorrhizal fungus Glomus intraradices sensu lato as estimated by real-time PCR of mitochondrial DNA

Real-time PCR in nuclear ribosomal DNA (nrDNA) is becoming a well-established tool for the quantification of arbuscular mycorrhizal (AM) fungi, but this genomic region does not allow the specific amplification of closely related genotypes. The large subunit of mitochondrial DNA (mtDNA) has a higher-resolution power, but mtDNA-based quantification has not been previously explored in AM fungi. We applied real-time PCR assays targeting the large subunit of mtDNA to monitor the DNA dynamics of two isolates of Glomus intraradices sensu lato coexisting in the roots of medic (Medicago sativa). The mtDNA-based quantification was compared to quantification in nrDNA. The ratio of copy numbers determined by the nrDNA- and mtDNA-based assays consistently differed between the two isolates. Within an isolate, copy numbers of the nuclear and the mitochondrial genes were closely correlated. The two quantification approaches revealed similar trends in the dynamics of both isolates, depending on whether they were inoculated alone or together. After 12 weeks of cultivation, competition between the two isolates was observed as a decrease in the mtDNA copy numbers of one of them. The coexistence of two closely related isolates, which cannot be discriminated by nrDNA-based assays, was thus identified as a factor influencing the dynamics of AM fungal DNA in roots. Taken together, the results of this study show that real-time PCR assays targeted to the large subunit of mtDNA may become useful tools for the study of coexisting AM fungi.

Polyploid evolution and Pleistocene glacial cycles: A case study from the alpine primrose Primula marginata (Primulaceae)

Background

Recent studies highlighted the role of Pleistocene climatic cycles in polyploid speciation and of southern Alpine refugia as reservoirs of diversity during glacial maxima. The polyploid Primula marginata, endemic to the southwestern Alps, includes both hexaploid and dodecaploid cytotypes that show no ecological or morphological differences. We used flow cytometry to determine variation and geographic distribution of cytotypes within and between populations and analyses of chloroplast (cp) and nuclear ribosomal (nr) DNA sequences from the Internal Transcribed Spacer (ITS) region to infer the evolutionary history of the two cytotypes and the auto- vs. allopolyploid origin of dodecaploid populations.

Results

We did not detect any intermediate cytotypes or variation of ploidy levels within populations. Hexaploids occur in the western and dodecaploids in the eastern part of the distributional range, respectively. The cpDNA and nrDNA topologies are in conflict, for the former supports shared ancestry between P. marginata and P. latifolia, while the latter implies common origins between at least some ITS clones of P. marginata and P. allionii.

Conclusions

Our results suggest an initial episode of chloroplast capture involving ancestral lineages of P. latifolia and P. marginata, followed by polyploidization between P. marginata-like and P. allionii-like lineages in a southern refugium of the Maritime Alps. The higher proportion of ITS polymorphisms in dodecaploid than in hexaploid accessions of P. marginata and higher total nucleotide diversity of ITS clones in dodecaploid vs. hexaploid individuals sequences are congruent with the allopolyploid hypothesis of dodecaploid origin.

Development of novel low-copy nuclear markers for Hieraciinae (Asteraceae) and their perspective for other tribes

PREMISE OF THE STUDY

The development of three low-copy nuclear markers for low taxonomic level phylogenies in Asteraceae with emphasis on the subtribe Hieraciinae is reported.

METHODS AND RESULTS

Marker candidates were selected by comparing a Lactuca complementary DNA (cDNA) library with public DNA sequence databases. Interspecific variation and phylogenetic signal of the selected genes were investigated for diploid taxa from the subtribe Hieraciinae and compared to a reference phylogeny. Their ability to cross-amplify was assessed for other Asteraceae tribes. All three markers had higher variation (2.1-4.5 times) than the internal transcribed spacer (ITS) in Hieraciinae. Cross-amplification was successful in at least seven other tribes of the Asteraceae. Only three cases indicating the presence of paralogs or pseudogenes were detected.

CONCLUSIONS

The results demonstrate the potential of these markers for phylogeny reconstruction in the Hieraciinae as well as in other Asteraceae tribes, especially for very closely related species.

Long-term tracing of Rhizophagus irregularis isolate BEG140 inoculated on Phalaris arundinacea in a coal mine spoil bank, using mitochondrial large subunit rDNA markers

During the last decade, the application of arbuscular mycorrhizal fungi (AMF) as bioenhancers has increased significantly. However, until now, it has been difficult to verify the inoculation success in terms of fungal symbiont establishment in roots of inoculated plants because specific fungal strains could not be detected within colonized roots. Using mitochondrial large subunit ribosomal DNA, we show that Rhizophagus irregularis (formerly known as Glomus intraradices) isolate BEG140 consists of two different haplotypes. We developed nested PCR assays to specifically trace each of the two haplotypes in the roots of Phalaris arundinacea from a field experiment in a spoil bank of a former coal mine, where BEG140 was used as inoculant. We revealed that despite the relatively high diversity of native R. irregularis strains, R. irregularis BEG140 survived and proliferated successfully in the field experiment and was found significantly more often in the inoculated than control plots. This work is the first one to show tracing of an inoculated AMF isolate in the roots of target plants and to verify its survival and propagation in the field. These results will have implications for basic research on the ecology of AMF at the intraspecific level as well as for commercial users of mycorrhizal inoculation.

Allopolyploid origin of highly invasive Centaurea stoebe s.l. (Asteraceae)

Spotted knapweed (Centaurea stoebe) occurs from Western Asia to Western Europe both as diploid and tetraploid cytotypes, predominantly in single-cytotype populations with higher frequency of diploid populations. Interestingly, only tetraploids have been recorded so far from its introduced range in North America where they became highly invasive. We performed phylogenetic and network analyses of more than 40 accessions of the C. stoebe and C. paniculata groups and other related taxa using cloned internal transcribed spacer (ITS) and sequences of the chloroplast trnT-trnL and atpBrbcL regions to (i) assess the evolutionary origin of tetraploid C. stoebe s.l., and (ii) uncover the phylogeny of the C. stoebe group. Both issues have not been studied so far and thus remained controversial. Cloned ITS sequences showed the presence of two slightly divergent ribotypes occurring in tetraploid cytotype, while only one major ribotype was present in diploid C. stoebe s.str. This pattern suggests an allopolyploid origin of tetraploids with contribution of the diploid C. stoebe s.str. genome. Although we were not able to detect the second parental taxon, we hypothesize that hybridization might have triggered important changes in morphology and life history traits, which in turn may explain the colonization success of the tetraploid taxon. Bayesian relaxed clock estimations indicate a relatively recent--Pleistocene origin of the tetraploid C. stoebe s.l. Furthermore, our analyses showed a deep split between the C. paniculata and C. stoebe groups, and a young diversification of the taxa within the C. stoebe group. In contrast to nrDNA analyses, the observed pattern based on two cpDNA regions was inconclusive with respect to the origin and phylogeny of the studied taxa, most likely due to shared ancient polymorphism and frequent homoplasies.

On the genome constitution and evolution of intermediate wheatgrass (Thinopyrum intermedium: Poaceae, Triticeae)

Background

The wheat tribe Triticeae (Poaceae) is a diverse group of grasses representing a textbook example of reticulate evolution. Apart from globally important grain crops, there are also wild grasses which are of great practical value. Allohexaploid intermediate wheatgrass, Thinopyrum intermedium (2n = 6x = 42), possesses many desirable agronomic traits that make it an invaluable source of genetic material useful in wheat improvement. Although the identification of its genomic components has been the object of considerable investigation, the complete genomic constitution and its potential variability are still being unravelled. To identify the genomic constitution of this allohexaploid, four accessions of intermediate wheatgrass from its native area were analysed by sequencing of chloroplast trnL-F and partial nuclear GBSSI, and genomic in situ hybridization.

Results

The results confirmed the allopolyploid origin of Thinopyrum intermedium and revealed new aspects in its genomic composition. Genomic heterogeneity suggests a more complex origin of the species than would be expected if it originated through allohexaploidy alone. While Pseudoroegneria is the most probable maternal parent of the accessions analysed, nuclear GBSSI sequences suggested the contribution of distinct lineages corresponding to the following present-day genera: Pseudoroegneria, Dasypyrum, Taeniatherum, Aegilops and Thinopyrum. Two subgenomes of the hexaploid have most probably been contributed by Pseudoroegneria and Dasypyrum, but the identity of the third subgenome remains unresolved satisfactorily. Possibly it is of hybridogenous origin, with contributions from Thinopyrum and Aegilops. Surprising diversity of GBSSI copies corresponding to a Dasypyrum-like progenitor indicates either multiple contributions from different sources close to Dasypyrum and maintenance of divergent copies or the presence of divergent paralogs, or a combination of both. Taeniatherum-like GBSSI copies are most probably pseudogenic, and the mode of their acquisition by Th. intermedium remains unclear.

Conclusions

Hybridization has played a key role in the evolution of the Triticeae. Transfer of genetic material via extensive interspecific hybridization and/or introgression could have enriched the species' gene pools significantly. We have shown that the genomic heterogeneity of intermediate wheatgrass is higher than has been previously assumed, which is of particular concern to wheat breeders, who frequently use it as a source of desirable traits in wheat improvement.

A next-generation sequencing method for overcoming the multiple gene copy problem in polyploid phylogenetics, applied to Poa grasses

BACKGROUND

Polyploidy is important from a phylogenetic perspective because of its immense past impact on evolution and its potential future impact on diversification, survival and adaptation, especially in plants. Molecular population genetics studies of polyploid organisms have been difficult because of problems in sequencing multiple-copy nuclear genes using Sanger sequencing. This paper describes a method for sequencing a barcoded mixture of targeted gene regions using next-generation sequencing methods to overcome these problems.

RESULTS

Using 64 3-bp barcodes, we successfully sequenced three chloroplast and two nuclear gene regions (each of which contained two gene copies with up to two alleles per individual) in a total of 60 individuals across 11 species of Australian Poa grasses. This method had high replicability, a low sequencing error rate (after appropriate quality control) and a low rate of missing data. Eighty-eight percent of the 320 gene/individual combinations produced sequence reads, and >80% of individuals produced sufficient reads to detect all four possible nuclear alleles of the homeologous nuclear loci with 95% probability.We applied this method to a group of sympatric Australian alpine Poa species, which we discovered to share an allopolyploid ancestor with a group of American Poa species. All markers revealed extensive allele sharing among the Australian species and so we recommend that the current taxonomy be re-examined. We also detected hypermutation in the trnH-psbA marker, suggesting it should not be used as a land plant barcode region. Some markers indicated differentiation between Tasmanian and mainland samples. Significant positive spatial genetic structure was detected at <100 km with chloroplast but not nuclear markers, which may be a result of restricted seed flow and long-distance pollen flow in this wind-pollinated group.

CONCLUSIONS

Our results demonstrate that 454 sequencing of barcoded amplicon mixtures can be used to reliably sample all alleles of homeologous loci in polyploid species and successfully investigate phylogenetic relationships among species, as well as to investigate phylogeographic hypotheses. This next-generation sequencing method is more affordable than and at least as reliable as bacterial cloning. It could be applied to any experiment involving sequencing of amplicon mixtures.

Similar patterns of rDNA evolution in synthetic and recently formed natural populations of Tragopogon (Asteraceae) allotetraploids

BACKGROUND

Tragopogon mirus and T. miscellus are allotetraploids (2n = 24) that formed repeatedly during the past 80 years in eastern Washington and adjacent Idaho (USA) following the introduction of the diploids T. dubius, T. porrifolius, and T. pratensis (2n = 12) from Europe. In most natural populations of T. mirus and T. miscellus, there are far fewer 35S rRNA genes (rDNA) of T. dubius than there are of the other diploid parent (T. porrifolius or T. pratensis). We studied the inheritance of parental rDNA loci in allotetraploids resynthesized from diploid accessions. We investigate the dynamics and directionality of these rDNA losses, as well as the contribution of gene copy number variation in the parental diploids to rDNA variation in the derived tetraploids.

RESULTS

Using Southern blot hybridization and fluorescent in situ hybridization (FISH), we analyzed copy numbers and distribution of these highly reiterated genes in seven lines of synthetic T. mirus (110 individuals) and four lines of synthetic T. miscellus (71 individuals). Variation among diploid parents accounted for most of the observed gene imbalances detected in F1 hybrids but cannot explain frequent deviations from repeat additivity seen in the allotetraploid lines. Polyploid lineages involving the same diploid parents differed in rDNA genotype, indicating that conditions immediately following genome doubling are crucial for rDNA changes. About 19% of the resynthesized allotetraploid individuals had equal rDNA contributions from the diploid parents, 74% were skewed towards either T. porrifolius or T. pratensis-type units, and only 7% had more rDNA copies of T. dubius-origin compared to the other two parents. Similar genotype frequencies were observed among natural populations. Despite directional reduction of units, the additivity of 35S rDNA locus number is maintained in 82% of the synthetic lines and in all natural allotetraploids.

CONCLUSIONS

Uniparental reductions of homeologous rRNA gene copies occurred in both synthetic and natural populations of Tragopogon allopolyploids. The extent of these rDNA changes was generally higher in natural populations than in the synthetic lines. We hypothesize that locus-specific and chromosomal changes in early generations of allopolyploids may influence patterns of rDNA evolution in later generations.

Molecular phylogenetic insights into the evolution of Octocorallia: a review

The anthozoan sub-class Octocorallia, comprising approximately 3000 species of soft corals, gorgonians, and sea pens, remains one of the most poorly understood groups of the phylum Cnidaria. Efforts to classify the soft corals and gorgonians at the suprafamilial level have long thwarted taxonomists, and the subordinal groups in current use are widely recognized to represent grades of colony forms rather than clades. Molecular phylogenetic analyses of the sub-class do not support either the current morphologically based subordinal or familial-level taxonomy. To date, however, the resolution necessary to propose an alternative, phylogenetic classification of Octocorallia or to elucidate patterns of morphological evolution within the group is lacking. Attempts to understand boundaries between species and interspecific or intraspecific phylogenetic relationships have been hampered by the very slow rate of mitochondrial gene evolution in Octocorallia, and a consequent dearth of molecular markers with variation sufficient to distinguish species (or sometimes genera). A review of the available ITS2 sequence data for octocorals, however, reveals a yet-unexplored phylogenetic signal both at sequence and secondary-structure levels. In addition, incongruence between mitochondrial and nuclear gene trees suggests that hybrid speciation and reticulate evolution may be an important mechanism of diversification in some genera. Emerging next-generation genomic-sequencing technologies offer the best hope for a breakthrough in our understanding of phylogenetic relationships and of evolution of morphological traits in Octocorallia. Genome and transcriptome sequencing may provide enough characters to resolve relationships at the deepest levels of the octocoral tree, while simultaneously offering an efficient means to screen for new genetic markers variable enough to distinguish species and populations.