Phylogenomic analyses reveal intractable evolutionary history of a temperate bamboo genus (Poaceae: Bambusoideae)

Shibataea is a genus of temperate bamboos(Poaceae:Bambusoideae)endemic to China,but little is known about its phylogenetic position and interspecific relationships.To elucidate the phylogenetic relationship of the bamboo genus Shibataea,we performed genome-scale phylogenetic analysis of all seven species and one variety of the genus using double digest restriction-site associated DNA sequencing(dd RAD-seq)and whole plastid genomes generated using genome skimming.Our phylogenomic analyses based on dd RAD-seq and plastome data congruently recovered Shibataea as monophyletic.The nuclear data resolved S.hispida as the earliest diverged species,followed by S.chinensis,while the rest of Shibataea can be further divided into two clades.However,the plastid and nuclear topologies conflict significantly.By comparing the results of network analysis and topologies reconstructed from different datasets,we identify S.kumasasa as the most admixed species,which may be caused by incomplete lineage sorting(ILS)or interspecific gene flow with four sympatric species.This study highlights the power of dd RAD and plastome data in resolving complex relationships in the intractable bamboo genus.     

Complete DNA sequences of the plastid genomes of two parasitic flowering plant species, <Emphasis Type="Italic">Cuscuta reflexa</Emphasis> and <Emphasis Type="Italic">Cuscuta gronovii</Emphasis>

Background  

The holoparasitic plant genus Cuscuta comprises species with photosynthetic capacity and functional chloroplasts as well as achlorophyllous and intermediate forms with restricted photosynthetic activity and degenerated chloroplasts. Previous data indicated significant differences with respect to the plastid genome coding capacity in different Cuscuta species that could correlate with their photosynthetic activity. In order to shed light on the molecular changes accompanying the parasitic lifestyle, we sequenced the plastid chromosomes of the two species Cuscuta reflexa and Cuscuta gronovii. Both species are capable of performing photosynthesis, albeit with varying efficiencies. Together with the plastid genome of Epifagus virginiana, an achlorophyllous parasitic plant whose plastid genome has been sequenced, these species represent a series of progression towards total dependency on the host plant, ranging from reduced levels of photosynthesis in C. reflexa to a restricted photosynthetic activity and degenerated chloroplasts in C. gronovii to an achlorophyllous state in E. virginiana.     

Complete Plastid Genome Sequence of the Basal Asterid Ardisia polysticta Miq. and Comparative Analyses of Asterid Plastid Genomes

Ardisia is a basal asterid genus well known for its medicinal values and has the potential for development of novel phytopharmaceuticals. In this genus of nearly 500 species, many ornamental species are commonly grown worldwide and some have become invasive species that caused ecological problems. As there is no completed plastid genome (plastome) sequence in related taxa, we sequenced and characterized the plastome of Ardisia polysticta to find plastid markers of potential utility for phylogenetic analyses at low taxonomic levels. The complete A. polysticta plastome is 156,506 bp in length and has gene content and organization typical of most asterids and other angiosperms. We identified seven intergenic regions as potentially informative markers with resolution for interspecific relationships. Additionally, we characterized the diversity of asterid plastomes with respect to GC content, plastome organization, gene content, and repetitive sequences through comparative analyses. The results demonstrated that the genome organizations near the boundaries between inverted repeats (IRs) and single-copy regions (SCs) are polymorphic. The boundary organization found in Ardisia appears to be the most common type among asterids, while six other types are also found in various asterid lineages. In general, the repetitive sequences in genic regions tend to be more conserved, whereas those in noncoding regions are usually lineage-specific. Finally, we inferred the whole-plastome phylogeny with the available asterid sequences. With the improvement in taxon sampling of asterid orders and families, our result highlights the uncertainty of the position of Gentianales within euasterids I.     

The first complete plastid genome of Burmannia disticha L. from the mycoheterotrophic monocot family Burmanniaceae

Burmanniaceae is one major group within the monocot order Dioscoreales that has not had its plastome sequenced. Members of Burmanniaceae are mostly achlorophyllous, although the genus Burmannia also includes autotrophs. Here, we report sequencing and analysis of the first Burmanniaceae plastid genome from Burmannia disticha L.. This plastome is 157,480 bp and was assembled as a circular sequence with the typical quadripartite structure of plant plastid genomes. This plastome has a regular number of potentially functional genes with a total of 111, including 78 protein coding genes, 4 ribosomal RNA (rRNA) genes, and 29 tRNA genes. The ratio of the total length of genic:intergenic DNA is 1.58:1, and the mean length of intergenic regions is 398 bp, the longest being 1918 bp. The overall GC content of the B. disticha plastome is 34.90%, and the IR regions in B. disticha are more GC rich (39.50%) than the LSC (32.30%) and SSC (28.80%) regions. Phylogenetic analysis of protein-coding sequences from plastomes of related species in the order Dioscoreales support a clade comprising Burmanniaceae and Dioscoreaceae. This phylogenetic placement is congruent with previous findings based on nuclear and mitochondrial evidence.     

Complete Plastid Genome Sequencing of Four Tilia Species (Malvaceae): A Comparative Analysis and Phylogenetic Implications

Tilia is an ecologically and economically important genus in the family Malvaceae. However, there is no complete plastid genome of Tilia sequenced to date, and the taxonomy of Tilia is difficult owing to frequent hybridization and polyploidization. A well-supported interspecific relationships of this genus is not available due to limited informative sites from the commonly used molecular markers. We report here the complete plastid genome sequences of four Tilia species determined by the Illumina technology. The Tilia plastid genome is 162,653 bp to 162,796 bp in length, encoding 113 unique genes and a total number of 130 genes. The gene order and organization of the Tilia plastid genome exhibits the general structure of angiosperms and is very similar to other published plastid genomes of Malvaceae. As other long-lived tree genera, the sequence divergence among the four Tilia plastid genomes is very low. And we analyzed the nucleotide substitution patterns and the evolution of insertions and deletions in the Tilia plastid genomes. Finally, we build a phylogeny of the four sampled Tilia species with high supports using plastid phylogenomics, suggesting that it is an efficient way to resolve the phylogenetic relationships of this genus.     

The complete nucleotide sequences of the five genetically distinct plastid genomes of Oenothera, subsection Oenothera: I. sequence evaluation and plastome evolution

The flowering plant genus Oenothera is uniquely suited for studying molecular mechanisms of speciation. It assembles an intriguing combination of genetic features, including permanent translocation heterozygosity, biparental transmission of plastids, and a general interfertility of well-defined species. This allows an exchange of plastids and nuclei between species often resulting in plastome–genome incompatibility. For evaluation of its molecular determinants we present the complete nucleotide sequences of the five basic, genetically distinguishable plastid chromosomes of subsection Oenothera (=Euoenothera) of the genus, which are associated in distinct combinations with six basic genomes. Sizes of the chromosomes range from 163 365 bp (plastome IV) to 165 728 bp (plastome I), display between 96.3% and 98.6% sequence similarity and encode a total of 113 unique genes. Plastome diversification is caused by an abundance of nucleotide substitutions, small insertions, deletions and repetitions. The five plastomes deviate from the general ancestral design of plastid chromosomes of vascular plants by a subsection-specific 56 kb inversion within the large single-copy segment. This inversion disrupted operon structures and predates the divergence of the subsection presumably 1 My ago. Phylogenetic relationships suggest plastomes I–III in one clade, while plastome IV appears to be closest to the common ancestor.     

Chloroplast nucleoids are highly dynamic in ploidy,number, and structure during angiosperm leaf development

Chloroplast nucleoids are large, compact nucleoprotein structures containing multiple copies of the plastid genome. Studies on structural and quantitative changes of plastid DNA (ptDNA) during leaf development are scarce and have produced controversial data. We have systematically investigated nucleoid dynamics and ptDNA quantities in the mesophyll of Arabidopsis, tobacco, sugar beet, and maize from the early post‐meristematic stage until necrosis. DNA of individual nucleoids was quantified by DAPI‐based supersensitive epifluorescence microscopy. Nucleoids occurred in scattered, stacked, or ring‐shaped arrangements and in recurring patterns during leaf development that was remarkably similar between the species studied. Nucleoids per organelle varied from a few in meristematic plastids to >30 in mature chloroplasts (corresponding to about 20–750 nucleoids per cell). Nucleoid ploidies ranged from haploid to >20‐fold even within individual organelles, with average values between 2.6‐fold and 6.7‐fold and little changes during leaf development. DNA quantities per organelle increased gradually from about a dozen plastome copies in tiny plastids of apex cells to 70–130 copies in chloroplasts of about 7 μm diameter in mature mesophyll tissue, and from about 80 plastome copies in meristematic cells to 2600–3300 copies in mature diploid mesophyll cells without conspicuous decline during leaf development. Pulsed‐field electrophoresis, restriction of high‐molecular‐weight DNA from chloroplasts and gerontoplasts, and CsCl equilibrium centrifugation of single‐stranded and double‐stranded ptDNA revealed no noticeable fragmentation of the organelle DNA during leaf development, implying that plastid genomes in mesophyll tissues are remarkably stable until senescence.     

Exploring the plastid genome disparity of liverworts

Sequencing the plastid genomes of land plants provides crucial improvements to our understanding of the plastome evolution of land plants. Although the number of available complete plastid genome sequences has rapidly increased in the recent years, only a few sequences have been yet released for the three bryophyte lineages, namely hornworts, liverworts, and mosses. Here, we explore the disparity of the plastome structure of liverworts by increasing the number of sequenced liverwort plastomes from five to 18. The expanded sampling included representatives of all major lineages of liverworts including the genus Haplomitrium. The disparity of the liverwort genomes was compared with other 2386 land plant plastomes with emphasis on genome size and GC‐content. We found evidence for structural conservatism of the plastid genomes in liverworts and a trend towards reduced plastome sequence length in liverworts and derived mosses compared to other land plants, including hornworts and basal lineages of mosses. Furthermore, Aneura and Haplomitrium were distinct from other liverworts by an increased GC content, with the one found in Haplomitrium only second to the lycophyte Selaginella. The results suggest the hypothesis that liverworts and other land plants inherited and conserved the plastome structure of their most recent algal ancestors.     

Passiflora plastome sequencing reveals widespread genomic rearrangements

Although past studies have included Passiflora among angiosperm lineages with highly rearranged plastid genomes (plastomes), knowledge about plastome organization in the genus is limited. So far only one draft and one complete plastome have been published. Expanded sampling of Passiflora plastomes is needed to understand the extent of the genomic rearrangement in the genus, which is also unusual in having biparental plastid inheritance and plastome‐genome incompatibility. We sequenced 15 Passiflora plastomes using either Illumina paired‐end or shotgun cloning and Sanger sequencing approaches. Assembled plastomes were annotated using Dual Organellar GenoMe Annotator (DOGMA) and tRNAscan‐SE. The Populus trichocarpa plastome was used as a reference to estimate genomic rearrangements in Passiflora by performing whole genome alignment in progressiveMauve. The phylogenetic distribution of rearrangements was plotted on the maximum likelihood tree generated from 64 plastid encoded protein genes. Inverted repeat (IR) expansion/contraction and loss of the two largest hypothetical open reading frames, ycf1 and ycf2, account for most plastome size variation, which ranges from 139 262 base pairs (bp) in P. biflora to 161 494 bp in P. pittieri. Passiflora plastomes have experienced numerous inversions, gene and intron losses along with multiple independent IR expansions and contractions resulting in a distinct organization in each of the three subgenera examined. Each Passiflora subgenus has a unique plastome structure in terms of gene content, order and size. The phylogenetic distribution of rearrangements shows that Passiflora has experienced widespread genomic changes, suggesting that such events may not be reliable phylogenetic markers.     

The evolution of chloroplast genes and genomes in ferns

Most of the publicly available data on chloroplast (plastid) genes and genomes come from seed plants, with relatively little information from their sister group, the ferns. Here we describe several broad evolutionary patterns and processes in fern plastid genomes (plastomes), and we include some new plastome sequence data. We review what we know about the evolutionary history of plastome structure across the fern phylogeny and we compare plastome organization and patterns of evolution in ferns to those in seed plants. A large clade of ferns is characterized by a plastome that has been reorganized with respect to the ancestral gene order (a similar order that is ancestral in seed plants). We review the sequence of inversions that gave rise to this organization. We also explore global nucleotide substitution patterns in ferns versus those found in seed plants across plastid genes, and we review the high levels of RNA editing observed in fern plastomes.     

Nucleomorph genome diversity and its phylogenetic implications in cryptomonad algae

The relationship between phylogeny and nucleomorph genome size was examined in 16 strains of cryptomonad algae using pulsed‐field gel electrophoresis, Southern hybridization and phylogenetic analyses. Our results suggest that all cryptomonads examined in this study contain three nucleomorph chromosomes and their total genome size ranges from 495 to 750 kb. In addition, we estimated the plastid genome size of the respective organisms. The plastid genomes of photosynthetic strains were approximately 120–160 kb in size, whereas the non‐photosynthetic Cryptomonas paramecium NIES715 possesses a genome of approximately 70 kb. Phylogenetic analysis of the nuclear small subunit ribosomal DNA (SSU rDNA) gene showed that nucleomorph genome size varies considerably within closely related strains. This result indicates that the reduction of nucleomorph genomes is a rapid phenomenon that occurred multiple times independently during cryptomonad evolution. The nucleomorph genome sizes of Cryptomonas rostratiformis NIES277 appeared to be approximately 495 kb. This is smaller than that of Guillardia theta CCMP327, which until now was thought to have the smallest known nucleomorph genome size among photosynthetic cryptomonads.     

Testing the phylogenetic position of a parasitic plant (Cuscuta, Convolvulaceae, asteridae): Bayesian inference and the parametric bootstrap on data drawn from three genomes

Previous findings on structural rearrangements in the chloroplast genome of Cuscuta (dodder), the only parasitic genus in the morning-glory family, Convolvulaceae, were attributed to its parasitic life style, but without proper comparison to related nonparasitic members of the family. Before molecular evolutionary questions regarding genome evolution can be answered, the phylogenetic problems within the family need to be resolved. However, the phylogenetic position of parasitic angiosperms and their precise relationship to nonparasitic relatives are difficult to infer. Problems are encountered with both morphological and molecular evidence. Molecular data have been used in numerous studies to elucidate relationships of parasitic taxa, despite accelerated rates of sequence evolution. To address the question of the position of the genus Cuscuta within Convolvulaceae, we generated a new molecular data set consisting of mitochondrial (atpA) and nuclear (RPB2) genes, and analyzed these data together with an existing chloroplast data matrix (rbcL, atpB, trnL-F, and psbE-J), to which an additional chloroplast gene (rpl2) was added. This data set was analyzed with an array of phylogenetic methods, including Bayesian analysis, maximum likelihood, and maximum parsimony. Further exploration of data was done by using methods of phylogeny hypothesis testing. At least two nonparasitic lineages are shown to diverge within the Convolvulaceae before Cuscuta. However, the exact sister group of Cuscuta could not be ascertained, even though many alternatives were rejected with confidence. Caution is therefore warranted when interpreting the causes of molecular evolution in Cuscuta. Detailed comparisons with nonparasitic Convolvulaceae are necessary before firm conclusions can be reached regarding the effects of the parasitic mode of life on patterns of molecular evolution in Cuscuta.     

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.     

Evolutionary processes in the genus Coreocarpus: insights from molecular phylogenetics

A molecular phylogenetic study of the plant genus Coreocarpus was conducted using nuclear (ITS) and plastid (rpl16 intron) DNA sequences, with phylogenies of the nuclear and plastid sequences highly congruent in defining a monophyletic group of six species (core Coreocarpus), although three other species often placed within the genus were excluded. Relationships within the genus are largely but not totally concordant with prior biosystematic studies. Despite strong molecular support, no morphological characters uniting the six species of core Coreocarpus have been identified; retention of plesiomorphic characters and the genetic lability of characters are two probable factors contributing to lack of consistent defining characters. The age of the core Coreocarpus is estimated at 1 million years because the basal species is endemic to a volcanic island that emerged in the past million years. Mapping the results of earlier breeding studies on the molecular phylogeny showed that use of cross-compatibility as a criterion for species delimitation would result in the recognition of paraphyletic species. Prior field, morphological, and biosystematic studies provided no indication of past hybridization in the evolution of Coreocarpus, and species in the genus appeared to be well defined morphologically. However, three instances of incongruence were observed. Two of these were between the nuclear and plastid partitions, and the third was between the morphological species assignment of one accession and the molecular data. If hybridization accounts for incongruence between the nuclear and plastid data, it occurred between species that now appear to be cross-incompatible and allopatric. The incongruence between morphological species assignment and the molecular data could be the result of parallel fixation of characters that have a simple genetic basis. This study suggests that the evolutionary history of Coreocarpus is much more complex than indicated from prior biosystematic investigations and that biosystematic and molecular phylogenetic studies may complement each other for elucidating the evolution and phylogeny of a group.     

Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae)

Complete plastid genome (plastome) sequences and nuclear ribosomal DNA (nrDNA) regions have been proposed as candidates for the next generation of DNA barcodes for plant species discrimination. However, the efficacy of this approach still lacks comprehensive evaluation. We carried out a case study in the economically important but phylogenetically and taxonomically difficult genus Panax (Araliaceae). We generated a large data set of plastomes and nrDNA sequences from multiple accessions per species. Our data improved the phylogenetic resolution and levels of species discrimination in Panax, compared to any previous studies using standard DNA barcodes. This provides new insights into the speciation, lineage diversification and biogeography of the genus. However, both plastome and nrDNA failed to completely resolve the phylogenetic relationships in the Panax bipinnatifidus species complex, and only half of the species within it were recovered as monophyletic units. The results suggest that complete plastome and ribosomal DNA sequences can substantially increase species discriminatory power in plants, but they are not powerful enough to fully resolve phylogenetic relationships and discriminate all species, particularly in evolutionarily young and complex plant groups. To gain further resolving power for closely related species, the addition of substantial numbers of nuclear markers is likely to be required.     

A study of 28 Elymus species using repetitive DNA sequences.

Four repetitive DNA sequences cloned from the barley (Hordeum vulgare) genome and common for different Triticeae species were used for a molecular study of phylogenetic relationships among 28 Elymus species. Two wild Hordeum species (H genome), two Pseudoroegneria species (S genome), Agropyron cristatum (P genome), and Australopyrum velutinum (W genome) were included as genomic representatives for the genomes that supposedly were involved in the evolution of the genus Elymus. Our results are essentially congruent with the genomic classification system. This study demonstrates that Elymus is not a monophyletic genus. Based on an analysis of Southern blot hybridization we could discriminate between SY and SH species owing to the strong specific hybridization pattern of the H genome. Hexaploid SYH species gave a hybridization pattern similar to SH species for the same reason. The results support the genomic composition of Elymus batalinii as SYP and also indicated the presence of at least one H genome in Elymus enysii with a hitherto unknown genomic constitution. Elymus erianthus had a hybridization pattern distinctly different from all other species in the investigation. Key words : Elymus, RFLP, phylogeny, repetitive DNA.     

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine proteincoding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.