Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III

Identification of Amazonian trees with DNA barcodes

BACKGROUND

Large-scale plant diversity inventories are critical to develop informed conservation strategies. However, the workload required for classic taxonomic surveys remains high and is particularly problematic for megadiverse tropical forests.

METHODOLOGY/PRINCIPAL FINDINGS

Based on a comprehensive census of all trees in two hectares of a tropical forest in French Guiana, we examined whether plant DNA barcoding could contribute to increasing the quality and the pace of tropical plant biodiversity surveys. Of the eight plant DNA markers we tested (rbcLa, rpoC1, rpoB, matK, ycf5, trnL, psbA-trnH, ITS), matK and ITS had a low rate of sequencing success. More critically, none of the plastid markers achieved a rate of correct plant identification greater than 70%, either alone or combined. The performance of all barcoding markers was noticeably low in few species-rich clades, such as the Laureae, and the Sapotaceae. A field test of the approach enabled us to detect 130 molecular operational taxonomic units in a sample of 252 juvenile trees. Including molecular markers increased the identification rate of juveniles from 72% (morphology alone) to 96% (morphology and molecular) of the individuals assigned to a known tree taxon.

CONCLUSION/SIGNIFICANCE

We conclude that while DNA barcoding is an invaluable tool for detecting errors in identifications and for identifying plants at juvenile stages, its limited ability to identify collections will constrain the practical implementation of DNA-based tropical plant biodiversity programs.

The evolution of antiherbivore defenses and their contribution to species coexistence in the tropical tree genus Inga

Plants and their herbivores constitute more than half of the organisms in tropical forests. Therefore, a better understanding of the evolution of plant defenses against their herbivores may be central for our understanding of tropical biodiversity. Here, we address the evolution of antiherbivore defenses and their possible contribution to coexistence in the Neotropical tree genus Inga (Fabaceae). Inga has >300 species, has radiated recently, and is frequently one of the most diverse and abundant genera at a given site. For 37 species from Panama and Peru we characterized developmental, ant, and chemical defenses against herbivores. We found extensive variation in defenses, but little evidence of phylogenetic signal. Furthermore, in a multivariate analysis, developmental, ant, and chemical defenses varied independently (were orthogonal) and appear to have evolved independently of each other. Our results are consistent with strong selection for divergent defensive traits, presumably mediated by herbivores. In an analysis of community assembly, we found that Inga species co-occurring as neighbors are more different in antiherbivore defenses than random, suggesting that possessing a rare defense phenotype increases fitness. These results imply that interactions with herbivores may be an important axis of niche differentiation that permits the coexistence of many species of Inga within a single site. Interactions between plants and their herbivores likely play a key role in the generation and maintenance of the conspicuously high plant diversity in the tropics.

Phylogeography of Camassia quamash in western North America: postglacial colonization and transport by indigenous peoples

Recent human activities have spread numerous plant species across the globe, yet it is unclear to what degree historical human activities influenced plant dispersal. In western North America, Camassia quamash was one of the most important food plants for indigenous peoples, who transported its propagules either intentionally or accidentally. We investigated how human and natural dispersal might have contributed to the current pattern of spatial genetic structure in C. quamash by performing phylogeographical surveys at two geographical scales. We sequenced two noncoding regions of chloroplast deoxyribonucleic acid (DNA) in 226 individuals from 53 populations of C. quamash as well as 126 individuals from 21 populations of the non-food plant Zigadenus venenosus. Contrary to the expectation of anthropogenic transport, C. quamash populations did not exhibit weaker genetic structure than Z. venenosus populations. We also failed to find convincing evidence for signatures of transport. Instead, our data showed strong effects of past glaciation and geographical barriers of the mountains in the Cascade Range, Olympic Peninsula and Vancouver Island. West of the Cascades, the species appears to have largely migrated northward from a southern refugium after deglaciation, whereas few populations having a highly divergent haplotype might have survived in southwestern Washington. Our data suggest that despite substantial ethnobotanical evidence for anthropogenic transport, the current pattern of genetic structure of C. quamash does not show any detectable signatures of transport by indigenous peoples and is better understood as the result of natural dispersal processes.

The evolution of self-compatible and self-incompatible populations in a hermaphroditic perennial, Trillium camschatcense (Melanthiaceae)

The evolution of selfing from outcrossing ancestors is known to have occurred repeatedly in angiosperms. Theoretical studies have argued that the transition from outcrossing to selfing is accomplished more easily than the reverse case, and phylogenetic analyses involving self-compatible (SC) and self-incompatible (SI) species has basically supported this assumption. The evolutionary direction of self-compatibility and self-incompatibility was examined in Trillium camschatcense, which contains geographically widespread SC populations, and restricted SI populations. Ecological surveys have revealed that the SC populations were suitable for outcrossing, and selfing in these populations did not confer any fitness advantage. Since reproductive fitness indicates the possibility of an evolutionary shift from self-compatibility to self-incompatibility, the phylogenetic relationships of SI and SC populations of T. camschatcense were investigated based on cpDNA variations and nuclear DNA microsatellite polymorphisms. Although phylogenetic analyses did not provide credible evidence to determine evolutionary direction, the SI populations turned out to be monophyletic with extremely low genetic differentiation. Based on these results, we proposed two possible scenarios for the evolutionary backgrounds of SI and SC populations in T. camschatcense. The plausibility of each scenario was evaluated based on the reproductive and geographical features of the mating systems.

The utility of nuclear conserved ortholog set II (COSII) genomic regions for species-level phylogenetic inference in Lycium (Solanaceae)

The identification of genomic regions with sufficient variation to elucidate fine-scale relationships among closely related species is a major goal of phylogenetic systematics. However, the accumulation of such multi-locus data sets brings its own challenges, given that gene trees do not necessarily represent the true species tree. Using genomic tools developed for Solanum (Solanaceae), we have evaluated the utility of nuclear conserved ortholog set II (COSII) regions for phylogenetic inference in tribe Lycieae (Solanaceae). Five COSII regions, with intronic contents ranging from 68% to 91%, were sequenced in 10 species. Their phylogenetic utility was assessed and compared with data from more commonly used nuclear (GBSSI, nrITS) and cpDNA spacer data. We compared the effectiveness of a traditional total evidence concatenation approach versus the recently developed Bayesian estimation of species trees (BEST) method to infer species trees given multiple independent gene trees. All of the sampled COSII regions had high numbers of parsimony-informative (PI) characters, and two of the COSII regions had more PI characters than the GBSSI, ITS, and cpDNA spacer data sets combined. COSII markers are a promising new tool for phylogenetic inference in Solanaceae, and should be explored in related groups. Both the concatenation and BEST approaches yielded similar topologies; however, when multiple individuals with polyphyletic alleles were included, BEST was clearly the more robust approach for inferring species trees in the presence of gene tree incongruence.

Pleistocene refugia and polytopic replacement of diploids by tetraploids in the Patagonian and Subantarctic plant Hypochaeris incana (Asteraceae, Cichorieae)

We report the phylogeographic pattern of the Patagonian and Subantarctic plant Hypochaeris incana endemic to southeastern South America. We applied amplified fragment length polymorphism (AFLP) and chloroplast DNA (cpDNA) analysis to 28 and 32 populations, respectively, throughout its distributional range and assessed ploidy levels using flow cytometry. While cpDNA data suggest repeated or simultaneous parallel colonization of Patagonia and Tierra del Fuego by several haplotypes and/or hybridization, AFLPs reveal three clusters corresponding to geographic regions. The central and northern Patagonian clusters (approximately 38-51 degrees S), which are closer to the outgroup, contain mainly tetraploid, isolated and highly differentiated populations with low genetic diversity. To the contrary, the southern Patagonian and Fuegian cluster (approximately 51-55 degrees S) contains mainly diploid populations with high genetic diversity and connected by high levels of gene flow. The data suggest that H. incana originated at the diploid level in central or northern Patagonia, from where it migrated south. All three areas, northern, central and southern, have similar levels of rare and private AFLP bands, suggesting that all three served as refugia for H. incana during glacial times. In southern Patagonia and Tierra del Fuego, the species seems to have expanded its populational system in postglacial times, when the climate became warmer and more humid. In central and northern Patagonia, the populations seem to have become restricted to favourable sites with increasing temperature and decreasing moisture and there was a parallel replacement of diploids by tetraploids in local populations.

Identification of highly variable chloroplast sequences and development of cpDNA-based molecular markers that distinguish four cytoplasm types in radish (Raphanus sativus L.)

Four types of cytoplasms (Ogura, DCGMS, DBRMF1, and DBRMF2) were identified in the previous studies using molecular markers based on mitochondrial genome variations in radish (Raphanus sativus L.). However, mtDNA markers have limitations in obtaining clear results due to complexity of radish mitochondrial genomes. To improve fidelity, molecular markers based on variation of chloroplast genome sequences were developed in this study. We searched for the sequence variations of chloroplast genome among the four cytoplasm types in 11 noncoding intergenic regions of ~8.7 kb. Highly variable intergenic regions between trnK and rps16 were identified, and a couple of 4-34 bp indels were used to develop a simple PCR-based marker that distinguished the four cytoplasm types based on the PCR product length polymorphism. Two additional cpDNA markers were developed by using a single nucleotide polymorphism and 17-bp insertion. Analysis of 90 accessions using both mtDNA and cpDNA markers showed the perfect match of results of both the markers, suggesting strict co-transmission of mitochondria and chloroplast in radish. Phylogenetic trees showed that two male-sterility inducing cytoplasms, Ogura and DCGMS, were closely related to DBRMF1 and DBRMF2, respectively. Analysis of 120 radish germplasms introduced from diverse countries showed that the frequency of male-sterility inducing mitotypes of Ogura and DCGMS was very low, and DCGMS was predominately detected in eastern European countries. Majority of accessions from Europe and Asia were shown to contain DBRMF2 and DBRMF1 mitotypes, respectively.

Chloroplast simple sequence repeats (cpSSRs): technical resources and recommendations for expanding cpSSR discovery and applications to a wide array of plant species

Chloroplast microsatellites, or simple sequence repeats (cpSSRs), are typically mononucleotide tandem repeats. When located in the noncoding regions of the chloroplast genome (cpDNA), they commonly show intraspecific variation in repeat number. Despite the growing number of studies applying cpSSRs, studies of economically important plants and their relatives remain over-represented. Thus, the potential of cpSSRs to offer unique insights into ecological and evolutionary processes in wild plant species has yet to be fully realized. This review provides an overview of the technical resources available to aid cpSSR discovery including a list of cpSSR primer sets available and cpDNA sequencing resources. Our updated analysis of 99 whole chloroplast genomes downloaded from GenBank confirms that potentially variable cpSSRs are abundant in the noncoding cpDNA of plants. Overall variation in the frequency of cpSSRs was extreme, ranging from one to 700 per genome (median = 93), while in 81 vascular plants, between 35 and 160 cpSSRs were detected per genome (median = 86). We offer five recommendations to aid wider development and application of cpSSRs: (i) When genus-specific cpSSR primers are available, cross-species amplification can often be fruitful. (ii) While potentially useful, universal cpSSR primers at best provide access to only a small number of variable markers. (iii) De novo sequencing of noncoding cpDNA is the most effective and efficient way to develop cpSSR markers in wild species. (iv) DNA sequencing of cpSSR alleles is essential, given the complex nature of the genetic variation associated with hypervariable cpDNA regions. (v) The reliability of cpSSR length based genetic assays need to be validated in all studies.

A new set of universal de novo sequencing primers for extensive coverage of noncoding chloroplast DNA: new opportunities for phylogenetic studies and cpSSR discovery

We present a new set of universal de novo sequencing primers targeting noncoding chloroplast DNA. The set of 107 polymerase chain reaction (PCR) primers span approximately 86% of the noncoding nucleotides in the large single copy region of Nicotiana tabacum, Oryza sativa and the orchid Phalaenopsis aphrodite. PCR tests confirmed the primers are effective in a wide range of monocots and dicots. More than 19.5 kb of cpDNA sequence was obtained across representative orchid genera with up to 82 chloroplast simple sequence repeats (cpSSRs) detected per genus. This primer set will facilitate both phylogenetic studies and rapid discovery of cpSSRs for plants, such as orchids, where there are limited genomic resources.

Multiplex sequencing of plant chloroplast genomes using Solexa sequencing-by-synthesis technology

Organellar DNA sequences are widely used in evolutionary and population genetic studies, however, the conservative nature of chloroplast gene and genome evolution often limits phylogenetic resolution and statistical power. To gain maximal access to the historical record contained within chloroplast genomes, we have adapted multiplex sequencing-by-synthesis (MSBS) to simultaneously sequence multiple genomes using the Illumina Genome Analyzer. We PCR-amplified approximately 120 kb plastomes from eight species (seven Pinus, one Picea) in 35 reactions. Pooled products were ligated to modified adapters that included 3 bp indexing tags and samples were multiplexed at four genomes per lane. Tagged microreads were assembled by de novo and reference-guided assembly methods, using previously published Pinus plastomes as surrogate references. Assemblies for these eight genomes are estimated at 88-94% complete, with an average sequence depth of 55x to 186x. Mononucleotide repeats interrupt contig assembly with increasing repeat length, and we estimate that the limit for their assembly is 16 bp. Comparisons to 37 kb of Sanger sequence show a validated error rate of 0.056%, and conspicuous errors are evident from the assembly process. This efficient sequencing approach yields high-quality draft genomes and should have immediate applicability to genomes with comparable complexity.

The D4 set: primers that target highly variable intron loops in plant chloroplast genomes

Chloroplast group II introns offer high-quality, rapidly evolving single-copy loci for comparative sequence analysis. These introns feature diagnostic secondary structures with loops that are among the least evolutionarily constrained sequence in plastomes. We exploited these structures to develop universal primers that amplify and sequence the large Domain IV (D4) loop in several angiosperm introns. With a single sequence read, we recover 300-600 nucleotides of highly variable sequence across angiosperms, with rates of change that are equal to or higher than many of the best known intergenic spacers in plant chloroplast genomes.

Evolution and phylogenetic utility of the PHOT gene duplicates in the Verbena complex (Verbenaceae): dramatic intron size variation and footprint of ancestral recombination

A well-resolved species level phylogeny is critically important in studying organismal evolution (e.g., hybridization, polyploidization, adaptive speciation). Lack of appropriate molecular markers that give sufficient resolution to gene trees is one of the major impediments to inferring species level phylogenies. In addition, sampling multiple independent loci is essential to overcome the lineage sorting problem. The availability of nuclear loci has often been a limiting factor in plant species-level phylogenetic studies. Here the two PHOT loci were developed as new sources of nuclear gene trees. The PHOT1 and PHOT2 gene trees of the Verbena complex (Verbenaceae) are well resolved and have good clade support. These gene trees are consistent with each other and previously generated chloroplast and nuclear waxy gene trees in most of the phylogenetic backbone as well as some terminal relationships, but are incongruent in some other relationships. Locus-specific primers were optimized for amplifying and sequencing these two loci in all Lamiales. Comparing intron size in the context of the gene trees shows dramatic variation within the Verbena complex, particularly at the PHOT1 locus. These variations are largely caused by invasions of short transposable elements and frequent long deletions and insertions of unknown causes. In addition, inspection of DNA sequences and phylogenetic analyses unmask a clear footprint of ancestral recombination in one species.

Intraspecific variation in Viola suavis in Europe: parallel evolution of white-flowered morphotypes

BACKGROUND AND AIMS

Viola species are commonly grown for their ornamental flowers, but their evolutionary history and taxonomy are often complicated and have been poorly explored so far. This is a study of the polymorphic, typically blue-flowered species Viola suavis, concentrating on the white-flowered populations of uncertain taxonomic assignment that occur in Spain and central and south-eastern Europe. The aim was to resolve their origin and taxonomic status and to study the intraspecific structure and (post)glacial history of this species.

METHODS

Viola suavis and five close relatives were sampled from multiple locations and subjected to molecular (AFLP, sequencing of nrDNA ITS) and morphometric analyses. Data on ploidy level and pollen fertility were also obtained, to address an assumed hybrid origin of the white-flowered populations.

KEY RESULTS

In V. suavis a strong intraspecific genetic split into two groups was observed, indicating that there has been a long-term isolation and survival in distinct glacial refugia. The white-flowered populations could be placed within the variation range of this species, and it is clear that they evolved independently in two distant areas. Their parallel evolution is supported by both morphological and genetic differentiation. The strongly reduced genetic variation and absence of unique AFLP fragments suggest their derived status and origin from the typical, blue-flowered populations.

CONCLUSIONS

These results suggest that intraspecific variation in V. suavis has been largely shaped by population isolations during the last glaciation and subsequent recolonizations, although cultivation and vegetative spread by humans have affected the present picture as well.

Phylogeny of Nolana (Solanaceae) of the Atacama and Peruvian deserts inferred from sequences of four plastid markers and the nuclear LEAFY second intron

The phylogeny of Nolana (Solanaceae), a genus primarily distributed in the coastal Atacama and Peruvian deserts with a few species in the Andes and one species endemic to the Galápagos Islands, was reconstructed using sequences of four plastid regions (ndhF, psbA-trnH, rps16-trnK and trnC-psbM) and the nuclear LEAFY second intron. The monophyly of Nolana was strongly supported by all molecular data. The LEAFY data suggested that the Chilean species, including Nolana sessiliflora, the N. acuminata group and at least some members of the Alona group, are basally diverged, supporting the Chilean origin of the genus. Three well-supported clades in the LEAFY tree were corroborated by the SINE (short interspersed elements) or SINE-like insertions. Taxa from Peru are grouped roughly into two clades. Nolana galapagensis from the Galápagos Island is most likely to have derived from a Peruvian ancestor. The monophyly of the morphologically well-diagnosed Nolana acuminata group (N. acuminata, N. baccata, N. paradoxa, N. parviflora, N. pterocarpa, N. rupicola and N. elegans) was supported by both plastid and LEAFY data. Incongruence between the plastid and the LEAFY data was detected concerning primarily the positions of N. sessiliflora, N. galapagensis, taxa of the Alona group and the two Peruvian clades. Such incongruence may be due to reticulate evolution or in some cases lineage sorting of plastid DNA. Incongruence between our previous GBSSI trees and the plastid-LEAFY trees was also detected concerning two well-supported major clades in the GBSSI tree. Duplication of the GBSSI gene may have contributed to this incongruence.

Multiple multilocus DNA barcodes from the plastid genome discriminate plant species equally well

A universal barcode system for land plants would be a valuable resource, with potential utility in fields as diverse as ecology, floristics, law enforcement and industry. However, the application of plant barcoding has been constrained by a lack of consensus regarding the most variable and technically practical DNA region(s). We compared eight candidate plant barcoding regions from the plastome and one from the mitochondrial genome for how well they discriminated the monophyly of 92 species in 32 diverse genera of land plants (N = 251 samples). The plastid markers comprise portions of five coding (rpoB, rpoC1, rbcL, matK and 23S rDNA) and three non-coding (trnH-psbA, atpF-atpH, and psbK-psbI) loci. Our survey included several taxonomically complex groups, and in all cases we examined multiple populations and species. The regions differed in their ability to discriminate species, and in ease of retrieval, in terms of amplification and sequencing success. Single locus resolution ranged from 7% (23S rDNA) to 59% (trnH-psbA) of species with well-supported monophyly. Sequence recovery rates were related primarily to amplification success (85-100% for plastid loci), with matK requiring the greatest effort to achieve reasonable recovery (88% using 10 primer pairs). Several loci (matK, psbK-psbI, trnH-psbA) were problematic for generating fully bidirectional sequences. Setting aside technical issues related to amplification and sequencing, combining the more variable plastid markers provided clear benefits for resolving species, although with diminishing returns, as all combinations assessed using four to seven regions had only marginally different success rates (69-71%; values that were approached by several two- and three-region combinations). This performance plateau may indicate fundamental upper limits on the precision of species discrimination that is possible with DNA barcoding systems that include moderate numbers of plastid markers. Resolution to the contentious debate on plant barcoding should therefore involve increased attention to practical issues related to the ease of sequence recovery, global alignability, and marker redundancy in multilocus plant DNA barcoding systems.

Phylogenetic utility of 141 low-copy nuclear regions in taxa at different taxonomic levels in two distantly related families of rosids

Angiosperm systematics has progressed to the point where it is now expected that multiple, independent markers be used in phylogenetic studies. Universal primers for amplifying informative regions of the chloroplast genome are readily available, but in the faster-evolving nuclear genome it is challenging to discover priming sites that are conserved across distantly related taxa. With goals including the identification of informative markers in rosids, and perhaps other angiosperms, we screened 141 nuclear primer combinations for phylogenetic utility in two distinct groups of rosids at different taxonomic levels-Psiguria (Cucurbitaceae) and Geraniaceae. We discovered three phylogenetically informative regions in Psiguria and two in Geraniaceae, but none that were useful in both groups. Extending beyond rosids, we combined our findings with those of another recent effort testing these primer pairs in Asteraceae, Brassicaceae, and Orchidaceae. From this comparison, we identified 32 primer combinations that amplified regions in representative species of at least two of the five distantly related angiosperm families, giving some prior indication about phylogenetic usefulness of these markers in other flowering plants. This reduced set of primer pairs for amplifying low-copy nuclear markers along with a recommended experimental strategy provide a framework for identifying phylogenetically informative regions in angiosperms.

Phylogeny of marattioid ferns (Marattiaceae): inferring a root in the absence of a closely related outgroup

Closely related outgroups are optimal for rooting phylogenetic trees; however, such ideal outgroups are not always available. A phylogeny of the marattioid ferns (Marattiaceae), an ancient lineage with no close relatives, was reconstructed using nucleotide sequences of multiple chloroplast regions (rps4 + rps4-trnS spacer, trnS-trnG spacer + trnG intron, rbcL, atpB), from 88 collections, selected to cover the broadest possible range of morphologies and geographic distributions within the extant taxa. Because marattioid ferns are phylogenetically isolated from other lineages, and internal branches are relatively short, rooting was problematic. Root placement was strongly affected by long-branch attraction under maximum parsimony and by model choice under maximum likelihood. A multifaceted approach to rooting was employed to isolate the sources of bias and produce a consensus root position. In a statistical comparison of all possible root positions with three different outgroups, most root positions were not significantly less optimal than the maximum likelihood root position, including the consensus root position. This phylogeny has several important taxonomic implications for marattioid ferns: Marattia in the broad sense is paraphyletic; the Hawaiian endemic Marattia douglasii is most closely related to tropical American taxa; and Angiopteris is monophyletic only if Archangiopteris and Macroglossum are included.

Reticulate or tree-like chloroplast DNA evolution in Sileneae (Caryophyllaceae)?

Despite sampling of up to 25kb of chloroplast DNA sequence from 24 species in Sileneae a number of nodes in the phylogeny remain poorly supported and it is not expected that additional sequence sampling will converge to a reliable phylogenetic hypothesis in these parts of the tree. The main reason for this is probably a combination of rapid radiation and substitution rate heterogeneity. Poor resolution among closely related species are often explained by low levels of variation in chloroplast data, but the problem with our data appear to be high levels of homoplasy. Tree-like cpDNA evolution cannot be rejected, but apparent incongruent patterns between different regions are evaluated with the possibility of ancient interspecific chloroplast recombination as explanatory model. However, several major phylogenetic relationships, previously not recognized, are confidently resolved, e.g. the grouping of the two SW Anatolian taxa S. cryptoneura and S. sordida strongly disagrees with previous studies on nuclear DNA sequence data, and indicate a possible case of homoploid hybrid origin. The closely related S. atocioides and S. aegyptiaca form a sister group to Lychnis and the rest of Silene, thus suggesting that Silene may be paraphyletic, despite recent revisions based on molecular data.

Complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera, and comparative analyses with other grass genomes

Comparisons of complete chloroplast genome sequences of Hordeum vulgare, Sorghum bicolor and Agrostis stolonifera to six published grass chloroplast genomes reveal that gene content and order are similar but two microstructural changes have occurred. First, the expansion of the IR at the SSC/IRa boundary that duplicates a portion of the 5' end of ndhH is restricted to the three genera of the subfamily Pooideae (Agrostis, Hordeum and Triticum). Second, a 6 bp deletion in ndhK is shared by Agrostis, Hordeum, Oryza and Triticum, and this event supports the sister relationship between the subfamilies Erhartoideae and Pooideae. Repeat analysis identified 19-37 direct and inverted repeats 30 bp or longer with a sequence identity of at least 90%. Seventeen of the 26 shared repeats are found in all the grass chloroplast genomes examined and are located in the same genes or intergenic spacer (IGS) regions. Examination of simple sequence repeats (SSRs) identified 16-21 potential polymorphic SSRs. Five IGS regions have 100% sequence identity among Zea mays, Saccharum officinarum and Sorghum bicolor, whereas no spacer regions were identical among Oryza sativa, Triticum aestivum, H. vulgare and A. stolonifera despite their close phylogenetic relationship. Alignment of EST sequences and DNA coding sequences identified six C-U conversions in both Sorghum bicolor and H. vulgare but only one in A. stolonifera. Phylogenetic trees based on DNA sequences of 61 protein-coding genes of 38 taxa using both maximum parsimony and likelihood methods provide moderate support for a sister relationship between the subfamilies Erhartoideae and Pooideae.