Why nuclear ribosomal DNA spacers (ITS) tell different stories in Quercus

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

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.

Lack of ITS sequence homogenization in Erysimum species (Brassicaceae) with different ploidy levels

The internal transcribed spacers (ITS) exhibit concerted evolution by the fast homogenization of these sequences at the intragenomic level. However, the rate and extension of this process are unclear and might be conditioned by the number and divergence of the different ITS copies. In some cases, such as hybrid species and polyploids, ITS sequence homogenization appears incomplete, resulting in multiple haplotypes within the same organism. Here, we studied the dynamics of concerted evolution in 85 individuals of seven plant species of the genus Erysimum (Brassicaceae) with multiple ploidy levels. We estimated the rate of concerted evolution and the degree of sequence homogenization separately for ITS1 and ITS2 and whether these varied with ploidy. Our results showed incomplete sequence homogenization, especially for polyploid samples, indicating a lack of concerted evolution in these taxa. Homogenization was usually higher in ITS2 than in ITS1, suggesting that concerted evolution operates more efficiently on the former. Furthermore, the hybrid origin of several species appears to contribute to the maintenance of high haplotype diversity, regardless of the level of ploidy. These findings indicate that sequence homogenization of ITS is a dynamic and complex process that might result in varying intra- and inter-genomic diversity levels.

Extensive nrDNA ITS polymorphism in Lycium: Non-concerted evolution and the identification of pseudogenes

The internal transcribed spacer (ITS) is one of the most extensively sequenced molecular markers in plant systematics due to its generally concerted evolution. While non-concerted evolution has been found in some plant taxa, such information is missing in Lycium. Molecular studies of six species and two variants of the genus Lycium revealed high levels of intra- and inter-individual polymorphism in the ITS, indicating non-concerted evolution. All genomic DNA ITS paralogues were identified as putative pseudogenes or functional paralogues through a series of comparisons of sequence features, including length and substitution variation, GC content, secondary structure stability, and the presence of conserved motifs in the 5.8S gene, and the rate of evolution. Approximately, 60% of ITS pseudogenes could be easily detected. Based on phylogenetic analysis, all pseudogenes were highly distinct from their corresponding functional copies, tended to evolve neutrally, and clustered randomly together in the evolutionary tree. The results probably suggest that this ITS non-concerted evolution is related to the recent divergence between tandem repeats within the Lycium genome and hybridization between species. Our study complements those of pseudogenes in plant taxa and provides a theoretical basis for the phylogeny and genetic origin of the genus Lycium while having important implications for the use of ITS molecular markers for phylogenetic reconstruction.

A new method for determining ribosomal DNA copy number shows differences between Saccharomyces cerevisiae populations

Ribosomal DNA genes (rDNA) encode the major ribosomal RNAs and in eukaryotes typically form tandem repeat arrays. Species have characteristic rDNA copy numbers, but there is substantial intra-species variation in copy number that results from frequent rDNA recombination. Copy number differences can have phenotypic consequences, however difficulties in quantifying copy number mean we lack a comprehensive understanding of how copy number evolves and the consequences. Here we present a genomic sequence read approach to estimate rDNA copy number based on modal coverage to help overcome limitations with existing mean coverage-based approaches. We validated our method using Saccharomyces cerevisiae strains with known rDNA copy numbers. Application of our pipeline to a global sample of S. cerevisiae isolates showed that different populations have different rDNA copy numbers. Our results demonstrate the utility of the modal coverage method, and highlight the high level of rDNA copy number variation within and between populations.

An Empirical Analysis Rejects the Hybrid Speciation Hypothesis of a Crucial Kiwifruit Species, Despite Genomic Evidence of Frequent Interspecific Gene Flow in the Genus

Hybrid speciation is an important way to generate species diversity. In general, however, interspecific hybridization is easily confused with the formation of hybrid species. Using the genomic resequencing data of the kiwifruit genus (Actinidia), at least ten species were documented recently as homoploid hybrid species, and thus a two-layer mode of species diversification has been proposed. As a crucial piece of evidence, Actinidia fulvicoma was identified as a hybrid derivative of Actinidia eriantha × Actinidia cylindrica, representing a rare case of hybrid species in kiwifruit that won the competition of ecological niches with one of its putative parental species, A. cylindrica. However, the hypothesized hybrid origin of A. fulvicoma is inconsistent with our specimen observations. Here, we present multiple lines of evidence to reject the hybrid speciation hypothesis for this species, despite genomic evidence for frequent interspecific gene flow. We collected the samples of A. fulvicoma in type locality and neighboring regions to contrast them with type specimen, and sequenced nuclear ribosomal DNA ITS, chloroplast trnL-trnF and mitochondrial nad2-i3, as well as four single-copy nuclear genes explored from kiwifruit genomes, to infer phylogenetic relationships among A. fulvicoma, its putative parental species, and their relatives. Our data definitely reveal that A. fulvicoma occupies an independent backbone lineage and it is not a hybrid. This study suggests that correct evolutionary applications on extensive surveys of the putative hybrid and its possible parents with strict criteria are necessary in the documentation of hybrid speciation to advance our understanding of the genomic basis of hybrid species.

Molecular phylogeny of mulberries reconstructed from ITS and two cpDNA sequences

Background

Species in the genus Morus (Moraceae) are deciduous woody plants of great economic importance. The classification and phylogenetic relationships of Morus, especially the abundant mulberry resources in China, is still undetermined. Internal transcribed spacer (ITS) regions are among the most widely used molecular markers in phylogenetic analyses of angiosperms. However, according to the previous phylogenetic analyses of ITS sequences, most of the mulberry accessions collected in China were grouped into the largest clade lacking for phylogenetic resolution. Compared with functional ITS sequences, ITS pseudogenes show higher sequence diversity, so they can provide useful phylogenetic information.

Methods

We sequenced the ITS regions and the chloroplast DNA regions TrnL-TrnF and TrnT-TrnL from 33 mulberry accessions, and performed phylogenetic analyses to explore the evolution of mulberry.

Results

We found ITS pseudogenes in 11 mulberry accessions. In the phylogenetic tree constructed from ITS sequences, clade B was separated into short-type sequence clades (clades 1 and 2), and a long-type sequence clade (clade 3). Pseudogene sequences were separately clustered into two pseudogroups, designated as pseudogroup 1 and pseudogroup 2. The phylogenetic tree generated from cpDNA sequences also separated clade B into two clades.

Conclusions

Two species were separated in clade B. The existence of three connection patterns and incongruent distribution patterns between the phylogenetic trees generated from cpDNA and ITS sequences suggested that the ITS pseudogene sequences connect with genetic information from the female progenitor. Hybridization has played important roles in the evolution of mulberry, resulting in low resolution of the phylogenetic analysis based on ITS sequences. An evolutionary pattern illustrating the evolution history of mulberry is proposed. These findings have significance for the conservation of local mulberry resources. Polyploidy, hybridization, and concerted evolution have all played the roles in the evolution of ITS sequences in mulberry. This study will expand our understanding of mulberry evolution.

Species Identification of Oaks (Quercus L., Fagaceae) from Gene to Genome

Species identification of oaks (Quercus) is always a challenge because many species exhibit variable phenotypes that overlap with other species. Oaks are notorious for interspecific hybridization and introgression, and complex speciation patterns involving incomplete lineage sorting. Therefore, accurately identifying Quercus species barcodes has been unsuccessful. In this study, we used chloroplast genome sequence data to identify molecular markers for oak species identification. Using next generation sequencing methods, we sequenced 14 chloroplast genomes of Quercus species in this study and added 10 additional chloroplast genome sequences from GenBank to develop a DNA barcode for oaks. Chloroplast genome sequence divergence was low. We identified four mutation hotspots as candidate Quercus DNA barcodes; two intergenic regions (matK-trnK-rps16 and trnR-atpA) were located in the large single copy region, and two coding regions (ndhF and ycf1b) were located in the small single copy region. The standard plant DNA barcode (rbcL and matK) had lower variability than that of the newly identified markers. Our data provide complete chloroplast genome sequences that improve the phylogenetic resolution and species level discrimination of Quercus. This study demonstrates that the complete chloroplast genome can substantially increase species discriminatory power and resolve phylogenetic relationships in plants.

Molecular evolution of the internal transcribed spacers in red oaks (Quercus sect. Lobatae)

Previous studies of the Internal Transcribed Spacers of the nuclear ribosomal DNA (ITS) in sections Quercus (white oaks), Protobalanus (intermediate or golden cup oaks), Cerris (Cerris oaks), and Ilex (Ilex oaks) suggest that ITS regions undergo full concerted evolution in oaks; however, ITS evolution patterns in red oaks (section Lobatae) are unknown due to scant representation in published work. To determine whether full concerted evolution occurs in red oaks, the purpose of this study was to examine ITS sequences from 40 red oak species. The results show incomplete concerted evolution and the presence of three ITS ribotypes of lengths 505, 609, 601 bp, hereafter referred to as ITS-S (small), I ITS-M (medium), and ITS-L (large), respectively. Thirty species had only one ribotype (ITS-M), nine species had two ribotypes (different combinations of ITS-L, ITS-M, and ITS-S), and only one species had all three ribotypes. Furthermore, examination of these three ribotypes showed that only ITS-M is putatively functional and ITS-L and ITS-S are pseudogenes. Bayesian analysis strongly supported (100%) two pseudogenes clades but provided weak support for the monophyly of a putative functional clade (ITS-M); moreover, within the "functional" clade, species relationships were uncertain and, in most cases, sequences from the same species failed to group together. The results of the current study suggest that ITS may not be appropriate for phylogeny reconstruction of red oaks due to low levels of interspecific variation and incomplete concerted evolution.

Intragenomic internal transcribed spacer 2 variation in a genus of parasitoid wasps (Hymenoptera: Braconidae): implications for accurate species delimitation and phylogenetic analysis

A recent DNA barcoding study of Australian microgastrines (Hymenoptera: Braconidae) sought to use next-generation sequencing of the cytochrome c oxidase subunit 1 (COI) barcoding gene region, the wingless (WG) gene and the internal transcribed spacer 2 (ITS2) to delimit molecular species in a highly diverse group of parasitic wasps. Large intragenomic distances between ITS2 variants, often larger than the average interspecific variation, caused difficulties in using ITS2 for species delimitation in both threshold and tree-based approaches, and the gene was not included in the reported results of the previous DNA barcoding study. We here report on the intragenomic, and the intra- and interspecies, variation in ITS2in the microgastrine genus Diolcogasterto further investigate the value of ITS2as a marker for species delimitation and phylogenetics of the Microgastrinae. Distinctive intragenomic variant patterns were found in different species of Diolcogaster, with some species possessing a single major variant, and others possessing many divergent variants. Characterizing intragenomic variation of ITS2is critical as it is a widely used marker in hymenopteran phylogenetics and species delimitation, and large intragenomic distances such as those found in this study may obscure phylogenetic signal.

Natural hybridisation among Quercus glabrescens, Q. rugosa and Q. obtusata (Fagaceae): Microsatellites and secondary metabolites markers

Natural hybridisation has significant ecological, genetic and evolutionary consequences altering morphological and chemical characters of individuals. Quercus glabrescens, Q. rugosa and Q. obtusata are white oak species well separated by their morphological characters when they occur in allopatry in Mexican temperate forests. However, in sympatry, individuals with atypical morphology have been observed, suggesting hybridisation events. In this study, we determined, with microsatellites and secondary metabolites, if interspecific gene flow occurs when these three oak species coexist in sympatry. In total, 180 individuals belonging to seven populations [three allopatric (one for each parental species) and four sympatric sites] were analysed. Allopatric populations represent well-defined genetic groups and the sympatric populations showed genetic evidence of hybridisation between Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata. The hybridisation percentage varied between sites and combination of involved species. We registered the presence of unique flavonoid compounds for Q. glabrescens (caffeic acid and flavonol 2), Q. rugosa (flavonol 5) and Q. obtusata (flavonol 1). Three compounds (quercetin rhamnoside, flavonol 3 and alkyl coumarate) were expressed in all taxa. Finally, the hybrid genotypes identified in this study (Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata) showed specific chemical profiles, resulting from a combination of those of their parental species. These results show that hybridisation events between these oak species alter chemical expression of secondary metabolites, creating a mosaic of resources and conditions that provide the substrate for different combinations of foliar-associated species such as herbivores, endophytic fungi or epiphyte plants.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.

Chloroplast variation is incongruent with classification of the Australian bloodwood eucalypts (genus Corymbia, family Myrtaceae)

Previous molecular phylogenetic analyses have resolved the Australian bloodwood eucalypt genus Corymbia (~100 species) as either monophyletic or paraphyletic with respect to Angophora (9-10 species). Here we assess relationships of Corymbia and Angophora using a large dataset of chloroplast DNA sequences (121,016 base pairs; from 90 accessions representing 55 Corymbia and 8 Angophora species, plus 33 accessions of related genera), skimmed from high throughput sequencing of genomic DNA, and compare results with new analyses of nuclear ITS sequences (119 accessions) from previous studies. Maximum likelihood and maximum parsimony analyses of cpDNA resolve well supported trees with most nodes having >95% bootstrap support. These trees strongly reject monophyly of Corymbia, its two subgenera (Corymbia and Blakella), most taxonomic sections (Abbreviatae, Maculatae, Naviculares, Septentrionales), and several species. ITS trees weakly indicate paraphyly of Corymbia (bootstrap support <50% for maximum likelihood, and 71% for parsimony), but are highly incongruent with the cpDNA analyses, in that they support monophyly of both subgenera and some taxonomic sections of Corymbia. The striking incongruence between cpDNA trees and both morphological taxonomy and ITS trees is attributed largely to chloroplast introgression between taxa, because of geographic sharing of chloroplast clades across taxonomic groups. Such introgression has been widely inferred in studies of the related genus Eucalyptus. This is the first report of its likely prevalence in Corymbia and Angophora, but this is consistent with previous morphological inferences of hybridisation between species. Our findings (based on continent-wide sampling) highlight a need for more focussed studies to assess the extent of hybridisation and introgression in the evolutionary history of these genera, and that critical testing of the classification of Corymbia and Angophora requires additional sequence data from nuclear genomes.

Unexpected High Intragenomic Variation in Two of Three Major Pest Thrips Species Does Not Affect Ribosomal Internal Transcribed Spacer 2 (ITS2) Utility for Thrips Identification

The mitochondrial cytochrome oxidase I gene (mtCO1) and the ribosomal internal transcribed spacer 2 region (ITS2) are among the most widely used molecular markers for insect taxonomic characterization. Three economically important species of thrips, Scirtothripsdorsalis, Thripspalmi, and Frankliniellaoccidentalis were selected to examine the extent of intragenomic variation within these two marker regions in the family Thripidae, and determine if this variation would affect the utility of markers in thrips molecular diagnostics. For each species, intragenomic (within individual) variation and intergenomic (among individuals) variation was assessed by cloning and sequencing PCR-amplified copies. Intergenomic variation was generally higher than intragenomic variation except in cases where intergenomic variation was very low, as in mtCO1 from S.dorsalis and F.occidentalis. Intragenomic variation was detected in both markers in all three of the thrips species, however, 2-3 times more intragenomic variation was observed for ITS2 than mtCO1 in both S.dorsalis and T.palmi. Furthermore, levels of intragenomic variation were low for both of the genes in F.occidentalis. In all of the three thrips species, no sex-based clustering of haplotypes was observed in either marker. Unexpected high intragenomic variation in ITS2 for two of three thrips species did not interfere with thrips diagnostics. However, caution should be taken in applying ITS2 to certain studies of S.dorsalis and T.palmi when high levels of intragenomic variation could be problematic or confounding. In such studies, mtCO1 may be a preferable marker. Possible reasons for discrepancies in intragenomic variation among genomic regions are discussed.

rRNA Pseudogenes in Filamentous Ascomycetes as Revealed by Genome Data

The nuclear ribosomal DNA (rDNA) is considered as a paradigm of concerted evolution. Components of the rDNA tandem repeats (45S) are widely used in phylogenetic studies of different organisms and the internal transcribed spacer (ITS) region was recently selected as a fungal DNA bar code. However, rRNA pseudogenes, as one kind of escape from concerted evolution, were reported in a wide range of organisms, especially in plants and animals. Moreover, large numbers of 5S rRNA pseudogenes were identified in several filamentous ascomycetes. To study whether rDNA evolves in a strict concerted manner and test whether rRNA pseudogenes exist in more species of ascomycetes, intragenomic rDNA polymorphisms were analyzed using whole genome sequences. Divergent rDNA paralogs were found to coexist within a single genome in seven filamentous ascomycetes examined. A great number of paralogs were identified as pseudogenes according to the mutation and secondary structure analyses. Phylogenetic analyses of the three rRNA coding regions of the 45S rDNA repeats, i.e., 18S, 5.8S, and 28S, revealed an interspecies clustering pattern of those different rDNA paralogs. The identified rRNA pseudogenic sequences were validated using specific primers designed. Mutation analyses revealed that the repeat-induced point (RIP) mutation was probably responsible for the formation of those rRNA pseudogenes.

ITS non-concerted evolution and rampant hybridization in the legume genus Lespedeza (Fabaceae)

The internal transcribed spacer (ITS) as one part of nuclear ribosomal DNA is one of the most extensively sequenced molecular markers in plant systematics. The ITS repeats generally exhibit high-level within-individual homogeneity, while relatively small-scale polymorphism of ITS copies within individuals has often been reported in literature. Here, we identified large-scale polymorphism of ITS copies within individuals in the legume genus Lespedeza (Fabaceae). Divergent paralogs of ITS sequences, including putative pseudogenes, recombinants, and multiple functional ITS copies were sometimes detected in the same individual. Thirty-seven ITS pseudogenes could be easily detected according to nucleotide changes in conserved 5.8S motives, the significantly lower GC contents in at least one of three regions, and the lost ability of 5.8S rDNA sequence to fold into a conserved secondary structure. The distribution patterns of the putative functional clones were highly different between the traditionally recognized two subgenera, suggesting different rates of concerted evolution in two subgenera which could be attributable to their different extents/frequencies of hybridization, confirmed by our analysis of the single-copy nuclear gene PGK. These findings have significant implications in using ITS marker for reconstructing phylogeny and studying hybridization.

Concerted evolution in the mitochondrial control region of the Amazon small-bodied frog Pseudopaludicola canga (Anura, Leiuperidae)

This study presents evidence of concerted evolution in the mitochondrial control region of the frog Pseudopaludicola canga. Four repeat units of 88 bp (as well as a fifth, incomplete unit) were observed in the 5' domain, with the duplicated segments of the same specimen being more related to one another than to the equivalent regions in other specimens, as a result of concerted evolution. We highlight that drawing conclusions from phylogeographical analysis using the control region containing VNTRs must be interpreted with caution, because it violated a basic assumption of phylogeny, since the regions cannot be treated as independent characters.

Efficient Detection of Novel Nuclear Markers for Brassicaceae by Transcriptome Sequencing

The lack of DNA sequence information for most non-model organisms impairs the design of primers that are universally applicable for the study of molecular polymorphisms in nuclear markers. Next-generation sequencing (NGS) techniques nowadays provide a powerful approach to overcome this limitation. We present a flexible and inexpensive method to identify large numbers of nuclear primer pairs that amplify in most Brassicaceae species. We first obtained and mapped NGS transcriptome sequencing reads from two of the distantly related Brassicaceae species, Cardamine hirsuta and Arabis alpina, onto the Arabidopsis thaliana reference genome, and then identified short conserved sequence motifs among the three species bioinformatically. From these, primer pairs to amplify coding regions (nuclear protein coding loci, NPCL) and exon-primed intron-crossing sequences (EPIC) were developed. We identified 2,334 universally applicable primer pairs, targeting 1,164 genes, which provide a large pool of markers as readily usable genomic resource that will help addressing novel questions in the Brassicaceae family. Testing a subset of the newly designed nuclear primer pairs revealed that a great majority yielded a single amplicon in all of the 30 investigated Brassicaceae taxa. Sequence analysis and phylogenetic reconstruction with a subset of these markers on different levels of phylogenetic divergence in the mustard family were compared with previous studies. The results corroborate the usefulness of the newly developed primer pairs, e.g., for phylogenetic analyses or population genetic studies. Thus, our method provides a cost-effective approach for designing nuclear loci across a broad range of taxa and is compatible with current NGS technologies.

Non-random expression of ribosomal DNA units in a grasshopper showing high intragenomic variation for the ITS2 region

We analyse intragenomic variation of the ITS2 internal transcribed spacer of ribosomal DNA (rDNA) in the grasshopper Eyprepocnemis plorans, by means of tagged PCR 454 amplicon sequencing performed on both genomic DNA (gDNA) and RNA-derived complementary DNA (cDNA), using part of the ITS2 flanking coding regions (5.8S and 28S rDNA) as an internal control for sequencing errors. Six different ITS2 haplotypes (i.e. variants for at least one nucleotide in the complete ITS2 sequence) were found in a single population, one of them (Hap4) being specific to a supernumerary (B) chromosome. The analysis of both gDNA and cDNA from the same individuals provided an estimate of the expression efficiency of the different haplotypes. We found random expression (i.e. about similar recovery in gDNA and cDNA) for three haplotypes (Hap1, Hap2 and Hap5), but significant underexpression for three others (Hap3, Hap4 and Hap6). Hap4 was the most extremely underexpressed and, remarkably, it showed the lowest sequence conservation for the flanking 5.8-28S coding regions in the gDNA reads but the highest conservation (100%) in the cDNA ones, suggesting the preferential expression of mutation-free rDNA units carrying this ITS2 haplotype. These results indicate that the ITS2 region of rDNA is far from complete homogenization in this species, and that the different rDNA units are not expressed at random, with some of them being severely downregulated.

Phylogeny and classification of the East Asian Amitostigma alliance (Orchidaceae: Orchideae) based on six DNA markers

BACKGROUND

Tribe Orchideae dominates the orchid flora of the temperate Northern Hemisphere but its representatives in East Asia had been subject to less intensive phylogenetic study than those in Eurasia and North America. Although this situation was improved recently by the molecular phylogenetic study of Jin et al., comparatively few species were analyzed from the species-rich and taxonomically controversial East Asian Amitostigma alliance. Here, we present a framework nrITS tree of 235 accessions of Orchideae plus an in-depth analysis of 110 representative accessions, encompassing most widely recognized species within the alliance, to elucidate their relationships.

RESULTS

We used parsimony, likelihood and Bayesian approaches to generate trees from data for two nuclear (nrITS, low-copy Xdh) and four chloroplast (matK, psbA-trnH, trnL-F, trnS-trnG) markers. Nuclear and plastid data were analyzed separately due to a few hard incongruences that most likely reflect chloroplast capture. Our results suggest key phylogenetic placements for Sirindhornia and Brachycorythis, and confirm previous assertions that the Amitostigma alliance is monophyletic and sister to the Eurasian plus European clades of subtribe Orchidinae. Seven robust clades are evident within the alliance, but none corresponds precisely with any of the traditional genera; the smaller and more morphologically distinct genera Tsaiorchis, Hemipilia, Neottianthe and Hemipiliopsis are monophyletic but each is nested within a polyphyletic plexus of species attributed to either Ponerorchis or the most plesiomorphic genus, Amitostigma. Two early-divergent clades that escaped analysis by Jin et al. undermine their attempt to circumscribe an expanded monophyletic genus Ponerorchis.

CONCLUSIONS

We provide a new framework on the complex phylogenetic relationships between Amitostigma and other genera traditionally included in its alliance; based on which, we combine the entire Amitostigma alliance into a morphologically and molecularly circumscribed Amitostigma sensu latissimo that also contains seven molecularly circumscribed sections. Our molecular trees imply unusually high levels of morphological homoplasy, but these will need to be quantified via a future group-wide review of the alliance based on living plants if morphology is to be fully integrated into our classification.

Taxonomic and ecological discrimination of Fagaceae species based on internal transcribed spacer polymerase chain reaction-restriction fragment length polymorphism

The internal transcribed spacer (ITS) of ribosomal DNA has been used to confirm taxonomic classifications and define phylogenies in several plant species following sequencing or polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) techniques. In this study, co-dominant ITS PCR-RFLP molecular markers were produced in 30 Fagaceae individuals belonging to the Castanea, Fagus and Quercus genera in order to assess the potential of this technique for taxonomic discrimination and determination of phylogenies. The complete ITS region (ITS1-5.8S rRNA-ITS2) was amplified in most of the Fagaceae individuals as a single fragment of ∼700 bp. The ITS amplified products were digested with nine restriction enzymes, but only four (HaeIII, HpaII, TaqI and Sau96I) produced polymorphic/discriminative patterns. The total expected heterozygosity (HE) was 20.31 % and the gene diversity (I), 32.97 %. The ITS polymorphism was higher within the Quercus genus (85.3 %). The ITS PCR-RFLP markers clustered the Fagaceae species according to genus or infrageneric group (in the case of Quercus sp. individuals). Five oaks did not cluster in line with the adopted infrageneric classification, but three of these were grouped according to their actual ecological distributions. The ITS PCR-RFLP markers indicated their potential for phylogenetic studies since all Fagaceae individuals were discriminated according to genus, and most of the oaks were clustered according to infrageneric group or ecological area.

SSU rDNA divergence in planktonic foraminifera: molecular taxonomy and biogeographic implications

The use of planktonic foraminifera in paleoceanography requires taxonomic consistency and precise assessment of the species biogeography. Yet, ribosomal small subunit (SSUr) DNA analyses have revealed that most of the modern morpho-species of planktonic foraminifera are composed of a complex of several distinct genetic types that may correspond to cryptic or pseudo-cryptic species. These genetic types are usually delimitated using partial sequences located at the 3'end of the SSUrDNA, but typically based on empirical delimitation. Here, we first use patristic genetic distances calculated within and among genetic types of the most common morpho-species to show that intra-type and inter-type genetic distances within morpho-species may significantly overlap, suggesting that genetic types have been sometimes inconsistently defined. We further apply two quantitative and independent methods, ABGD (Automatic Barcode Gap Detection) and GMYC (General Mixed Yule Coalescent) to a dataset of published and newly obtained partial SSU rDNA for a more objective assessment of the species status of these genetic types. Results of these complementary approaches are highly congruent and lead to a molecular taxonomy that ranks 49 genetic types of planktonic foraminifera as genuine (pseudo)cryptic species. Our results advocate for a standardized sequencing procedure allowing homogenous delimitations of (pseudo)cryptic species. On the ground of this revised taxonomic framework, we finally provide an integrative taxonomy synthesizing geographic, ecological and morphological differentiations that can occur among the genuine (pseudo)cryptic species. Due to molecular, environmental or morphological data scarcities, many aspects of our proposed integrative taxonomy are not yet fully resolved. On the other hand, our study opens up the potential for a correct interpretation of environmental sequence datasets.

Mitochondrial and nuclear phylogenetic analysis with Sanger and next-generation sequencing shows that, in Área de Conservación Guanacaste, northwestern Costa Rica, the skipper butterfly named Urbanus belli (family Hesperiidae) comprises three morphologically cryptic species

Background

Skipper butterflies (Hesperiidae) are a relatively well-studied family of Lepidoptera. However, a combination of DNA barcodes, morphology, and natural history data has revealed several cryptic species complexes within them. Here, we investigate three DNA barcode lineages of what has been identified as Urbanus belli (Hesperiidae, Eudaminae) in Área de Conservación Guanacaste (ACG), northwestern Costa Rica.

Results

Although no morphological traits appear to distinguish among the three, congruent nuclear and mitochondrial lineage patterns show that “Urbanus belli” in ACG is a complex of three sympatric species. A single strain of Wolbachia present in two of the three cryptic species indicates that Urbanus segnestami Burns (formerly Urbanus belliDHJ01), Urbanus bernikerni Burns (formerly Urbanus belliDHJ02), and Urbanus ehakernae Burns (formerly Urbanus belliDHJ03) may be biologically separated by Wolbachia, as well as by their genetics. Use of parallel sequencing through 454-pyrosequencing improved the utility of ITS2 as a phylogenetic marker and permitted examination of the intra- and interlineage relationships of ITS2 variants within the species complex. Interlineage, intralineage and intragenomic compensatory base pair changes were discovered in the secondary structure of ITS2.

Conclusion

These findings corroborate the existence of three cryptic species. Our confirmation of a novel cryptic species complex, initially suggested by DNA barcode lineages, argues for using a multi-marker approach coupled with next-generation sequencing for exploration of other suspected species complexes.

Multiple ITS copies reveal extensive hybridization within Rheum (Polygonaceae), a genus that has undergone rapid radiation

Background

During adaptive radiation events, characters can arise multiple times due to parallel evolution, but transfer of traits through hybridization provides an alternative explanation for the same character appearing in apparently non-sister lineages. The signature of hybridization can be detected in incongruence between phylogenies derived from different markers, or from the presence of two divergent versions of a nuclear marker such as ITS within one individual.

Methodology/Principal Findings

In this study, we cloned and sequenced ITS regions for 30 species of the genus Rheum, and compared them with a cpDNA phylogeny. Seven species contained two divergent copies of ITS that resolved in different clades from one another in each case, indicating hybridization events too recent for concerted evolution to have homogenised the ITS sequences. Hybridization was also indicated in at least two further species via incongruence in their position between ITS and cpDNA phylogenies. None of the ITS sequences present in these nine species matched those detected in any other species, which provides tentative evidence against recent introgression as an explanation. Rheum globulosum, previously indicated by cpDNA to represent an independent origin of decumbent habit, is indicated by ITS to be part of clade of decumbent species, which acquired cpDNA of another clade via hybridization. However decumbent and glasshouse morphology are confirmed to have arisen three and two times, respectively.

Conclusions

These findings suggested that hybridization among QTP species of Rheum has been extensive, and that a role of hybridization in diversification of Rheum requires investigation.

Deciduous trees and the application of universal DNA barcodes: a case study on the circumpolar Fraxinus

The utility of DNA barcoding for identifying representative specimens of the circumpolar tree genus Fraxinus (56 species) was investigated. We examined the genetic variability of several loci suggested in chloroplast DNA barcode protocols such as matK, rpoB, rpoC1 and trnH-psbA in a large worldwide sample of Fraxinus species. The chloroplast intergenic spacer rpl32-trnL was further assessed in search for a potentially variable and useful locus. The results of the study suggest that the proposed cpDNA loci, alone or in combination, cannot fully discriminate among species because of the generally low rates of substitution in the chloroplast genome of Fraxinus. The intergenic spacer trnH-psbA was the best performing locus, but genetic distance-based discrimination was moderately successful and only resulted in the separation of the samples at the subgenus level. Use of the BLAST approach was better than the neighbor-joining tree reconstruction method with pairwise Kimura's two-parameter rates of substitution, but allowed for the correct identification of only less than half of the species sampled. Such rates are substantially lower than the success rate required for a standardised barcoding approach. Consequently, the current cpDNA barcodes are inadequate to fully discriminate Fraxinus species. Given that a low rate of substitution is common among the plastid genomes of trees, the use of the plant cpDNA "universal" barcode may not be suitable for the safe identification of tree species below a generic or sectional level. Supplementary barcoding loci of the nuclear genome and alternative solutions are proposed and discussed.

Complex evolution in Arundinarieae (Poaceae: Bambusoideae): incongruence between plastid and nuclear GBSSI gene phylogenies

The monophyly of tribe Arundinarieae (the temperate woody bamboos) has been unequivocally recovered in previous molecular phylogenetic studies. In a recent phylogenetic study, 10 major lineages in Arundinarieae were resolved based on eight non-coding plastid regions, which conflicted significantly with morphological classifications both at the subtribal and generic levels. Nevertheless, relationships among and within the 10 lineages remain unclear. In order to further unravel the evolutionary history of Arundinarieae, we used the nuclear GBSSI gene sequences along with those of eight plastid regions for phylogenetic reconstruction, with an emphasis on Chinese species. The results of the plastid analyses agreed with previous studies, whereas 13 primary clades revealed in the GBSSI phylogeny were better resolved at the generic level than the plastid phylogeny. Our analyses also revealed many inconsistencies between the plastid DNA and the nuclear GBSSI trees. These results implied that the nuclear genome and the plastid genome had different evolutionary trajectories. The patterns of incongruence suggested that lack of informative characters, incomplete lineage sorting, and/or hybridization (introgression) could be the causes. Seven putative hybrid species were hypothesized, four of which are discussed in detail on the basis of topological incongruence, chromosome numbers, morphology, and distribution patterns, and those taxa probably resulted from homoploid hybrid speciation. Overall, our study indicates that the tribe Arundinarieae has undergone a complex evolution.

Early evolutionary colocalization of the nuclear ribosomal 5S and 45S gene families in seed plants: evidence from the living fossil gymnosperm Ginkgo biloba

In seed plants, the colocalization of the 5S loci within the intergenic spacer (IGS) of the nuclear 45S tandem units is restricted to the phylogenetically derived Asteraceae family. However, fluorescent in situ hybridization (FISH) colocalization of both multigene families has also been observed in other unrelated seed plant lineages. Previous work has identified colocalization of 45S and 5S loci in Ginkgo biloba using FISH, but these observations have not been confirmed recently by sequencing a 1.8 kb IGS. In this work, we report the presence of the 45S-5S linkage in G. biloba, suggesting that in seed plants the molecular events leading to the restructuring of the ribosomal loci are much older than estimated previously. We obtained a 6.0 kb IGS fragment showing structural features of functional sequences, and a single copy of the 5S gene was inserted in the same direction of transcription as the ribosomal RNA genes. We also obtained a 1.8 kb IGS that was a truncate variant of the 6.0 kb IGS lacking the 5S gene. Several lines of evidence strongly suggest that the 1.8 kb variants are pseudogenes that are present exclusively on the satellite chromosomes bearing the 45S-5S genes. The presence of ribosomal IGS pseudogenes best reconciles contradictory results concerning the presence or absence of the 45S-5S linkage in Ginkgo. Our finding that both ribosomal gene families have been unified to a single 45S-5S unit in Ginkgo indicates that an accurate reassessment of the organization of rDNA genes in basal seed plants is necessary.

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.

Testing DNA barcoding in closely related groups of Lysimachia L. (Myrsinaceae)

It has been suggested that rbcL and matK are the core barcodes in plants, but they are not powerful enough to distinguish between closely related plant groups. Additional barcodes need to be evaluated to improve the level of discrimination between plant species. Because of their well-studied taxonomy and extreme diversity, we used Chinese Lysimachia (Myrsinaceae) species to test the performance of core barcodes (rbcL and matK) and two additional candidate barcodes (trnH-psbA and the nuclear ribosomal ITS); 97 accessions from four subgenus representing 34 putative Lysimachia species were included in this study. And many closely related species pairs in subgen. Lysimachia were covered to detect their discriminatory power. The inefficiency of rbcL and matK alone or combined in closely related plant groups was validated in this study. TrnH-psbA combined with rbcL + matK did not yet perform well in Lysimachia groups. In contrast, ITS, alone or combined with rbcL and/or matK, revealed high resolving ability in Lysimachia. We support ITS as a supplementary barcode on the basis of core barcode rbcL and matK. Besides, this study also illustrates several mistakes or underlying evolutionary events in Lysimachia detected by DNA barcoding.

The 5S rDNA family evolves through concerted and birth-and-death evolution in fish genomes: an example from freshwater stingrays

Background

Ribosomal 5S genes are well known for the critical role they play in ribosome folding and functionality. These genes are thought to evolve in a concerted fashion, with high rates of homogenization of gene copies. However, the majority of previous analyses regarding the evolutionary process of rDNA repeats were conducted in invertebrates and plants. Studies have also been conducted on vertebrates, but these analyses were usually restricted to the 18S, 5.8S and 28S rRNA genes. The recent identification of divergent 5S rRNA gene paralogs in the genomes of elasmobranches and teleost fishes indicate that the eukaryotic 5S rRNA gene family has a more complex genomic organization than previously thought. The availability of new sequence data from lower vertebrates such as teleosts and elasmobranches enables an enhanced evolutionary characterization of 5S rDNA among vertebrates.

Results

We identified two variant classes of 5S rDNA sequences in the genomes of Potamotrygonidae stingrays, similar to the genomes of other vertebrates. One class of 5S rRNA genes was shared only by elasmobranches. A broad comparative survey among 100 vertebrate species suggests that the 5S rRNA gene variants in fishes originated from rounds of genome duplication. These variants were then maintained or eliminated by birth-and-death mechanisms, under intense purifying selection. Clustered multiple copies of 5S rDNA variants could have arisen due to unequal crossing over mechanisms. Simultaneously, the distinct genome clusters were independently homogenized, resulting in the maintenance of clusters of highly similar repeats through concerted evolution.

Conclusions

We believe that 5S rDNA molecular evolution in fish genomes is driven by a mixed mechanism that integrates birth-and-death and concerted evolution.

A molecular phylogenetic approach to western North America endemic Artemisia and allies (Asteraceae): untangling the sagebrushes

PREMISE OF THE STUDY

Artemisia subgenus Tridentatae plants characterize the North American Intermountain West. These are landscape-dominant constituents of important ecological communities and habitats for endemic wildlife. Together with allied species and genera (Picrothamnus and Sphaeromeria), they make up an intricate series of taxa whose limits are uncertain, likely the result of reticulate evolution. The objectives of this study were to resolve relations among Tridentatae species and their near relatives by delimiting the phylogenetic positions of subgenus Tridentatae species with particular reference to its New World geographic placement and to provide explanations for the relations of allied species and genera with the subgenus with an assessment of their current taxonomic placement.

METHODS

Bayesian inference and maximum parsimony analysis were based on 168 newly generated sequences (including the nuclear ITS and ETS and the plastid trnS(UGA)-trnfM(CAU) and trnS(GCU)-trnC(GCA)) and 338 previously published sequences (ITS and ETS). Genome size by flow cytometry of species from Sphaeromeria was also determined.

KEY RESULTS

The results support an expanded concept and reconfiguration of Tridentatae to accommodate additional endemic North American Artemisia species. The monotypic Picrothamnus and all Sphaeromeria species appear nested within subgenus Tridentatae clade.

CONCLUSIONS

A redefinition of subgenus Tridentatae to include other western North American endemics is supported. We propose a new circumscription of the subgenus and divide it into three sections: Tridentatae, Filifoliae, and Nebulosae. The position of the circumboreal and other North American species suggests that subgenus Artemisia is the ancestral stock for the New World endemics, including those native to South America.

The ITS1-5.8S-ITS2 sequence region in the Musaceae: structure, diversity and use in molecular phylogeny

Genes coding for 45S ribosomal RNA are organized in tandem arrays of up to several thousand copies and contain 18S, 5.8S and 26S rRNA units separated by internal transcribed spacers ITS1 and ITS2. While the rRNA units are evolutionary conserved, ITS show high level of interspecific divergence and have been used frequently in genetic diversity and phylogenetic studies. In this work we report on the structure and diversity of the ITS region in 87 representatives of the family Musaceae. We provide the first detailed information on ITS sequence diversity in the genus Musa and describe the presence of more than one type of ITS sequence within individual species. Both Sanger sequencing of amplified ITS regions and whole genome 454 sequencing lead to similar phylogenetic inferences. We show that it is necessary to identify putative pseudogenic ITS sequences, which may have negative effect on phylogenetic reconstruction at lower taxonomic levels. Phylogenetic reconstruction based on ITS sequence showed that the genus Musa is divided into two distinct clades – Callimusa and Australimusa and Eumusa and Rhodochlamys. Most of the intraspecific banana hybrids analyzed contain conserved parental ITS sequences, indicating incomplete concerted evolution of rDNA loci. Independent evolution of parental rDNA in hybrids enables determination of genomic constitution of hybrids using ITS. The observation of only one type of ITS sequence in some of the presumed interspecific hybrid clones warrants further study to confirm their hybrid origin and to unravel processes leading to evolution of their genomes.

Barcoding diatoms: exploring alternatives to COI-5P

Diatoms are a diverse lineage with species that can be difficult to identify or cryptic, but DNA barcoding, a molecular technique, can assist identification and facilitate studies of speciation and biogeography. The most common region used for DNA barcoding, COI-5P, can distinguish diatom species, but has not displayed universality (i.e., successful PCR amplification from diverse taxa). Therefore, we have assessed the following alternative markers: ∼1400bp of rbcL; 748bp at the 3' end of rbcL (rbcL-3P); LSU D2/D3 and UPA. Sellaphora isolates were used to determine each marker's ability to discriminate among closely related species and culture collection material was utilized to explore further marker universality. All of the alternative markers investigated have greater universality than COI-5P. Both full and partial (3P) rbcL regions had the power to discriminate between all species, but rbcL-3P can be sequenced more easily. LSU D2/D3 could distinguish between all but the most closely related species (96%), whereas UPA only distinguished 20% of species. Our observations suggest that rbcL-3P should be used as the primary marker for diatom barcoding, while LSU D2/D3 should be sequenced as a secondary marker to facilitate environmental surveys.

Biogeography and divergence times of genus Macroptilium (Leguminosae)

BACKGROUND AND AIMS

Macroptilium is a herbaceous legume genus with 18 currently accepted species, seven of them with economic importance due to their use as forage, green fertilizer and in medicine. The genus is strictly American, with an unknown biogeographic history. The aim of this study was to infer a biogeographic pattern of Macroptilium and to estimate its divergence times, using sequences from the nuclear ribosomal DNA internal transcribed spacers.

METHODOLOGY

To study the historical biogeography of Macroptilium, two approaches were used: area optimization on a previously obtained phylogeny and a dispersal-vicariance analysis. Divergence times were calculated by Bayesian methods.

PRINCIPAL RESULTS

The analyses revealed that Macroptilium has its origin in the middle Pliocene, with an estimated age that ranges from 2.9 to 4 million years. The biogeographic analyses placed its origin in South America, specifically on the Chaquean sub-region, where most of the cladogenetic events of the genus took place.

CONCLUSIONS

Macroptilium constitutes a further example of the geographic pattern displayed by numerous Neotropical taxa that moved north from South America to dominate the Central American lowlands after the land connection across the Isthmus of Panama was established.

The evolution pattern of rDNA ITS in Avena and phylogenetic relationship of the Avena species (Poaceae: Aveneae)

Ribosomal ITS sequences are commonly used for phylogenetic reconstruction because they are included in rDNA repeats, and these repeats often undergo rapid concerted evolution within and between arrays. Therefore, the rDNA ITS copies appear to be virtually identical and can sometimes be treated as a single gene. In this paper we examined ITS polymorphism within and among 13 diploid (A and C genomes), seven tetraploid (AB, AC and CC genomes) and four hexaploid (ACD genome) to infer the extent and direction of concerted evolution, and to reveal the phylogenetic and genome relationship among species of Avena. A total of 170 clones of the ITS1-5.8S-ITS2 fragment were sequenced to carry out haplotype and phylogenetic analysis. In addition, 111 Avena ITS sequences retrieved from GenBank were combined with 170 clones to construct a phylogeny and a network. We demonstrate the major divergence between the A and C genomes whereas the distinction among the A and B/D genomes was generally not possible. High affinity among the A(d) genome species A. damascena and the ACD genome species A. fatua was found, whereas the rest of the ACD genome hexaploids and the AACC tetraploids were highly affiliated with the A(l) genome diploid A. longiglumis. One of the AACC species A. murphyi showed the closest relationship with most of the hexaploid species. Both C(v) and C(p) genome species have been proposed as paternal donors of the C-genome carrying polyploids. Incomplete concerted evolution is responsible for the observed differences among different clones of a single Avena individual. The elimination of C-genome rRNA sequences and the resulting evolutionary inference of hexaploid species are discussed.