Plastome phylogenomic analysis of Torreya (Taxaceae)

Torreya Arn., a small genus of Taxaceae, consists of six species occurring in North America and eastern Asia. Several phylogenetic studies have previously been undertaken to reveal relationships within this genus, although only a few DNA segments or species were used. In the present study, we sequenced five Torreya plastomes and combined these with two existing plastomes from the genus to investigate plastome evolution and phylogenetic relationships within Torreya. All sequenced Torreya plastomes shared the same complement of 82 protein‐coding genes, 4 ribosomal RNA genes, and 31 transfer RNA genes. Phylogenetic inference using a maximum likelihood framework consisted of an 82‐gene, 17‐taxon dataset, including all species of Torreya, resolved Torreya as a monophyletic clade. Strongly supported relationships within the genus include the position of the early diverging T. jackii Chun, the two sister pairs T. fargesii Franch.–T. nucifera (L.) Siebold & Zucc. and T. grandis Fortune ex Lindl.–T. californica Torr., and the monophyly of the clade including T. fargesii var. yunnanensis, T. fargesii, and T. nucifera. In addition to the inference of species relationships, divergence time estimation and biogeographical analysis were carried out. The diversification of Torreya was estimated to be approximately 8.9 Ma. Ancestral state reconstruction of the geographical area suggested China/eastern North America as the most likely ancestral region for the six extant Torreya species.     

Phylogenomics and biogeography of Wisteria: Implications on plastome evolution among inverted repeat-lacking clade (IRLC) legumes

The genus Wisteria (Fabaceae) is disjunctly distributed in eastern Asian and eastern North American temperate deciduous forests, and it is widely cultivated around the world as spectacular garden plants. It is a member of inverted repeat-lacking clade (IRLC). The IRLC Species are characterized by the loss of an IR region in their plastomes, which has long been of great interest. In this research, we report whole plastome sequences from all four Wisteria species and a Wisteriopsis japonica, combining these with existing data to explore phylogenetic relationships and biogeography of Wisteria, as well as plastome evolution of IRLC species. Phylogenetic analyses recognized a clade containing Glycyrrhiza–Wisteriopsis–Wisteria as sister to the remaining genera of IRLC. North American Wisteria frutescens and the three Asian species formed reciprocal clades, and Wisteria brachybotrys was sister to Wisteria floribunda and Wisteria sinensis. Wisteria may have originated in Japan near the boundary of the Oligocene and Miocene. The disappearance of Bering Land Bridge in the late Miocene might lead to the Eastern Asian–Eastern North American disjunction of Wisteria. Allopatric speciation of Wisteria between the Japanese archipelago and the Asian continent in the Quaternary increased the species richness of eastern Asia in comparison with eastern North America. Synonymous substitution rates (d_S) of protein-coding genes in the IRLC species were around 2-fold (SC genes) or 11-fold (IR genes) higher than those of non-IRLC species. For both SC and IR genes, herbaceous legumes have around 3-fold higher d_S than woody ones. Both loss of one IR region and herbaceous habit elevated substitution rates of the plastomes.     

Plastomes of nine hornbeams and phylogenetic implications

Poor phylogenetic resolution and inconsistency of gene trees are major complications when attempting to construct trees of life for various groups of organisms. In this study, we addressed these issues in analyses of the genus Carpinus (hornbeams) of the Betulaceae. We assembled and annotated the chloroplast (cp) genomes (plastomes) of nine hornbeams representing main clades previously distinguished in this genus. All nine plastomes are highly conserved, with four regions, and about 158–160 kb long, including 121–123 genes. Phylogenetic analyses of whole plastome sequences, noncoding sequences, and the well‐aligned coding genes resulted in high resolution of the sampled species in contrast to the failure based on a few cpDNA markers. Phylogenetic relationships in a few clades based only on the coding genes are slightly inconsistent with those based on the noncoding and total plastome datasets. Moreover, these plastome trees are highly incongruent with those based on bi‐parentally inherited internal transcribed spacer (ITS) sequence variations. Such high inconsistencies suggest widespread occurrence of incomplete lineage sorting and hybrid introgression during diversification of these hornbeams.     

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.     

青藏高原特有种-藏扇穗茅叶绿体基因组测序及序列分析

藏扇穗茅（Littledalea tibetica）是禾本科（Poaceae）雀麦族（Bromeae）中一个具有重要生态价值的多年生高山特有种,主要分布于青藏高原及其毗邻地区。本文采用基于第二代高通量测序平台的Illumina MiSeq技术,对青藏高原特有种—藏扇穗茅进行了叶绿体基因组测序,首次建立了雀麦族物种的标准测序流程;同时,以其近缘物种—黑麦草（Lolium perenne）的叶绿体基因组序列作为参考,组装获得它的叶绿体基因组序列。结果表明,藏扇穗茅叶绿体基因组序列全长136 852 bp,GC含量为38.5%,呈典型的四段式结构,其中大（LSC）、小（SSC）单拷贝区大小分别为80 970和12 876 bp,反向互补重复区（IR）大小为21 503 bp,共注释得到141个基因,包含95个蛋白编码基因、38个tRNA基因和8个rRNA基因,主要分布于大单拷贝区和小单拷贝区。同时,基于藏扇穗茅和其它30种禾本科植物叶绿体基因全序列构建的系统发育树显示,藏扇穗茅与早熟禾亚科中小麦族植物亲缘关系较近。     

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 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.     

Phylogenomic relationships and species identification of the olive genus Olea (Oleaceae)

The olive genus Olea includes c. 30–40 taxa in three subgenera (Olea, Tetrapilus, and Paniculatae) within the family Oleaceae. Historically, the Olea genus was classified into four groups that were overall well supported by reconstructed phylogenies, despite incomplete sampling of subgenus Tetrapilus and poor resolution within clades. These analyses also showed that the genus was not monophyletic. Reliable identification of Olea species is important for both their conservation and utilization of this economically important genus. In this study, we used phylogenomic data from genome skimming to resolve relationships within Olea and to identify molecular markers for species identification. We assembled the complete plastomes, and nrDNA of 26 individuals representing 13 species using next-generation sequencing and added 18 publicly available accessions of Olea. We also developed nuclear SNPs using the genome skimming data to infer the phylogenetic relationships of Olea. Large-scale phylogenomic analyses of 138 samples of tribe Oleeae supported the polyphyly of Olea, with Olea caudatilimba and Olea subgenus Tetrapilus not sharing their most recent common ancestor with the main Olea clade (subgenus Paniculatae and subgenus Olea). The interspecific phylogenetic resolution was poor owing to a possible rapid radiation. By comparing with the plastome data, we identified the markers ycf1b and psbE-petL as the best Olea-specific chloroplast DNA barcodes. Compared with universal barcodes, specific DNA barcodes and super-barcode exhibited higher discriminatory power. Our results demonstrated the power of phylogenomics to improve phylogenetic relationships of intricate groups and provided new insights into barcodes that allow for accurate identification of Olea species.     

Plastome characteristics of Cannabaceae

Cannabaceae is an economically important family that includes ten genera and ca.117 accepted species. To explore the structure and size variation of their plastomes,we sequenced ten plastomes representing all ten genera of Cannabaceae.Each plastome possessed the typical angiosperm quadripartite structure and contained a total of 128 genes.The Inverted Repeat (IR) regions in five plastomes had experienced small expansions (330-983 bp) into the Large Single-Copy (LSC) region.The plastome of Chaetachme aristata has experienced a 942-bp IR contraction and lost rpl22 and rps19 in its IRs.The substitution rates of rps19 and rpl22 decreased after they shifted from the LSC to IR.A 270-bp inversion was detected in the Parasponia rugosa plastome,which might have been mediated by 18-bp inverted repeats.Repeat sequences,simple sequence repeats,and nucleotide substitution rates varied among these plastomes. Molecular markers with more than 13% variable sites and 5% parsimony-informative sites were identified,which may be useful for further phylogenetic analysis and species identification.Our results show strong support for a sister relationship between Gironniera and Lozanell (BS=100).Celtis,Cannabis-Humulus,Chaetachme-Pteroceltis,and Trema-Parasponia formed a strongly supported clade,and their relationships were well resolved with strong support (BS=100).The availability of these ten plastomes provides valuable genetic information for accurately identifying species,clarifying taxonomy and reconstructing the intergeneric phylogeny of Cannabaceae.     

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.     

Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics

Here, we report a completely sequenced plastome using Illumina/Solexa sequencing-by-synthesis (SBS) technology. The plastome of Magnolia kwangsiensis Figlar & Noot. is 159?667 bp in length with a typical quadripartite structure: 88?030 bp large single-copy (LSC) and 18?669 bp small single-copy (SSC) regions, separated by two 26?484 bp inverted repeat (IR) regions. The overall predicted gene number is 129, among which 17 genes are duplicated in IR regions. The plastome of M. kwangsiensis is identical in its gene order to previously published plastomes of magnoliids. Furthermore, the C-to-U type RNA editing frequency of 114 seed plants is positively correlated with plastome GC content and plastome length, whereas plastome length is not correlated with GC content. A total of 16 potential putative barcoding or low taxonomic level phylogenetic study markers in Magnoliaceae were detected by comparing the coding and noncoding regions of the plastome of M. kwangsiensis with that of Liriodendron tulipifera L. At least eight markers might be applied not only to Magnoliaceae but also to other taxa. The 86 mononucleotide cpSSRs that distributed in single-copy noncoding regions are highly valuable to study population genetics and conservation genetics of this endangered rare species.     

Total duplication of the small single copy region in the angiosperm plastome: Rearrangement and inverted repeat instability in Asarum

Premise of the Study

As more plastomes are assembled, it is evident that rearrangements, losses, intergenic spacer expansion and contraction, and syntenic breaks within otherwise functioning plastids are more common than was thought previously, and such changes have developed independently in disparate lineages. However, to date, the magnoliids remain characterized by their highly conserved plastid genomes (plastomes).

Methods

Illumina HiSeq and MiSeq platforms were used to sequence the plastomes of Saruma henryi and those of representative species from each of the six taxonomic sections of Asarum. Sequenced plastomes were compared in a phylogenetic context provided by maximum likelihood and parsimony inferences made using an additional 18 publicly available plastomes from early‐diverging angiosperm lineages.

Key Results

In contrast to previously published magnoliid plastomes and the newly sequenced Saruma henryi plastome published here, Asarum plastomes have undergone extensive disruption and contain extremely lengthy AT‐repeat regions. The entirety of the small single copy region (SSC) of A. canadense and A. sieboldii var. sieboldii has been incorporated into the inverted repeat regions (IR), and the SSC of A. delavayi is only 14 bp long. All sampled Asarum plastomes share an inversion of a large portion of the large single copy region (LSC) such that trnE‐UUC is adjacent to the LSC‐IR boundary.

Conclusions

Plastome divergence in Asarum appears to be consistent with trends seen in highly rearranged plastomes of the monocots and eudicots. We propose that plastome instability in Asarum is due to repetitive motifs that serve as recombinatory substrates and reduce genome stability.     

The first complete plastome sequence of the basal asterid family Styracaceae (Ericales) reveals a large inversion

Plastome sequences are rich sources of information for resolving difficult phylogenetic relationships and provide genomic data for conservation studies. Here, the complete plastome sequence of Alniphyllum eberhardtii Guillaumin is reported, representing the first plastome of the basal asterid family Styracaceae (Ericales). The plastome is 155,384 bp in length and contains 79 protein-coding genes, 30 tRNA genes and 4 rRNA genes, totaling 113 unique genes with 19 genes in the inverted repeat region. Unusual features of the plastome include the presence a large 20-kb inversion in the Large Single-Copy region, the pseudogenization of the accD gene, and the loss of the second intron from clpP. The 20-kb inversion includes 14 genes and has not been previously reported in other Ericales plastomes. Thirty-nine plastid simple sequence repeats (SSRs) that may provide genetic resources for the conservation of this economically import timber plant are characterized. Phylogenetic results inferred from ML and MP analyses of 66 plastid genes and 26 taxa reveal that the Styracaceae are sister to a clade including Actinidiaceae and Ericaceae and suggest that complete plastomes are likely to be very helpful in resolving the basal relationships among Ericales families, which have resisted resolution in smaller phylogenetic data sets.     

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 Complete Chloroplast Genome of <Emphasis Type="Italic">Coix lacryma-jobi</Emphasis> and a Comparative Molecular Evolutionary Analysis of Plastomes in Cereals

Graminoid molecular evolution was investigated by chloroplast genome (plastome) scale analyses. A complete plastome from Coix lacryma-jobi (Poaceae) and a draft plastome from Joinvillea plicata (Joinvilleaceae) were sequenced and analyzed. The draft plastome included conserved protein-coding loci routinely analyzed in previous studies plus one additional locus of demonstrated phylogenetic utility. The methodological approach was to directly sequence overlapping amplicons from known plastome regions. Over 100 pairs of amplification and sequencing primers were designed and positioned to flank overlapping 1,200-base pair fragments around the entire plastome. Newly determined sequences were analyzed with published plastomes from representatives of Panicoideae, Ehrhartoideae, and Pooideae. Considerable variation was found for studies within the family and even within Andropogoneae. Readily interpreted mutation patterns were observed, such as small inversions in hairpin-loop regions and indels, which were common in intergenic spacers. Maximum or near-maximum bootstrap support was observed in all analyses resolving relationships between subfamilies. However, the addition of characters from noncoding regions increased the number of parsimony-informative characters and lengthened short internal branches (Andropogoneae), better defining intergeneric relationships. Thus, characters in complete plastomes can be used over a wide scope of phylogenetic studies.     

Phylogenomics and rapid diversification of the genus Eutrema on the Qinghai–Tibet Plateau and adjacent regions

Both geographic isolation and polyploidization are assumed to play an important role in driving species diversification. However, this is rarely illustrated through phylogenomic analyses. The genus Eutrema (Brassicaceae), which also includes the salt-resistant species, are distributed mainly in Asia with extensive species diversification in the Qinghai–Tibet Plateau (QTP) and adjacent regions. In this study, we revealed almost fully resolved backbone relationships of the genus with genome re-sequencing data for genomes of 168 individuals from 28 species. Phylogenetic analyses of both plastomes and single-copy nuclear genes from the whole genome recovered six well-supported clades with almost consistent relationships. The first two clades are mainly distributed in central China and central Asia, while the other four in the QTP and adjacent regions. All of them diversified within 12 million years. Within each clade, we recovered numerous conflicts in the interspecific relationships between nuclear and plastome phylogenies, likely suggesting hybridization and incomplete lineage sorting during species diversification. Our estimation of genome size and comparison of the number of the single-copy nuclear genes demonstrated frequent occurrences of polyploids in the genus. Except for an establishment of the backbone phylogeny, our phylogenomic analyses suggest that in addition to strong geographic isolation, polyploidization may have played an important role in species diversification of this genus.     

Plastomes of Mimosoideae: structural and size variation,sequence divergence,and phylogenetic implication

Plastomes of Fabaceae show both significant structural and size variation; however, most published plastomes are from subfamily Papilionoideae and only a few are from the other two subfamilies. In order to address the plastome structural and size variation of subfamily Mimosoideae, we integrated 11 newly sequenced plastomes from representing genera with three previously published ones. Each mimosoid plastome presented a typical quadripartite structure and contained 111 unique genes. Their inverted repeats (IRs) experienced multiple expansion/contraction; a ca. 13-kb IR expansion into small single copy (SSC) was detected in plastomes of a clade formed by tribe Ingeae and Acacia sensu stricto (s.s.), and a ca. 1.7-kb IR expansion into and a ca. 1.9-kb contraction out of large single copy (LSC) were found in Pithecellobium flexicaule and Acacia dealbata, respectively. Linear regression analysis showed decreased synonymous substitution rates of genes relocating from SSC into IR. A loss of both introns of clpP occurred in A. dealbata and Faidherbia albida, and a duplicated clpP copy was detected in A. dealbata. Furthermore, a 421-bp inversion that containing rps18 was found in A. dealbata. The size of mimosoid plastomes was found significantly affected by a IR-SC boundary shift, and also associated with repeat content. Plastome coding and noncoding regions with variable sequence divergence may supply valuable markers for molecular evolutionary and phylogenetic studies at different taxonomic levels. Plastid phylogenomics well resolved relationships among sampled mimosoid species.     

长柄扁桃叶绿体基因组比较

蔷薇科桃属植物长柄扁桃Amygdalus pedunculata Pall.是我国重点发展的12种木本油料作物之一。将自然生长在陕西省榆林市毛乌素沙地的长柄扁桃进行高通量测序、组装得到了GenBank登录号为MG602257的长柄扁桃叶绿体基因组。与NCBI中其他两个长柄扁桃叶绿体基因组 (MG869261和KY101153) 进行序列比对发现：长柄扁桃叶绿体基因组均呈现典型的4分区域结构;基因组序列GC含量均是36.8%;总长度MG602257相似文献