Patterns of nucleotide substitution in angiosperm cpDNA trnL (UAA)-trnF (GAA) regions

Comparative Genomics and Phylogenetic Analysis of the Chloroplast Genomes in Three Medicinal Salvia Species for Bioexploration

To systematically determine their phylogenetic relationships and develop molecular markers for species discrimination of Salvia bowleyana, S. splendens, and S. officinalis, we sequenced their chloroplast genomes using the Illumina Hiseq 2500 platform. The chloroplast genomes length of S. bowleyana, S. splendens, and S. officinalis were 151,387 bp, 150,604 bp, and 151,163 bp, respectively. The six genes ndhB, rpl2, rpl23, rps7, rps12, and ycf2 were present in the IR regions. The chloroplast genomes of S. bowleyana, S. splendens, and S. officinalis contain 29 tandem repeats; 35, 29, 24 simple-sequence repeats, and 47, 49, 40 interspersed repeats, respectively. The three specific intergenic sequences (IGS) of rps16-trnQ-UUG, trnL-UAA-trnF-GAA, and trnM-CAU-atpE were found to discriminate the 23 Salvia species. A total of 91 intergenic spacer sequences were identified through genetic distance analysis. The two specific IGS regions (trnG-GCC-trnM-CAU and ycf3-trnS-GGA) have the highest K2p value identified in the three studied Salvia species. Furthermore, the phylogenetic tree showed that the 23 Salvia species formed a monophyletic group. Two pairs of genus-specific DNA barcode primers were found. The results will provide a solid foundation to understand the phylogenetic classification of the three Salvia species. Moreover, the specific intergenic regions can provide the probability to discriminate the Salvia species between the phenotype and the distinction of gene fragments.

Evolutionary directions of single nucleotide substitutions and structural mutations in the chloroplast genomes of the family Calycanthaceae

BACKGROUND

Chloroplast genome sequence data is very useful in studying/addressing the phylogeny of plants at various taxonomic ranks. However, there are no empirical observations on the patterns, directions, and mutation rates, which are the key topics in chloroplast genome evolution. In this study, we used Calycanthaceae as a model to investigate the evolutionary patterns, directions and rates of both nucleotide substitutions and structural mutations at different taxonomic ranks.

RESULTS

There were 2861 polymorphic nucleotide sites on the five chloroplast genomes, and 98% of polymorphic sites were biallelic. There was a single-nucleotide substitution bias in chloroplast genomes. A → T or T → A (2.84%) and G → C or C → G (3.65%) were found to occur significantly less frequently than the other four transversion mutation types. Synonymous mutations kept balanced pace with nonsynonymous mutations, whereas biased directions appeared between transition and transversion mutations and among transversion mutations. Of the structural mutations, indels and repeats had obvious directions, but microsatellites and inversions were non-directional. Structural mutations increased the single nucleotide mutations rates. The mutation rates per site per year were estimated to be 0.14-0.34 × 10- 9 for nucleotide substitution at different taxonomic ranks, 0.64 × 10- 11 for indels and 1.0 × 10- 11 for repeats.

CONCLUSIONS

Our direct counts of chloroplast genome evolution events provide raw data for correctly modeling the evolution of sequence data for phylogenetic inferences.

Development of Chloroplast Genomic Resources in Chinese Yam (Dioscorea polystachya)

Chinese yam has been used both as a food and in traditional herbal medicine. Developing more effective genetic markers in this species is necessary to assess its genetic diversity and perform cultivar identification. In this study, new chloroplast genomic resources were developed using whole chloroplast genomes from six genotypes originating from different geographical locations. The Dioscorea polystachya chloroplast genome is a circular molecule consisting of two single-copy regions separated by a pair of inverted repeats. Comparative analyses of six D. polystachya chloroplast genomes revealed 141 single nucleotide polymorphisms (SNPs). Seventy simple sequence repeats (SSRs) were found in the six genotypes, including 24 polymorphic SSRs. Forty-three common indels and five small inversions were detected. Phylogenetic analysis based on the complete chloroplast genome provided the best resolution among the genotypes. Our evaluation of chloroplast genome resources among these genotypes led us to consider the complete chloroplast genome sequence of D. polystachya as a source of reliable and valuable molecular markers for revealing biogeographical structure and the extent of genetic variation in wild populations and for identifying different cultivars.

Population and phylogenomic decomposition via genotyping-by-sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic-scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.

Coexistence of minicircular and a highly rearranged mtDNA molecule suggests that recombination shapes mitochondrial genome organization

Recombination has been proposed as a possible mechanism to explain mitochondrial (mt) gene rearrangements, although the issue of whether mtDNA recombination occurs in animals has been controversial. In this study, we sequenced the entire mt genome of the megaspilid wasp Conostigmus sp., which possessed a highly rearranged mt genome. The sequence of the A+T-rich region contained a number of different types of repeats, similar to those reported previously in the nematode Meloidogyne javanica, in which recombination was discovered. In Conostigmus, we detected the end products of recombination: a range of minicircles. However, using isolated (cloned) fragments of the A+T-rich region, we established that some of these minicircles were found to be polymerase chain reaction (PCR) artifacts. It appears that regions with repeats are prone to PCR template switching or PCR jumping. Nevertheless, there is strong evidence that one minicircle is real, as amplification primers that straddle the putative breakpoint junction produce a single strong amplicon from genomic DNA but not from the cloned A+T-rich region. The results provide support for the direct link between recombination and mt gene rearrangement. Furthermore, we developed a model of recombination which is important for our understanding of mtDNA evolution.

Evolution of symbiosis in the legume genus Aeschynomene

Legumes in the genus Aeschynomene form nitrogen-fixing root nodules in association with Bradyrhizobium strains. Several aquatic and subaquatic species have the additional capacity to form stem nodules, and some of them can symbiotically interact with specific strains that do not produce the common Nod factors synthesized by all other rhizobia. The question of the emergence and evolution of these nodulation characters has been the subject of recent debate. We conducted a molecular phylogenetic analysis of 38 different Aeschynomene species. The phylogeny was reconstructed with both the chloroplast DNA trnL intron and the nuclear ribosomal DNA ITS/5.8S region. We also tested 28 Aeschynomene species for their capacity to form root and stem nodules by inoculating different rhizobial strains, including nodABC-containing strains (ORS285, USDA110) and a nodABC-lacking strain (ORS278). Maximum likelihood analyses resolved four distinct phylogenetic groups of Aeschynomene. We found that stem nodulation may have evolved several times in the genus, and that all Aeschynomene species using a Nod-independent symbiotic process clustered in the same clade. The phylogenetic approach suggested that Nod-independent nodulation has evolved once in this genus, and should be considered as a derived character, and this result is discussed with regard to previous experimental studies.

Barcoding the major Mediterranean leguminous crops by combining universal chloroplast and nuclear DNA sequence targets

The ability to discriminate all species is the ultimate target in barcoding. The Mediterranean basin is a center of origin for legumes and thus they have played a key role in feeding the Mediterranean population. It is also a region with important protected designation of origin and protected geographical indication legumes that provide income in rural areas. We evaluated the use of two chloroplast regions, trnL and rpoC1, and a nuclear internal transcriber region, ITS2, for their efficiency to barcode the main Mediterranean leguminous crops. Twenty-five legume species were studied. Plant material of pasture and legumes was obtained from the Greek GenBank and the Fodder Crops and Pastures Institute (National Agricultural Research Foundation). DNA was extracted with the Qiagen DNeasy plant mini-kit and PCR amplification was performed using the Kapa Taq DNA polymerase using primers amplifying the chloroplast trnL and rpoC1 regions or the nuclear region ITS2. PCR products were sequenced and the sequences were aligned using CLUSTAL W. Species identification based on the sequence similarity approach was performed using the GenBank database. In order to evaluate intraspecific and interspecific divergence in legumes we used Molecular Evolutionary Genetics Analysis 5 and for pairwise Kimura 2-parameter distance calculations for all 3 DNA regions (2 chloroplast regions, trnL and rpoC1, and the nuclear region ITS2). Four tree-based methods (neighbor joining and maximum parsimony, maximum likelihood, and Bayesian inference analyses) were used to exhibit the molecular identification results to represent differences as an uprooted dendrogram. Additionally, the sequence character-based method was used with DnaSP and the information from each site was treated as a character to distinguish the species from one another. The DNA regions trnL and ITS2 successfully (100%) discriminated the Mediterranean crop legume species used, while rpoC1 identified only 72% of them. Furthermore, the use of the trnL region enabled the discrimination of even very closely related species, like Phaseolus lunatus and P. coccineus or Vicia faba subsp major with V. faba subsp minor, which are so closely related that even in NCBI they were both referred as Phaseolus vulgaris and V. faba, respectively. We conclude that trnL and ITS2 are efficient DNA barcoding target regions in order to discriminate Mediterranean leguminous crops and provide a reliable and efficient tool for the scientific, agricultural and industrial community.

Rapid plant identification using species- and group-specific primers targeting chloroplast DNA

Plant identification is challenging when no morphologically assignable parts are available. There is a lack of broadly applicable methods for identifying plants in this situation, for example when roots grow in mixture and for decayed or semi-digested plant material. These difficulties have also impeded the progress made in ecological disciplines such as soil- and trophic ecology. Here, a PCR-based approach is presented which allows identifying a variety of plant taxa commonly occurring in Central European agricultural land. Based on the trnT-F cpDNA region, PCR assays were developed to identify two plant families (Poaceae and Apiaceae), the genera Trifolium and Plantago, and nine plant species: Achillea millefolium, Fagopyrum esculentum, Lolium perenne, Lupinus angustifolius, Phaseolus coccineus, Sinapis alba, Taraxacum officinale, Triticum aestivum, and Zea mays. These assays allowed identification of plants based on size-specific amplicons ranging from 116 bp to 381 bp. Their specificity and sensitivity was consistently high, enabling the detection of small amounts of plant DNA, for example, in decaying plant material and in the intestine or faeces of herbivores. To increase the efficacy of identifying plant species from large number of samples, specific primers were combined in multiplex PCRs, allowing screening for multiple species within a single reaction. The molecular assays outlined here will be applicable manifold, such as for root- and leaf litter identification, botanical trace evidence, and the analysis of herbivory.

Reconstruction of structural evolution in the trnL intron P6b loop of symbiotic Nostoc (Cyanobacteria)

In this study we reconstruct the structural evolution of the hyper-variable P6b region of the group I trnLeu intron in a monophyletic group of lichen-symbiotic Nostoc strains and establish it as a useful marker in the phylogenetic analysis of these organisms. The studied cyanobacteria occur as photosynthetic and/or nitrogen-fixing symbionts in lichen species of the diverse Nephroma guild. Phylogenetic analyses and secondary structure reconstructions are used to improve the understanding of the replication mechanisms in the P6b stem-loop and to explain the observed distribution patterns of indels. The variants of the P6b region in the Nostoc clade studied consist of different combinations of five sequence modules. The distribution of indels together with the ancestral character reconstruction performed enables the interpretation of the evolution of each sequence module. Our results indicate that the indel events are usually associated with single nucleotide changes in the P6b region and have occurred several times independently. In spite of their homoplasy, they provide phylogenetic information for closely related taxa. Thus we recognize that features of the P6b region can be used as molecular markers for species identification and phylogenetic studies involving symbiotic Nostoc cyanobacteria.

Molecular phylogeny of tribe Forsythieae (Oleaceae) based on nuclear ribosomal DNA internal transcribed spacers and plastid DNA trnL-F and matK gene sequences

The tribe Forsythieae comprises 2 genera (Forsythia and Abeliophyllum) and 14 species distributed mostly in the Far East. Although Forsythieae is considered monophyletic, with many symplesiomorphic characters, the phylogenetic status of Abeliophyllum remains controversial. We assessed the phylogenetic relationships of Forsythieae, based on a 3.3-kb plastid fragment (trnL-F region and matK gene) and nuclear internal transcribed spacer (ITS) region DNA sequences. We obtained a highly resolved and strongly supported topology with possible outgroups. The topology of the combined tree was congruent with those of the ITS region and matK gene. Maximum parsimony, maximum likelihood, and Bayesian inference tree analyses for the combined data also yielded identical relationships. Combined sequence data strongly supported the monophyly of Forsythieae and the close relationship between Fontanesia and Jasminum. Oleaceae, not Fontanesia, was found to be a sister group to Forsythieae. Moreover, the genus Abeliophyllum was distinctly independent of Forsythia. Three Forsythia lineages were suggested: (a) ONJ (ovata-nakaii-japonica clade), (b) VGE (viridissima-giraldiana-europaea), and (c) KISS (koreana-intermedia-saxatilis-suspensa). Our results indicated that F. × intermedia is not a hybrid between F. suspensa and F. viridissima, but further studies are needed to determine its taxonomic identity. Furthermore, the diverse fruit shapes in Oleaceae are assumed to be the result of parallelism or convergence.

Molecular identification of Mango, Mangifera indica L.var. totupura

Mango (>Mangifera indica) belonging to Anacardiaceae family is a fruit that grows in tropical regions. It is considered as the King of fruits. The present work was taken up to identify a tool in identifying the mango species at the molecular level. The chloroplast trnL-F region was amplified from extracted total genomic DNA using the polymerase chain reaction (PCR) and sequenced. Sequence of the dominant DGGE band revealed that Mangifera indica in tested leaves was Mangifera indica (100% similarity to the ITS sequences of Mangifera indica). This sequence was deposited in NCBI with the accession no. GQ927757.

Molecular phylogeny of the genus Luzula DC. (Juncaceae, Monocotyledones) based on plastome and nuclear ribosomal regions: a case of incongruence, incomplete lineage sorting and hybridisation

The genus Luzula consists of 115 species distributed throughout the world. Luzula is monophyletic, but species relationships within the genus are difficult to determine primarily due to the similar morphology even within geographically remote taxa (especially within the section Luzula). The plastome trnL intron, trnL-F intergenic spacer and the nuclear ribosomal ITS1-5.8S-ITS2 regions were analysed using maximum parsimony and maximum likelihood reconstruction in 93 species of Luzula. The incongruent phylogenetic signals obtained from the chloroplast and the nuclear genomes point to incomplete lineage sorting as well as recent hybridisation in this group. Although tree-building analyses revealed several well-supported lineages, the outcomes for many groups were ambiguous. In the total evidence tree, Luzula species were grouped within six main clades (1. subgenus Marlenia, 2. subgenus Pterodes except for L. pilosa, 3. sections Anthelaea and Nodulosae, 4. sections Diprophyllatae and Thyrsanochlamydeae, 5. section Alpinae except for a few species and 6. section Luzula). The subgenus Marlenia occupies the early derived lineage within the genus Luzula. The traditionally accepted subgenera Pterodes and Luzula (and its sections) appear to be non-monophyletic. A statistical parsimony network approach showed that ancient haplotypes and ribotypes co-occur with their descendants in Luzula. Furthermore, many haplotypes are shared among different species. Within the Luzula section Luzula, both recent hybridisation and incomplete lineage sorting of ancestral polymorphisms may represent potential sources of the incongruence between chloroplast and nuclear data.

Molecular analysis of Stemona plants in China based on sequences of four chloroplast DNA regions

Stemona sessilifolia, S. japonica and S. tuberosa are the three original sources of Stemonae Radix specified in the Chinese Pharmacopoeia (CP), and have been traditionally used for antitussive and insecticidal remedy. Significant variations in alkaloids composition and content, as well as different degrees of antitussive activities were found among them. In order to identify the genuine sources of Stemonae Radix accurately in genetic level, we determined the nucleotide sequences of chloroplast DNA trnL-trnF, trnH-psbA, petB-petD and trnK-rps16 regions of the species recorded in CP and S. parviflora, as well as the common counterfeits of Stemonae Radix, Asparagus species. The results revealed that the sequences of petB-petD and trnK-rps16 regions, showing relatively high substitution rate, were more informative than those of trnL-trnF and trnH-psbA regions. The sequences from all the four regions provided useful information to discriminate the three CP species from each other and from S. parviflora and the counterfeits. A phylogenetic tree reconstructed by the trnH-psbA sequences for 9 Stemona species distributed in China and Thailand showed that the three CP species belonged to the same clade, among which S. japonica and S. sessillifolia formed a sister group, showing closer relations to each other than to S. tuberosa. By contrast, S. parviflora was genetically far from the three CP species. Intra-species variations were observed in the three CP species. Especially, in S. tuberosa two types of petB-petD sequence and four types each of trnL-trnF, trnK-rps16 and trnH-psbA sequences resulted in 6 haplotypes; whereas, these differences had no relation with the different chemical types, but seemed to be consistent with geographical distribution.

Paraphyly of organelle DNAs in Cycas Sect. Asiorientales due to ancient ancestral polymorphisms

Background

This study addresses the apportionment of genetic diversity between Cycas revoluta and C. taitungensis, species that constitute the section Asiorientales and represent a unique, basal lineage of the Laurasian genus Cycas. Fossil evidence indicates divergence of the section from the rest of Cycas at least 30 million years ago. Geographically, C. taitungensis is limited to Taiwan whereas C. revoluta is found in the Ryukyu Archipelago and on mainland China.

Results

The phylogenies of ribosomal ITS region of mtDNA and the intergenic spacer between atpB and rbcL genes of cpDNA were reconstructed. Phylogenetic analyses revealed paraphyly of both loci in the two species and also in the section Asiorientales. The lack of reciprocal monophyly between these long isolated sections is likely due to persistent shared ancestral polymorphisms. Molecular dating estimated that mt- and cp DNA lineages coalesced to the most recent common ancestors (TMRCA) about 327 (mt) and 204 MYA (cp), corresponding with the divergence of cycad sections in the Mesozoic.

Conclusion

Fates of newly derived mutations of cycads follow Klopfstein et al.'s surfing model where the majority of new mutations do not spread geographically and remain at low frequencies or are eventually lost by genetic drift. Only successful 'surfing mutations' reach very high frequencies and occupy a large portion of a species range. These mutations exist as dominant cytotypes across populations and species. Geographical subdivision is lacking in both species, even though recurrent gene flow by both pollen and seed is severely limited. In total, the contrasting levels between historical and ongoing gene flow, large population sizes, a long lifespan, and slow mutation rates in both organelle DNAs have all likely contributed to the unusually long duration of paraphyly in cycads.

Combined genotypic–phenotypic phylogeny of the tribe Spathicarpeae (Araceae) with reference to independent events of invasion to Andean regions

The Neotropical tribe Spathicarpeae is noteworthy within the family Araceae for having members occurring in very dry and/or very high (above 3000m) areas in South America. Considering the small size of previously included genera in the tribe (10 geophytic genera, none with more than 15 described species), generic circumscription has been considered troublesome. To address these issues, we sequenced two cpDNA markers (matK and trnL-F) and reconstructed the phylogeny of the tribe using two approaches: Maximum parsimony and Bayesian Inference. We also evaluated the effect of adding phenotypic characters in a "total evidence" analysis. Our results confirm the previous suggestion that the two non-geophytic genera Dieffenbachia and Bognera should be included in Spathicarpeae, but the precise placement of Bognera within the tribe is still unresolved. Our results also support the segregation of Asterostigma pavonii and Asterostigma integrifolium, recently described as new monospecific genera (Incarum and Croatiella, respectively). The reconstructed phylogenies suggest that the invasion of the Andean range by the Spathicarpeae occurred in at least two independent events. The inclusion of phenotypic characters in a combined analysis has improved both internal resolution and the number of strongly supported clades.

Molecular evidence for an Asian origin and a unique westward migration of species in the genus Castanea via Europe to North America

The genus Castanea (Fagaceae) is widely distributed in the deciduous forests of the Northern Hemisphere. The striking similarity between the floras of eastern Asia and those of eastern North America and the difference in chestnut blight resistance among species has been of interest to botanists for a century. To infer the biogeographical history of the genus, the phylogeny of Castanea was estimated using DNA sequence data from different regions of the chloroplast genome. Sequencing results support the genus Castanea as a monophyletic group with Castanea crenata as basal. The three Chinese species form a strongly supported sister clade to the North American and European clade. A unique westward expansion of extant Castanea species is hypothesized with Castanea originating in eastern Asia, an initial diversification within Asia during the Eocene followed by intercontinental dispersion and divergence between the Chinese and the European/North American species during the middle Eocene and a split between the European and the North American species in the late Eocene. The differentiation within North America and China might have occurred in early or late Miocene. The North America species are supported as a clade with C. pumila var. ozarkensis, the Ozark chinkapin, as the basal lineage, sister to the group comprising C. pumila var. pumila, the Allegheny chinkapin, and Castanea dentata, the American chestnut. Morphological evolution of one nut per bur in the genus may have occurred independently on two continents.

Species identification using sequences of the trnL intron and the trnL-trnF IGS of chloroplast genome among popular plants in Taiwan

Forensic botanical comparison can be hampered by the lack of appropriate DNA databases. While DNA sequence databases for many mitochondrial loci have been established for the identification of animal species, less is known regarding the genomes of plants. We report on the use of the trnL intron and the trnL-trnF intergenic spacer (IGS) in the chloroplast genome and establish a DNA sequence database for plant species identification. The DNA sequences at these two loci from commonly encountered plants, including monocots and dicots, were aligned to establish a DNA database of local plants. The database comprises 373 individual sequences representing 80 families, 206 genera and 269 species. These plant species can be grouped to species level using both sequence and length polymorphisms at these loci. To validate the database for future forensic purposes, we sequenced 20 blind samples and searched the local database and the databases of GenBank and EMBL. Fifteen of these 20 samples used in blind trial testing matched their respective species from our local DNA database but only 6 matched species registered in the GenBank and EMBL databases. The sequences of two species used in the blind trial did not match any sequence registered in any of these databases. Cluster analysis was performed to demonstrate the family and genus distribution of samples. Neighbor-joining trees of the two DNA regions from 70 samples of the local database and 10 of the species used in the blind trials were constructed and clustered to both family and genus. The bootstrap values of the trnL intron were higher than most of those of the trnL-trnF IGS. The sequence database described in this study can be used to identify plant species using DNA sequences of the trnL intron and trnL-trnF IGS of chloroplast genome and illustrates its value in plant species identification.

Molecular evolution of the trnL(UAA) intron in bryophytes

Structure, variability, and molecular evolution of the trnL(UAA) intron in bryophytes (mosses and liverworts) is analyzed based on more than 1000 sequences representing all classes, including comparisons of lengths and GC-contents, sequence similarities, evolutionary rates and ti/tv ratios of the major lineages and selected genera. Secondary structure analyses of the more variable stem-loop regions facilitated recognition of sequence repeats and minute inversions that often occurred independently in non-related lineages, thus supporting alignment construction and homology assessment. The most length-variable stem-loop region P8 does not share a common evolutionary history across all major bryophyte lineages. Independent nucleotide additions such as internally repeated sequence segments resulted in non-homologous P8 sequences that cannot be folded into a common P8 secondary structure, neither for all bryophytes nor for liverworts or mosses. To address evolutionary patterns, separate analyses of P6/P8 and the remaining intron (core) have to be performed, as overall values of the complete intron are misleading. It is argued that a transition bias observed above the genus level in the core structure is caused by structural constraints, not by its higher GC-content in comparison to the more AT-rich P6 and P8. Compensating base pair changes detected in highly conserved elements are often characteristic of the major bryophyte lineages (classes). Sequence divergence and evolutionary rates are generally higher in liverworts than in mosses, resulting in ambiguous alignments of P6 and P8 even within classes. In mosses, trends towards length reduction of P8 and lower evolutionary rates of the intron are observed. Average intraspecific variation is less than 1%, corresponding to 2-3 mutations in the complete intron.

The chloroplast trnT-trnF region in the seed plant lineage Gnetales

The trnT-trnF region is located in the large single-copy region of the chloroplast genome. It consists of the trnL intron, a group I intron, and the trnT-trnL and trnL-trnF intergenic spacers. We analyzed the evolution of the region in the three genera of the gymnosperm lineage Gnetales (Gnetum, Welwitschia, and Ephedra), with especially dense sampling in Gnetum for which we sequenced 41 accessions, representing most of the 25-35 species. The trnL intron has a conserved secondary structure and contains elements that are homologous across land plants, while the spacers are so variable in length and composition that homology cannot be found even among the three genera. Palindromic sequences that form hairpin structures were detected in the trnL-trnF spacer, but neither spacer contained promoter elements for the tRNA genes. The absence of promoters, presence of hairpin structures in the trnL-trnF spacer, and high sequence variation in both spacers together suggest that trnT and trnF are independently transcribed. Our model for the expression and processing of the genes tRNA(Thr)(UGU), tRNA(Leu)(UAA), and tRNA(Phe) (GAA) therefore attributes the seemingly neutral evolution of the two spacers to their escape from functional constraints.

Molecular phylogeny, recent radiation and evolution of gross morphology of the rhubarb genus Rheum (Polygonaceae) inferred from chloroplast DNA trnL-F sequences

BACKGROUND AND AIMS

Rheum, a highly diversified genus with about 60 species, is mainly confined to the mountainous and desert regions of the Qinghai-Tibetan plateau and adjacent areas. This genus represents a good example of the extensive diversification of the temperate genera in the Qinghai-Tibetan plateau, in which the forces to drive diversification remain unknown. To date, the infrageneric classification of Rheum has been mainly based on morphological characters. However, it may have been subject to convergent evolution under habitat pressure, and the systematic position of some sections are unclear, especially Sect. Globulosa, which has globular inflorescences, and Sect. Nobilia, which has semi-translucent bracts. Recent palynological research has found substantial contradictions between exine patterns and the current classification of Rheum. Two specific objectives of this research were (1) to evaluate possible relationships of some ambiguous sections with a unique morphology, and (2) to examine possible occurrence of the radiative speciation with low genetic divergence across the total genus and the correlation between the extensive diversification time of Rheum and past geographical events, especially the recent large-scale uplifts of the Qinghai-Tibetan Plateau.

METHODS

The chloroplast DNA trnL-F region of 29 individuals representing 26 species of Rheum, belonging to seven out of eight sections, was sequenced and compared. The phylogenetic relationships were further constructed based on the sequences obtained.

KEY RESULTS

Despite the highly diversified morphology, the genetic variation in this DNA fragment is relatively low. The molecular phylogeny is highly inconsistent with gross morphology, pollen exine patterns and traditional classifications, except for identifying all samples of Sect. Palmata, three species of Sect. Spiciformia and a few species of Sect. Rheum as corresponding monophyletic groups. The monotypic Sect. Globulosa showed a tentative position within the clade comprising five species of Sect. Rheum. All of the analyses revealed the paraphyly of R. nobile and R. alexandrae, the only two species of Sect. Nobilia circumscribed by the possession of large bracts. The crude calibration of lineages based on trnL-F sequence differentiation implied an extensive diversification of Rheum within approx. 7 million years.

CONCLUSIONS

Based on these results, it is suggested that the rich geological and ecological diversity caused by the recent large-scale uplifts of the Qinghai-Tibetan Plateau since the late Tertiary, coupled with the oscillating climate of the Quaternary stage, might have promoted rapid speciation in small and isolated populations, as well as allowing the fixation of unique or rare morphological characters in Rheum. Such a rapid radiation, combined with introgressive hybridization and reticulate evolution, may have caused the transfer of cpDNA haplotypes between morphologically dissimilar species, and might account for the inconsistency between morphological classification and molecular phylogeny reported here.

Phylogeny of Eleusine (Poaceae: Chloridoideae) based on nuclear ITS and plastid trnT–trnF sequences

Phylogenetic relationships in the genus Eleusine (Poaceae: Chloridoideae) were investigated using nuclear ITS and plastid trnT-trnF sequences. Separate and combined data sets were analyzed using parsimony, distance, and likelihood based methods, including Bayesian. Data congruence was examined using character and topological measures. Significant data heterogeneity was detected, but there was little conflict in the topological substructure measures for triplets and quartets, and resolution and clade support increased in the combined analysis. Data incongruence may be a result of noise and insufficient information in the slower evolving trnT-trnF. Monophyly of Eleusine is strongly supported in all analyses, but basal relationships in the genus remain uncertain. There is good support for a CAIK clade (E. coracana subsp. coracana and africana, E. indica, and E. kigeziensis), with E. tristachya as its sister group. Two putative ITS homeologues (A and B loci) were identified in the allotetraploid E. coracana; the 'B' locus sequence type was not found in the remaining species. Eleusine coracana and its putative 'A' genome donor, the diploid E. indica, are confirmed close allies, but sequence data contradicts the hypothesis that E. floccifolia is its second genome donor. The 'B' genome donor remains unidentified and may be extinct.

TrnL- trnF intergenic spacer and trnL intron define major clades within Luzula and Juncus (Juncaceae): importance of structural mutations

Seven hundred fifty-two to one thousand ninety-seven base pairs of the trnL intron and trnL- trnF intergenic spacer of the chloroplast DNA of 55 Juncaceae taxa (Juncus, Luzula, Rostkovia, and Oxychloë) was sequenced. Seventeen structural mutations (13 indels marked A to M, 3 parts of the trnF pseudogene, and insertion "o" within a pseudogene) within the chloroplast trnL- trnF region were examined as possible indicators for phylogenetic relationships in Juncaceae. Juncus trifidus (section Steirochloa) was clearly separated from the other taxa by two large (>80 bp) indels. The "Southern Hemisphere clade" was strongly supported by a unique insertion (334 bp) in the trnL intron. The monophyly of Luzula was supported by three small (<10 bp) indels in the trnL-F spacer. They were found in all 22 examined members that represent the taxonomic and geographical diversity of the genus Luzula. A tandemly duplicated tRNA pseudogene was found in the Juncus subgenus Juncus species and is supported by four small unique indels too. The acceptor stem and D-domain-encoding regions are separated by a unique 8-bp insertion. The T-domain and acceptor stem-encoding regions were not found in the pseudogene repeats. Only the Juncus sections Ozophyllum and Iridifolii contain the 5' acceptor stem, D-domain, and anticodon domain of the tRNAF encoding DNA. The structural mutations in the trnL intron and the trnL- trnF intergenic spacer are useful for phylogenetic reconstruction in the Juncaceae.

Evolution of the trnF(GAA) Gene in Arabidopsis Relatives and the Brassicaceae Family: Monophyletic Origin and Subsequent Diversification of a Plastidic Pseudogene

Recently, we used the 5'-trnL(UAA)-trnF(GAA) region of the chloroplast DNA for phylogeographic reconstructions and phylogenetic analysis among the genera Arabidopsis, Boechera, Rorippa, Nasturtium, and Cardamine. Despite the fact that extensive gene duplications are rare among the chloroplast genome of higher plants, within these taxa the anticodon domain of the trnF(GAA) gene exhibit extensive gene duplications with one to eight tandemly repeated copies in close 5' proximity of the functional gene. Interestingly, even in Arabidopsis thaliana we found six putative pseudogenic copies of the functional trnF gene within the 5'-intergenic trnL-trnF spacer. A reexamination of trnL(UAA)-trnF(GAA) regions from numerous published phylogenetic studies among halimolobine, cardaminoid, and other cruciferous taxa revealed not only extensive trnF gene duplications but also favor the hypothesis about a single origin of trnF pseudogene formation during evolution of the Brassicaceae family 16-21 MYA. Conserved sequence motifs from this tandemly repeated region are codistributed nonrandomly throughout the plastome, and we found some similarities with a DNA sequence duplication in the rps7 gene and its adjacent spacer. Our results demonstrate the potential evolutionary dynamics of a plastidic region generally regarded as highly conserved and probably cotranscribed and, as shown here for several genera among cruciferous plants, greatly characterized by parallel gains and losses of duplicated trnF copies.

Extensive intraindividual variation in plastid rDNA sequences from the holoparasite Cynomorium coccineum (Cynomoriaceae)

Ribosomal genes are considered to have a high degree of sequence conservation between species and also at higher taxonomic levels. In this paper we document a case where a single individual of Cynomorium coccineum (Cynomoriaceae), a nonphotosynthetic holoparasitic plant, contains highly divergent plastid ribosomal genes. PCR amplification a nearly complete ribosomal DNA cistron was performed using genomic DNA, the products cloned, and the 23S rDNA genes were sequenced from 19 colonies. Of these, five distinct types were identified. Fifteen of the sequences were nearly identical (11 or fewer differences) and these were designated Type I. The remaining types (II-V) were each represented by a single clone and differed from Type I by 93 to 255 changes. Compared with green vascular plants, we found that there are more substitutional differences in the 23S rDNA sequences within a single individual of Cynomorium than among all sequenced photosynthetic vascular plants. Several trends of molecular evolution observed in 16S rDNA from other holoparasitic angiosperms and heterotrophic green algae have been also observed in Cynomorium 23S rDNA. Higher-order structures were constructed for representatives of the five clone types, and in all cases these possessed complete complements of the major structural elements present in functional plastid 23S rRNAs. These data indicate that such molecules may be subject to purifying selection, thus providing indirect evidence that they have retained some degree of functionality. This intraindividual polymorphism is probably a case of plastid heteroplasmy but translocation of ribosomal cistrons to the nucleus or mitochondria has not been tested and therefore cannot be ruled out.

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 evolution of the trnTUGU-trnFGAA region in Bryophytes

Structure, variability, and molecular evolution of the trnT-F region in the Bryophyta (mosses and liverworts) is analyzed based on about 200 sequences of the trnT-L spacer and trnL 5' exon, 1000 sequences of the trnL intron, and 800 sequences of the trnL 3' exon and trnL-F spacer, including comparisons of lengths, GC contents, sequence similarities, and functional elements. Mutations occurring in the trnL 5' and 3' exons, including compensatory base pair changes, and a transition in the trnL anticodon in Takakia lepidozioides, are discussed. All three non-coding regions display a mosaic structure of highly variable elements (V1 - V3 in the trnT-L spacer, V4/V5 corresponding to stem-loop regions P6/P8 in the trnL intron, and V6/V7 in the trnL-F spacer) and more conserved elements. In the trnL intron this structure is a consequence of the defined secondary structure necessary for correct splicing, whereas in both spacers conserved regions are restricted to promoter elements. At least the highly variable regions in the trnT-L spacer and stem-loop region P8 of the trnL intron seem to evolve independently in the major bryophyte lineages and are therefore not suitable for high taxonomic level phylogenetic reconstructions. In mosses, a trend of length reduction towards the more derived lineages is observed in all three non-coding regions. GC contents are mostly linked to sequence variability, with the conserved regions being more GC rich and the more variable AT rich. The lowest GC values (< 10 %) are found in the trnT-L spacer of mosses. In addition to two putative sigma (70)-type promoters in the trnT-L spacer, a third putative promoter is present in the trnL-F spacer, although trnL and trnF are assumed to be co-transcribed. Consensus sequences are provided for the -35 and -10 sequences of the major bryophyte lineages. The third promoter is part of a hairpin secondary structure, whose loop region is highly homoplastic in mosses due to an inversion occurring independently in non-related taxa, even at the intraspecific level.

Noncoding plastid trnT-trnF sequences reveal a well resolved phylogeny of basal angiosperms

Recent contributions from DNA sequences have revolutionized our concept of systematic relationships in angiosperms. However, parts of the angiosperm tree remain unclear. Previous studies have been based on coding or rDNA regions of relatively conserved genes. A phylogeny for basal angiosperms based on noncoding, fast-evolving sequences of the chloroplast genome region trnT-trnF is presented. The recognition of simple direct repeats allowed a robust alignment. Mutational hot spots appear to be confined to certain sectors, as in two stem-loop regions of the trnL intron secondary structure. Our highly resolved and well-supported phylogeny depicts the New Caledonian Amborella as the sister to all other angiosperms, followed by Nymphaeaceae and an Austrobaileya-Illicium-Schisandra clade. Ceratophyllum is substantiated as a close relative of monocots, as is a monophyletic eumagnoliid clade consisting of Piperales plus Winterales sister to Laurales plus Magnoliales. Possible reasons for the striking congruence between the trnT-trnF based phylogeny and phylogenies generated from combined multi-gene, multi-genome data are discussed.

Chloroplast and nuclear evidence for multiple origins of polyploids and diploids of Hedera (Araliaceae) in the Mediterranean basin

Chloroplast (trnT-L) and nuclear rDNA (ITS) sequence analyses of the Araliaceae provide strong molecular evidence for the monophyly of the genus Hedera. Phylogenetic reconstructions suggest multiple origins and an active polyploidization process not only in the formation of tetraploids (2n = 96), hexaploids (2n = 144), and octoploids (2n = 192), but also of diploids (2n = 48). A high basic chromosome number of x = 24, extensive polyphyly in widespread diploids, and terminal placement of Hedera in phylogenies of the Araliaceae reveal that extant diploid taxa may be, in fact, assemblages of ancestral polyploids from plants of n = 12. Four major lineages containing four types of chloroplast (chlorotypes I, II, III, and IV), which are defined by different trnT-L nucleotide substitutions and two large insertions (50- and 30-bp), provide evidence for evolutionary processes and historical biogeography in Hedera. We propose a scenario where an initial colonization in the Mediterranean basin by Asian ancestors (carrying the ancestral Araliaceae chlorotype I) is followed by differentiation into the four chlorotypes of the Mediterranean region, and then recolonization of Asia and northern Europe only by chlorotype III. The Macaronesian taxa (Hedera azorica, Hedera maderensis ssp. maderensis, and Hedera canariensis) appear to have originated from a single-colonization event to each archipelago with no further contact either with continental or insular species.