Inferring phylogeny at low taxonomic levels: utility of rapidly evolving cpDNA and nuclear ITS loci

Plant molecular systematic studies of closely related taxa have relied heavily on sequence data from nuclear ITS and cpDNA. Positive attributes of using ITS sequence data include the rapid rate of evolution compared to most plastid loci and availability of universal primers for amplification and sequencing. On the other hand, ITS sequence data may not adequately track organismal phylogeny if concerted evolution and high rDNA array copy number do not permit identification of orthologous copies. Shaw et al. (American Journal of Botany 92: 142-166) recently identified nine plastid regions that appear to provide more potentially informative characters than many other plastid loci. In the present study, sequences of these loci and ITS were obtained for six taxonomic groups in which phylogenetic relationships have been difficult to establish using other data. The relative utility of these regions was compared by assessing the number of parsimony informative characters, character congruence, resolution of inferred trees, clade support, and accuracy. No single locus emerged as the best in all lineages for any of these measures of utility. Results further indicated that in preliminary studies, sampling strategy should include at least four exemplar taxa. The importance of sampling data from independent distributions is also discussed.     

Phylogeny of basal eudicots: Insights from non-coding and rapidly evolving DNA

Sequence data of the trnL group I intron, the petD group II intron, the trnL-F and petB-D spacers, and the rapidly evolving matK gene were analysed from all families of the basal eudicot grade and from representatives of 19 core eudicot orders. The dataset comprised 5654 positions of aligned sequence plus a matrix of 1087 binary indel characters. Mutational hotspots correspond in number and extension to hotspots already known from basal angiosperms and, with respect to secondary structure, are generally located in terminal parts of stem-loop regions. Parsimony, Bayesian, and likelihood analyses depict Ranunculales as sister to all remaining eudicots with maximum support. The branching order in the basal eudicot grade is further resolved as Sabiales, Proteales, Trochodendrales, and Buxales. Nearly all of the backbone nodes gain high confidence, except for the node showing Proteales diverging before Trochodendrales, which is only moderately supported (83% JK). In Ranunculales, the woody Eupteleaceae are first-branching, with Papaveraceae plus Fumariaceae coming next. Within Proteales, Nelumbo is clearly resolved as sister to a Platanaceae–Proteaceae clade. Gunnerales are found as the first branch in core eudicots, with maximum support in the combined analysis. This node is also resolved with matK alone, but unsupported. It appears that the combined analysis of sequence data from rapidly evolving and non-coding genomic regions leads to significantly improved statistical support values in comparison to earlier studies of basal eudicots using multiple conserved genes.See also Electronic Supplement at doi:10.1016/j.ode.2006.08.001     

Development of DNA barcodes for selected Acacia species by using rbcL and matK DNA markers

Acacia species are very important tree species in tropical and subtropical countries of the World for their economic and medicinal benefits. Precise identification of Acacia is very important to distinguish the invasive species from rare species however, it is difficult to differentiate Acacia species based on morphological charcters. In addition, precise identification is also important for wood charcterization in the forest industry as these species are declining due to illegal logging and deforestation. To overcome thsese limitations of morphological identification, DNA barcoding is being used as an efficient and quick approach for precise identification of tree species. In this study, we selected two chloroplast and plastid base DNA markers (rbcL and matK) for the identification of five selected tree species of Acacia (A. albida, A. ampliceps, A. catechu, A. coriacea and A. tortilis). The genomic DNA of the selected Acacia species was extracted, amplified through PCR using specific primers and subsequently sequenced through Sanger sequencing. In matK DNA marker the average AT nucleotide contents were higher (59.46%) and GC contents were lower (40.44%) as compared to the AT (55.40%) and GC content (44.54%) in rbcL marker. The means genetic distance K2P between the Acacia species was higher in matK (0.704%) as compared to rbcL (0.230%). All Acacia species could be identified based on unique SNPs profile. Based on SNP data profiles, DNA sequence based scannable QR codes were developed for accurate identification of Acacia species. The phylogenetic analysis based on both markers (rbcL and matK) showed that both A. coriacea and A. tortilis were closely related with each other and clustered in the same group while other two species A. albida and A. catechu were grouped together. The specie A. ampliceps remained ungrouped distantly, compared with other four species. These finding highlights the potential of DNA barcoding for efficient and reproducible identification of Acacia species.     

Phylogenetic utility of ycf1 in orchids: a plastid gene more variable than matK

Plastid DNA sequences have been widely used by systematists for reconstructing plant phylogenies. The utility of any DNA region for phylogenetic analysis is determined by ease of amplification and sequencing, confidence of assessment in phylogenetic character alignment, and by variability across broad taxon sampling. Often, a compromise must be made between using relatively highly conserved coding regions or highly variable introns and intergenic spacers. Analyses of a combination of these types of DNA regions yield phylogenetic structure at various levels of a tree (i.e., along the spine and at the tips of the branches). Here, we demonstrate the phylogenetic utility of a heretofore unused portion of a plastid protein-coding gene, hypothetical chloroplast open reading frame 1 (ycf1), in orchids. All portions of ycf1 examined are highly variable, yet alignable across Orchidaceae, and are phylogenetically informative at the level of species. In Orchidaceae, ycf1 is more variable than matK both in total number of parsimony informative characters and in percent variability. The nrITS region is more variable than ycf1, but is more difficult to align. Although we only demonstrate the phylogenetic utility of ycf1 in orchids, it is likely to be similarly useful among other plant taxa.     

Molecular phylogenetics of Caryophyllales based on nuclear 18S rDNA and plastid rbcL, atpB, and matK DNA sequences

To study the inter- and infrafamilial phylogenetic relationships in the order Caryophyllales sensu lato (s.l.), ～930 base pairs of the matK plastid gene have been sequenced and analyzed for 127 taxa. In addition, these sequences have been combined with the rbcL plastid gene for 53 taxa and with the rbcL and atpB plastid genes as well as the nuclear 18S rDNA for 26 taxa to provide increased support for deeper branches. The red pigments of Corbichonia, Lophiocarpus, and Sarcobatus have been tested and shown to belong to the betacyanin class of compounds. Most taxa of the order are clearly grouped into two main clades (i.e., "core" and "noncore" Caryophyllales) which are, in turn, divided into well-defined subunits. Phytolaccaceae and Molluginaceae are polyphyletic, and Portulacaceae are paraphyletic, whereas Agdestidaceae, Barbeuiaceae, Petiveriaceae, and Sarcobataceae should be given familial recognition. Two additional lineages are potentially appropriate to be elevated to the family level in the future: the genera Lophiocarpus and Corbichonia form a well-supported clade on the basis of molecular and chemical evidence, and Limeum appears to be separated from other Molluginaceae based on both molecular and ultrastructural data.     

Molecular systematics of Malpighiaceae: evidence from plastid rbcL and matK sequences

Phylogenetic analyses of DNA nucleotide sequences from the plastid genes rbcL and matK were employed to investigate intergeneric relationships within Malpighiaceae. Cladistic relationships generated from the independent data matrices for the family are generally in agreement with those from the combined matrix. At the base of Malpighiaceae are several clades mostly representing genera from a paraphyletic subfamily Byrsonimoideae. Intergeneric relationships among these byrsonimoid malpighs are well supported by the bootstrap, and the tribe Galphimeae is monophyletic. There is also a well-supported clade of genera corresponding to tribes Banisterieae plus Gaudichaudieae present in all trees, and many of the relationships among these banisterioid malpighs are well supported by the bootstrap. However, tribes Hiraeae and Tricomarieae (the hiraeoid malpighs) are paraphyletic and largely unresolved. Species of Mascagnia are distributed throughout these hiraeoid clades, confirming the suspected polyphyly of this large genus. Optimization of selected morphological characters on these trees demonstrates clear phylogenetic trends such as the evolution of globally symmetrical from radially symmetrical pollen, increased modification and sterilization of stamens, and switch from base chromosome number n = 6 to n = 10.     

Phylogenetic analyses of Zostera species based on rbcL and matK nucleotide sequences: implications for the origin and diversification of seagrasses in Japanese waters

Seagrasses are composed of four families belonging to angiosperms and they are thought to become adaptive to aquatic life independently. Zosteraceae is one such family and because of the relatively high species diversity around Japan and Korea coast areas, the family might have arisen therefrom. To elucidate the origin and evolution of Zosteraceae which consists of three genera, Phyllospadix, Zostera, and Heterozostera, 2.8 kb nucleotide sequences of rbcL and matK genes in the chloroplast genome were examined for various species, including cosmopolitan Z. marina and endemic Z. caulescens. The phylogenetic analysis reveals the following three features. First, based on the synonymous nucleotide substitution rate of the rice chloroplast genome, we estimated the divergence times between Zosteraceae and its closest relative, Potamogetonaceae, and between different genera, Zostera and Phyllospadix, as approximately 100 million years (myr) and 36 myr, respectively, suggesting that Zosteraceae emerged somewhere in the period from 36 myr ago to 100 myr ago. Second, two subgenera of Zostera, Zostera and Zosterella, exhibit their reciprocal monophyly and appear to have differentiated from each other approximately 33 myr ago. However, the third genus Heterozostera branched off only 5 myr ago from the stem lineage leading to Zosterella and this seems too recent in comparison with the ancient divergence of the two subgenera. Third, we estimated the most recent common ancestor of subgenus Zostera as 6 myr. In Z. marina four haplotypes were found in the sample and have diversified in the past 1.5 myr. One haplotype is shared by both sides of the Japan Archipelago and its closely related haplotypes occur also in eastern Pacific Ocean. Based on these phylogeographic analyses, we propose a provisional age related classification of Zosteraceae to argue the origin and evolution.     

Phylogenetic studies of Ophioglossaceae: evidence from rbcL and trnL-F plastid DNA sequences and morphology

Ophioglossaceae are a putatively ancient lineage of ferns in which the aerial portion of the plant is composed of a single leaf. The simplicity of foliar morphology has limited the number of characters available for constructing classifications and contributed to taxonomic difficulties at nearly every level of classification within the family. Analysis of plastid DNA rbcL sequences from 36 species representing the diversity of Ophioglossaceae supported the monophyly of the family. Intrafamilial relationships were examined using rbcL and trnL-F plastid DNA sequences and morphological data. Individual and combined analyses of the three data sets revealed two main clades within the family, here termed ophioglossoid and botrychioid. In the botrychioid clade, Helminthostachys was sister to a broadly defined Botrychium, within which Botrychium in the narrow sense of some authors and Sceptridium were sister. Botrypus was paraphyletic, with Botrypus virginianus sister to Botrychium plus Sceptridium, and with Botrypus strictus sister to all other botrychioid species except Helminthostachys. In the ophioglossoid clade, Ophioglossum in the narrow sense was sister to Cheiroglossa plus Ophioderma, but relationships within Ophioglossum were not well supported.     

Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms

Sequences of spacers and group I introns in plant chloroplast genomes have recently been shown to be very effective in phylogenetic reconstruction at higher taxonomic levels and not only for inferring relationships among species. Group II introns, being more frequent in those genomes than group I introns, may be further promising markers. Because group II introns are structurally constrained, we assumed that sequences of a group II intron should be alignable across seed plants. We designed universal amplification primers for the petD intron and sequenced this intron in a representative selection of 47 angiosperms and three gymnosperms. Our sampling of taxa is the most representative of major seed plant lineages to date for group II introns. Through differential analysis of structural partitions, we studied patterns of molecular evolution and their contribution to phylogenetic signal. Nonpairing stretches (loops, bulges, and interhelical nucleotides) were considerably more variable in both substitutions and indels than in helical elements. Differences among the domains are basically a function of their structural composition. After the exclusion of four mutational hotspots accounting for less than 18% of sequence length, which are located in loops of domains I and IV, all sequences could be aligned unambiguously across seed plants. Microstructural changes predominantly occurred in loop regions and are mostly simple sequence repeats. An indel matrix comprising 241 characters revealed microstructural changes to be of lower homoplasy than are substitutions. In showing Amborella first branching and providing support for a magnoliid clade through a synapomorphic indel, the petD data set proved effective in testing between alternative hypotheses on the basal nodes of the angiosperm tree. Within angiosperms, group II introns offer phylogenetic signal that is intermediate in information content between that of spacers and group I introns on the one hand and coding sequences on the other.     

Molecular identification and phylogenetic analysis of Aloe shadensis from Saudi Arabia based on matK,rbcL and ITS DNA barcode sequence

The Kingdom of Saudi Arabia thrives with great plant diversity, including rare plants of the family Asphodelaceae that have multiple benefits and are still being studied. Aloe shadensis is one of these plants that must be preserved and documented in its natural environment. The most appropriate molecular approach currently approved for documentation is the sequencing of some genomic markers. The current study is the first to use genomic markers to record this rare plant. In this study, the plastid genes matK (Maturase K), rbcL (Ribulose-bisphosphate carboxylase/oxygenase large subunit), and the nuclear region ITS (Internal transcribed spacer) were used to reveal their efficiency in identifying the plant under study. This study is the first to deal with this plant and document it using these genetic markers. The study showed a promising result concerning identifying the sequence of the matK gene and ITS region, while the rbcL gene did not give a good indicator through the used primers. The obtained sequences of the matK gene and the ITS region were determined through two different sets of primers in each case then deposited in GenBank. The evolutionary relatedness of Aloe shadensis was established with the different species of Aloe. The study showed that the closest species is Aloe vera with a similarity of more than 99 %. The study concludes with the possibility of using these genes to correctly identify, distinguish and document the species of Aloe shadensis.     

薯蓣属根状茎组基于3个cpDNA序列的分子系统研究

Seventeen species, one subspecies and one variety of Dioscorea sect. Stenophora Uline were investigated for their phylogenetic relationships based on a sequence analysis of chloroplast matK and rbcL genes and trnL-F intergenic spacer by maximum parsimony and maximum likelihood methods. The results showed that (a) sect. Stenophora was a strongly supported monophyletic group; (b) D. rockii, D. membranacea, D. banzhuana, and D. simulans formed a moderately supported monophyletic group, and D. prazeri was weakly supported to be sister to this group; (c) D. althaeoides and D. nipponica ssp. nipponica formed a moderately supported clade, and D. nipponica ssp. rosthornii was not a member of this clade; (d) D. zingiberensis and D. sinoparviflora showed a moderate to strong sister relationship; and (e) D. collettii var. hypoglauca and D. collettii var. collettii were sister to each other, but with only weak support.     

Specialized structures in the leaf epidermis of basal angiosperms: morphology, distribution, and homology

The morphology of specialized structures in the leaf epidermis of 32 species of basal (ANITA: Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileyaceae) angiosperms, representing all seven families and 11 of 14 genera, was investigated using light and scanning electron microscopy. Distribution, density, and size of structures were also measured, and character evolution was analyzed. Hydropotes are a synapomorphy of Nymphaeales and ethereal oil cells are a synapomorphy of Austrobaileyales, but uniseriate nonglandular trichomes appear to have arisen independently several times. Specialized structures are frequently characterized by adjacent epidermal cells that have striking similarities in their form and arrangement (i.e., architecture) to subsidiary cells of certain types of stomatal complexes. Additionally, forms intermediate to oil cells and stomata, to trichomes and stomata, and to hydropotes and oil cells are present in some taxa. Thus, all of these specialized structures and their adjacent epidermal cells form complexes that may be homologous with, and evolutionarily derived from stomatal complexes, and the specialized structure, or portion thereof, may be homologous to the stoma or guard mother cell. Improved knowledge of the morphology and evolution of these structures in the earliest branching extant angiosperm lineages has a bearing on many diverse areas of botany.     

Phylogenetic relationships of Empetraceae inferred from sequences of chloroplast gene matK and nuclear ribosomal DNA ITS region

We used sequences of nrDNA ITS and chloroplast gene matK to evaluate the monophyly of Empetrum and Corema and to examine phylogenetic relationships of the Empetraceae. Sequences of these two DNA markers were obtained for 11 plant samples, representing species of Empetrum from both the Southern and Northern Hemispheres, species and subspecies of Corema, and the monotypic Ceratiola. Sequences of four species of Rhododendron were used for rooting purposes. Our results show that species of Empetrum form a clade sister to the clade containing both Corema and Ceratiola. These two clades are strongly supported in both the matK and ITS trees, suggesting that Ceratiola is more closely related to Corema than to Empetrum, and is not of a hybrid origin between the ancestors of the latter two genera. In the matK tree, Corema conradii is more closely related to Ceratiola than to Corema album and C. album subsp. azoricum, whereas in the ITS tree, Ceratiola is allied with Corema album and C. album subsp. azoricum. This suggests that C. conradii might be a hybrid between ancestral populations of Ceratiola and C. album. The monophyly of Empetrum rejects the hypothesis of its independent origin in the two Hemispheres. Our trees also suggest the fact that the modern amphitropical distribution of Empetrum is the result of long distance dispersal, not of the vicarious events.     

Nuclear DNA amounts in angiosperms: targets, trends and tomorrow

BACKGROUND AND AIMS: The amount of DNA in an unreplicated gametic chromosome complement is known as the C-value and is a key biodiversity character of fundamental significance with many practical and predictive uses. Since 1976, Bennett and colleagues have assembled eight compilations of angiosperm C-values for reference purposes and subsequently these have been pooled into the Angiosperm DNA C-values Database (http://data.kew.org/cvalues/). Since the last compilation was published in 2005, a large amount of data on angiosperm genome size has been published. It is therefore timely to bring these data together into a ninth compilation of DNA amounts. Scope The present work lists DNA C-values for 2221 species from 151 original sources (including first values for 1860 species not listed in previous compilations). Combining these data with those published previously shows that C-values are now available for 6287 angiosperm species. KEY FINDINGS: Analysis of the dataset, which is by far the largest of the nine compilations published since 1976, shows that angiosperm C-values are now being generated at the highest rate since the first genome sizes were estimated in the 1950s. The compilation includes new record holders for the smallest (1C = 0·0648 pg in Genlisea margaretae) and largest (1C = 152·23 pg in Paris japonica) genome sizes so far reported, extending the range encountered in angiosperms to nearly 2400-fold. A review of progress in meeting targets set at the Plant Genome Size meetings shows that although representation for genera, geographical regions and some plant life forms (e.g. island floras and parasitic plants) has improved, progress to increase familial representation is still slow. In terms of technique it is now clear that flow cytometry is soon likely to become the only method available for plant genome size estimations. Fortunately, this has been accompanied by numerous careful studies to improve the quality of data generated using this technique (e.g. design of new buffers, increased awareness and understanding of problems caused by cytosolic inhibitors). It is also clear that although the speed of DNA sequencing continues to rise dramatically with the advent of next-generation and third-generation sequencing technologies, 'complete genome sequencing' projects are still unable to generate accurate plant genome size estimates.     

Nuclear DNA amounts in angiosperms: progress, problems and prospects

BACKGROUND: The nuclear DNA amount in an unreplicated haploid chromosome complement (1C-value) is a key diversity character with many uses. Angiosperm C-values have been listed for reference purposes since 1976, and pooled in an electronic database since 1997 (http://www.kew.org/cval/homepage). Such lists are cited frequently and provide data for many comparative studies. The last compilation was published in 2000, so a further supplementary list is timely to monitor progress against targets set at the first plant genome size workshop in 1997 and to facilitate new goal setting. SCOPE: The present work lists DNA C-values for 804 species including first values for 628 species from 88 original sources, not included in any previous compilation, plus additional values for 176 species included in a previous compilation. CONCLUSIONS: 1998-2002 saw striking progress in our knowledge of angiosperm C-values. At least 1700 first values for species were measured (the most in any five-year period) and familial representation rose from 30 % to 50 %. The loss of many densitometers used to measure DNA C-values proved less serious than feared, owing to the development of relatively inexpensive flow cytometers and computer-based image analysis systems. New uses of the term genome (e.g. in 'complete' genome sequencing) can cause confusion. The Arabidopsis Genome Initiative C-value for Arabidopsis thaliana (125 Mb) was a gross underestimate, and an exact C-value based on genome sequencing alone is unlikely to be obtained soon for any angiosperm. Lack of this expected benchmark poses a quandary as to what to use as the basal calibration standard for angiosperms. The next decade offers exciting prospects for angiosperm genome size research. The database (http://www.kew.org/cval/homepage) should become sufficiently representative of the global flora to answer most questions without needing new estimations. DNA amount variation will remain a key interest as an integrated strand of holistic genomics.     

Is informed consent possible in the rapidly evolving world of DNA sampling?

Ethical concerns about the human genome diversity project were discussed in Montreal last year during the 1st International Conference on DNA Sampling and Banking. This, the first article in a 2-part series on the conference, examines issues related to informed consent.     

Low levels of intraspecific genetic variation at a rapidly evolving chloroplast dna locus in North American duckweeds (Lemnaceae)

Although most previous studies on chloroplast (cp) DNA variation in plants have concentrated on systematics and evolution above the species level, intraspecific variation in cpDNA is common and has provided useful insights into population-level evolutionary processes. Polymerase chain reaction methods were used to examine restriction site and sequence variation in the chloroplast rpLI6 gene within and among populations of duckweed species (Spirodela and Lemna) from the southern and eastern United States. To our knowledge, the rpL16 region has not previously been used to investigate cpDNA variation in nature. While considerable restriction site and sequence variation were detected among species, no variation was found within populations of either of the two species (S. punctata and L. minor) selected for sequence analysis, and S. punctata showed no interpopulational variation. Two cpDNA haplotypes were identified in L. minor, with one haplotype restricted to three sites in Louisiana and the other found in all other populations sampled. This paucity of variation cannot be readily explained as the result of a low mutation rate. In general, group II introns appear to be subject to very little functional constraint, and extensive sequence differences have been found between species in the chloroplast rpL16 intron in particular. However, factors such as historical range expansions and contractions, founding effects, fluctuations in local population size, and natural selection may play a role in reducing cpDNA sequence variability in these species.