Analysis of Acorus calamus chloroplast genome and its phylogenetic implications

Determining the phylogenetic relationships among the major lines of angiosperms is a long-standing problem, yet the uncertainty as to the phylogenetic affinity of these lines persists. While a number of studies have suggested that the ANITA (Amborella-Nymphaeales-Illiciales-Trimeniales-Aristolochiales) grade is basal within angiosperms, studies of complete chloroplast genome sequences also suggested an alternative tree, wherein the line leading to the grasses branches first among the angiosperms. To improve taxon sampling in the existing chloroplast genome data, we sequenced the chloroplast genome of the monocot Acorus calamus. We generated a concatenated alignment (89,436 positions for 15 taxa), encompassing almost all sequences usable for phylogeny reconstruction within spermatophytes. The data still contain support for both the ANITA-basal and grasses-basal hypotheses. Using simulations we can show that were the ANITA-basal hypothesis true, parsimony (and distance-based methods with many models) would be expected to fail to recover it. The self-evident explanation for this failure appears to be a long-branch attraction (LBA) between the clade of grasses and the out-group. However, this LBA cannot explain the discrepancies observed between tree topology recovered using the maximum likelihood (ML) method and the topologies recovered using the parsimony and distance-based methods when grasses are deleted. Furthermore, the fact that neither maximum parsimony nor distance methods consistently recover the ML tree, when according to the simulations they would be expected to, when the out-group (Pinus) is deleted, suggests that either the generating tree is not correct or the best symmetric model is misspecified (or both). We demonstrate that the tree recovered under ML is extremely sensitive to model specification and that the best symmetric model is misspecified. Hence, we remain agnostic regarding phylogenetic relationships among basal angiosperm lineages.     

The chloroplast genome of the “basal” angiosperm Calycanthus fertilis – structural and phylogenetic analyses

The nucleotide sequence of the complete chloroplast genome of a basal angiosperm, Calycanthus fertilis, has been determined. The circular 153337 bp long cpDNA is colinear with those of tobacco, Arabidopsis and spinach. A total of 133 predicted genes (115 individual gene species, 18 genes duplicated in the inverted repeats) including 88 potential protein-coding genes (81 gene species), 8 ribosomal RNA genes (4 gene species) and 37 tRNA genes (30 gene species) representing 20 amino acids were identified based on similarity to their homologs from other chloroplast genomes. This is the highest gene number ever registered in an angiosperm plastome. Calycanthus fertilis cpDNA also contains a homolog of the recently discovered mitochondrial ACRS gene. Since no gene transfer from mitochondria to the chloroplast has ever been documented, we investigated the evolutionary affinity of this gene in detail. Phylogenetic analysis of the protein-coding subset of the plastome suggests that the ancient line of Laurales emerged after the split of the angiosperms into monocots and dicots. Calycanthus fertilis Walter var. ferax (Michy.) Rehder is a synonym of C. floridus L. var. glaucus (Willd.) Torr. & A. Gray.Data deposition: The sequence reported in this paper has been deposited in the EMBL database (accession no. AJ428413).     

Complete nucleotide sequence of the sugarcane (Saccharum officinarum) chloroplast genome: a comparative analysis of four monocot chloroplast genomes.

The complete nucleotide sequence of the chloroplast genome of sugarcane (Saccharum officinarum) has been determined. It is a circular double-stranded DNA molecule, 141,182 bp in size, and is composed of a large single copy of 83,048 bp, a small single copy of 12,544 bp, and a pair of inverted repeat regions of 22,795 bp each. A comparative analysis among monocots showed that the sugarcane chloroplast genome was very similar to maize but not to rice or wheat. Between sugarcane and maize at the rps16-trnQ (UUG) region, however, a length polymorphism was identified. With regard to insertions/deletions equal to or longer than 5 bp, a total of 53 insertion and 31 deletion events were identified in the sugarcane chloroplast genome. Of the 84 loci identified, a pair of direct repeat sequences was located side by side in a tandem fashion in 47 loci (56.0%). A recombination event during plant evolution is discussed at two sites between the sugarcane and tobacco chloroplast genomes.     

Analysis of the Amborella trichopoda chloroplast genome sequence suggests that amborella is not a basal angiosperm

Phylogenetic analyses based on comparison of a limited number of genes recently suggested that Amborella trichopoda is the most ancient angiosperm. Here we present the complete sequence of the chloroplast genome of this plant. It does not display any of the genes characteristic of chloroplast DNA of the gymnosperm Pinus thunbergii (chlB, chlL, chlN, psaM, and ycf12). The majority of phylogenetic analyses of protein-coding genes of this chloroplast DNA suggests that Amborella is not the basal angiosperm and not even the most basal among dicots.     

The chloroplast genome of Nymphaea alba: whole-genome analyses and the problem of identifying the most basal angiosperm

Angiosperms (flowering plants) dominate contemporary terrestrial flora with roughly 250,000 species, but their origin and early evolution are still poorly understood. In recent years, molecular evidence has accumulated suggesting a dicotyledonous origin of monocots. Phylogenetic reconstructions have suggested that several dicotyledonous groups that include taxa such as Amborella, Austrobaileya, and Nymphaea branch off as the most basal among angiosperms. This has led to the concept of monocots, "eudicots," "basal dicots," and "ANITA" groupings. Here, we present the sequence and phylogenetic analyses of the chloroplast DNA of Nymphaea alba. Phylogenetic analyses of our 14-species data set, consisting of 29,991 aligned nucleotide positions per chloroplast genome, revealed consistent support for Nymphaea being a divergent member of a monophyletic dicot assemblage. Three distinct angiosperm lineages were supported in the majority of our phylogenetic analyses-eudicots, Magnoliopsida, and monocots. However, the monocot lineage leading to the grasses was the deepest branching. Although analyses of only one individual gene alignment (out of 61) is consistent with some recently proposed hypotheses for the paraphyly of dicots, we also report observations that nine genes do not support paraphyly of dicots. Instead, they support the basal monocot-dicot split. Consistent with this finding, we also report observations suggesting that the monocot lineage leading to the grasses has the strongest phylogenetic affinity to gymnosperms. Our findings have general implications for studies of substitution model specification and analyses of concatenated genome data.     

The complete nucleotide sequence of the coffee (Coffea arabica L.) chloroplast genome: organization and implications for biotechnology and phylogenetic relationships amongst angiosperms

The chloroplast genome sequence of Coffea arabica L., the first sequenced member of the fourth largest family of angiosperms, Rubiaceae, is reported. The genome is 155 189 bp in length, including a pair of inverted repeats of 25 943 bp. Of the 130 genes present, 112 are distinct and 18 are duplicated in the inverted repeat. The coding region comprises 79 protein genes, 29 transfer RNA genes, four ribosomal RNA genes and 18 genes containing introns (three with three exons). Repeat analysis revealed five direct and three inverted repeats of 30 bp or longer with a sequence identity of 90% or more. Comparisons of the coffee chloroplast genome with sequenced genomes of the closely related family Solanaceae indicated that coffee has a portion of rps19 duplicated in the inverted repeat and an intact copy of infA . Furthermore, whole-genome comparisons identified large indels (> 500 bp) in several intergenic spacer regions and introns in the Solanaceae, including trnE (UUC)– trnT (GGU) spacer, ycf4 – cemA spacer, trnI (GAU) intron and rrn5 – trnR (ACG) spacer. Phylogenetic analyses based on the DNA sequences of 61 protein-coding genes for 35 taxa, performed using both maximum parsimony and maximum likelihood methods, strongly supported the monophyly of several major clades of angiosperms, including monocots, eudicots, rosids, asterids, eurosids II, and euasterids I and II. Coffea (Rubiaceae, Gentianales) is only the second order sampled from the euasterid I clade. The availability of the complete chloroplast genome of coffee provides regulatory and intergenic spacer sequences for utilization in chloroplast genetic engineering to improve this important crop.     

Nelumbonaceae:Systematic position and species diversification revealed by the complete chloroplast genome

<正>Nelumbonaceae is a morphologically unique family of angiosperms and was traditionally placed in Nymphaeales;more recently,it was placed in Proteales based on molecular data,or in an order of its own,Nelumbonales. To determine the systematic position of the family and to date the divergence time of the family and the divergence time of its two intercontinentally disjunct species,we sequenced the entire chloroplast genome of Nelumbo lutea and most of the chloroplast genes of,N.nucifera.We carried out phylogenetic and molecular dating analyses of the two species and representatives of 47 other plant families,representing the major lineages of angiosperms, using 83 plastid genes.The N.lutea genome was 163 510 bp long,with a total of 130 coding genes and an overall GC content of 38%.No significant structural differences among the genomes of N.lutea,Nymphaea alba, and Platanus occidentalis were observed.The phylogenetic relationships based on the 83 plastid genes revealed a close relationship between Nelumbonaceae and Platanaceae.The divergence times were estimated to be 109 Ma between the two families and 1.5 Ma between the two Nelumbo species.The estimated time was only slightly longer than the age of known Nelumbo fossils,suggesting morphological stasis within Nelumbonaceae.We conclude that Nelumbonaceae holds a position in or close to Proteales.We further conclude that the two species of Nelumbo diverged recently from a common ancestor and do not represent ancient relicts on different continents.     

The complete chloroplast genome sequence of the Chinese endemic species Sorbus setschwanensis (Rosaceae) and its phylogenetic analysis

Sorbus setschwanensis Koehne is a pinnate-leaved Sorbus s.str. species endemic to China with narrow distribution and intriguing phylogeny that needs wider attention. In this paper, the complete chloroplast (cp) genome of S. setschwanensis is reported, and its phylogenetic position is analyzed. The complete cp genome of S. setschwanensis is 160 064 bp in size with 36.50% GC content. It has a typical quadripartite structure including a pair of inverted repeat regions (IRs) of 26 378 bp that separates a large single copy (LSC) region of 86 013 bp and a small single copy (SSC) region of 19 295 bp. The cp genome encodes 108 genes, comprising 76 protein-coding genes, 28 tRNA genes and 4 rRNA genes. Additionally, 52 simple sequence repeats (SSRs) and 43 dispersed repeats were identified. Comparison of the whole cp genome with those of other Sorbus species showed an overall high degree of sequence similarity, but there are six highly variable regions (trnR-atpA, petN-psbM, ndhC-trnV, trnE-trnT, trnT-trnL and rpl32-trnL) located in intergenic spacers that may be useful as molecular markers in future population genetic and phylogenetic studies in the genus. Phylogenetic analyses based on 108 coding genes from 25 species in Rosaceae revealed that S. setschwanensis is nested within Sorbus sect. Sorbus together with other pinnately leaved species, but does not form a sister lineage to S. rufopilosa belonging to the same series Multijugae. Thus, the systematic position of S. setschwanensis and relationships with other species in the genus needs to be further studied.     

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next‐generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch‐Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.     

Genome size and karyotype evolution in the slipper orchids (Cypripedioideae: Orchidaceae)

Nuclear DNA contents (4C) were estimated by Feulgen microdensitometry in 27 species of slipper orchids. These data and recent information concerning the molecular systematics of Cypripedioideae allow an interesting re-evaluation of karyotype and genome size variation among slipper orchids in a phylogenetic context. DNA amounts differed 5.7-fold, from 24.4 pg in Phragmipedium longifolium to 138.1 pg in Paphiopedilum wardii. The most derived clades of the conduplicate-leaved slipper orchids have undergone a radical process of genome fragmentation that is most parsimoniously explained by Robertsonian changes involving centric fission. This process seems to have occurred independently of genome size variation. However, it may reflect environmental or selective pressures favoring higher numbers of linkage groups in the karyotype.     

红边龙血树叶绿体基因组特征及其系统发育分析

红边龙血树(Dracaena marginata)是一种在全球广泛种植的龙血树属园艺植物,具有较高的观赏价值和药用价值。本研究首次利用高通量测序技术对红边龙血树叶片进行全基因组测序,组装得到完整的叶绿体基因组序列,并进行注释、序列特征比较和系统发育分析。结果表明,红边龙血树叶绿体基因组包含一个典型的四分体结构,长度为154926 bp,是目前已报道的龙血树属中叶绿体基因组最小的物种;共拥有132个基因,包含86个编码蛋白基因、38个转运RNA基因和8个核糖体RNA基因;密码子偏好性分析发现存在偏好使用A/U碱基结尾的现象,整体上密码子偏好性较低;共鉴定出46个简单重复序列位点和54个长重复序列,分别在大单拷贝区与反向重复区有最大检出率;种间边界分析发现边界区域基因存在相对位置差异,扩张收缩情况总体较为相似;与近缘种进行系统发育分析,红边龙血树与细枝龙血树聚为一类,关系最近,符合形态学分类特征。对红边龙血树叶绿体基因组的解析为龙血树属植物的物种鉴定、遗传多样性和叶绿体转基因工程等提供了重要数据基础。     

Tissue distribution and biosynthesis of 1,2-saturated pyrrolizidine alkaloids in Phalaenopsis hybrids (Orchidaceae)

Phalaenopsis hybrids contain two 1,2-saturated pyrrolizidine monoesters, T-phalaenopsine (necine base trachelanthamidine) and its stereoisomer Is-phalaenopsine (necine base isoretronecanol). T-Phalaenopsine is the major alkaloid accounting for more than 90% of total alkaloid. About equal amounts of alkaloid were genuinely present as free base and its N-oxide. The structures were confirmed by GC-MS. The quantitative distribution of phalaenopsine in various organs and tissues of vegetative rosette plants and flowering plants revealed alkaloid in all tissues. The highest concentrations were found in young and developing tissues (e.g., root tips and young leaves), peripheral tissues (e.g., of flower stalks) and reproductive organs (flower buds and flowers). Within flowers, parts that usually attract insect visitors (e.g., labellum with colorful crests as well as column and pollinia) show the highest alkaloid levels. Tracer feeding experiments with (14)C-labeled putrecine revealed that in rosette plants the aerial roots were the sites of phalaenopsine biosynthesis. However active biosynthesis was only observed in roots still attached to the plant but not in excised roots. There is a slow but substantial translocation of newly synthesized alkaloid from the roots to other plant organs. A long-term tracer experiment revealed that phalaenopsine shows neither turnover nor degradation. The results are discussed in the context of a polyphyletic molecular origin of the biosynthetic pathways of pyrrolizidine alkaloids in various scattered angiosperm taxa. The ecological role of the so called non-toxic 1,2-saturated pyrrolizidine alkaloids is discussed in comparison to the pro-toxic 1,2-unsaturated pyrrolizidine alkaloids. Evidence from the plant-insect interphase is presented indicating a substantial role of the 1,2-saturated alkaloids in plant and insect defense.     

Taxonomy and biogeography of Diapensia (Diapensiaceae) based on chloroplast genome data

Diapensia L. is the second largest genus of Diapensiaceae. The taxonomic treatment within Diapensia and relationships within Diapensiaceae have been disputed. Chloroplast genome sequence data have proved to be useful for plant phylogenetic analyses and species delimitation. In this study, we de novo sequenced and assembled 22 chloroplast genomes of 15 species of Diapensiaceae, including all recognized species of Diapensia with multiple samples. A super‐matrix containing a total of 107 genes and 18 taxa was constructed for phylogenetic analyses to resolve phylogenetic relationships among genera of the family and within Diapensia. The resulting phylogenetic tree showed the following strongly supported relationships: (Galax, (Pyxidanthera, (Berneuxia, ((Schizocodon, Diapensia), and Shortia s.s.)))). The dated phylogeny and reconstructed lineage‐through‐time plot for the family indicated rapid diversification in the Neogene and an acceleration of diversification rate after c. 8 Ma. Biogeographic analysis suggested that Diapensia originated in the Northeast Asian mountains approximately 6.06 Ma, followed by northward dispersal to the Arctic and southwestward dispersal to the Himalaya–Hengduan Mountains. Phylogenetic relationships within Diapensia were well resolved. Based on the phylogenetic results, we proposed to reinstate the species status of Diapensia bulleyana Forrest ex Diels, and raised D. purpurea f. albida to the species rank (D. albida [W. E. Evans] J. F. Ye comb. & stat. nov.). The distribution ranges of all species delineated based on the phylogenetic results were revised accordingly based on specimen occurrences. Our study adds new examples for the power of plastid genome data for resolving phylogenetic relationships and clarifying taxonomic disputes among closely related species.     

Chloroplast genome (cpDNA) of Cycas taitungensis and 56 cp protein-coding genes of Gnetum parvifolium: insights into cpDNA evolution and phylogeny of extant seed plants

Phylogenetic relationships among the 5 groups of extant seed plants are presently unsettled. To reexamine this long-standing debate, we determine the complete chloroplast genome (cpDNA) of Cycas taitungensis and 56 protein-coding genes encoded in the cpDNA of Gnetum parvifolium. The cpDNA of Cycas is a circular molecule of 163,403 bp with 2 typical large inverted repeats (IRs) of 25,074 bp each. We inferred phylogenetic relationships among major seed plant lineages using concatenated 56 protein-coding genes in 37 land plants. Phylogenies, generated by the use of 3 independent methods, provide concordant and robust support for the monophylies of extant seed plants, gymnosperms, and angiosperms. Within the modern gymnosperms are 2 highly supported sister clades: Cycas-Ginkgo and Gnetum-Pinus. This result agrees with both the "gnetifer" and "gnepines" hypotheses. The sister relationships in Cycas-Ginkgo and Gnetum-Pinus clades are further reinforced by cpDNA structural evidence. Branch lengths of Cycas-Ginkgo and Gnetum were consistently the shortest and the longest, respectively, in all separate analyses. However, the Gnetum relative rate test revealed this tendency only for the 3rd codon positions and the transversional sites of the first 2 codon positions. A PsitufA located between psbE and petL genes is here first detected in Anthoceros (a hornwort), cycads, and Ginkgo. We demonstrate that the PsitufA is a footprint descended from the chloroplast tufA of green algae. The duplication of ycf2 genes and their shift into IRs should have taken place at least in the common ancestor of seed plants more than 300 MYA, and the tRNAPro-GGG gene was lost from the angiosperm lineage at least 150 MYA. Additionally, from cpDNA structural comparison, we propose an alternative model for the loss of large IR regions in black pine. More cpDNA data from non-Pinaceae conifers are necessary to justify whether the gnetifer or gnepines hypothesis is valid and to generate solid structural evidence for the monophyly of extant gymnosperms.     

The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications

We determined the nucleotide sequence of the complete mitochondrial genome of the nematode species Anisakis simplex. The genome is circular, 13,916 bp in size and conforms to the general characteristics of nematode mitochondrial DNAs. The gene arrangement of A. simplex is the same as that of Ascaris suum and almost identical to those of rhabditid species with a minor exception concerning the relative position of the AT-rich and non-coding regions and radically different from those of spirurid species. Along with comparisons of gene arrangement, phylogenetic analyses (maximum parsimony, neighbour joining and maximum likelihood methods) based on concatenated amino acid sequences of 12 protein-coding genes from 13 nematode species provided strong support for the sister-group relationship between Ascaridida and Rhabditida. The Shimodaira-Hasegawa and Templeton's tests both rejected the alternative hypothesis of a closer relationship between Ascaridida and Spirurida. These results contradicted the traditional view of nematode classification and a recent molecular phylogenetic study of 18S rDNA data that assigned Ascaridida and Spirurida as being a sister-group. Mapping of gene arrangement across the phylogenetic tree lead to the assumption that the conserved gene arrangement found in Ascaridida-Rhabditida members might have been acquired after the most recent common ancestor of ascaridid/rhabditid members branched off from the basal stock of the rhabditid lineage.     

An updated genus‐wide phylogenetic analysis of Arisaema (Araceae) with reference to sections

Arisaema has a long and complicated taxonomic history regarding its infrageneric classification. In the latest system, 14 sections were recognized, based on an unpublished, tentative phylogenetic analysis; in addition the type species and nomenclatural priority for each section were confirmed. Here, we present an updated, genus‐wide phylogenetic analysis, based on four plastid non‐coding regions (3′trnL–trnF, rpl20–5′rps12, psbB–psbH and rpoC2–rps2) for > 150 accessions. The maximum parsimony and maximum likelihood phylogenetic analyses identified eight major clades and one branch with unique sequence variation, although the relationships were unclear due to a polytomy and weak support. In the phylogenetic trees, most of the sections proposed in the latest system were distinct and corresponded to the major clades, but some sections are not monophyletic. On the basis of the phylogenetic relationship: (1) A. schimperianum is treated as a member of section Arisaema, rather than section Tenuipistillata or section Sinarisaema; and (2) section Fimbriata, which was synonymised into section Attenuata in the latest system but is morphologically distinct from the other species, is redefined as a monotypic section. In conclusion, we recognize 15 sections of Arisaema and species‐level classifications are discussed in a phylogenetic context.     

Complete chloroplast genomes of Liliaceae (s.l.) species: comparative genomic and phylogenetic analyses

We first report the complete chloroplast (cp) genome of Fritillaria taipaiensis and determine its characteristics, sequence divergence and phylogenetic relationships by comparing it with complete cp genomes of Liliaceae s.l. (including e.g. Nartheciaceae, Amaryllidaceae and Asparagaceae) species obtained from NCBI Genbank. We show that the ycf1, ycf15 and infA genes have become pseudogenes or are lost in some of the seventeen Liliaceae species, and that dispersed repeats are prevailing among the four types of repeats (dispersed, palindromic, complement and tandem repeats). The number of simple sequence repeats ranged from 53 to 84 in the seventeen species, with mononucleotide repeats being the most abundant, followed by dinucleotides. A total of nine genes with positive selection sites were identified (atpB, atpE, ndhF, ndhH, petB, rpl2, rpl20, rpl22 and ycf2). Furthermore, we examined 19 mutational hotspot regions, including three coding regions (rps16, infA and rpl22) and sixteen non-coding regions. A phylogenetic analysis of the complete cp genomes and protein-coding sequences showed that Fritillaria is most closely related to Lilium. Moreover, Asparagus and Polygonatum, Hosta and Yucca are closely related to the Liliaceae. These results will contribute to further study of evolutionary patterns and phylogenetic relationships in Liliaceae s.l.