Complete plastid genome sequence of Primula sinensis (Primulaceae): structure comparison, sequence variation and evidence for accD transfer to nucleus

Comparison of the Complete Chloroplast Genomes of Astilbe: Two Korean Endemic Plant Species

BACKGROUND

Astilbe, consisting of about 18 species, is distributed throughout East Asia and Northeastern America, and most Astilbe species are widely cultivated as ornamental plants. A total of four species of Astilbe have been confirmed to be distributed throughout Korea, two of which are endemic to Korea.

METHODS

In this study, we sequenced and assembled the complete chloroplast genomes of two endemic Korean plants using Illumina sequencing technology, identified simple sequence repeats (SSRs) and repetitive sequences, and compared them with three previously reported chloroplast genomes.

RESULTS

The chloroplast genomes of the two species were 156,968 and 57,142 bp in length and had a four-part circular structure. They consisted of a large single-copy region of 87,223 and 87,272 bp and a small single-copy region of 18,167 and 18,138 bp, separated by a pair of inverted repeats (IRa and IRb, 25,789 and 25,866 bp). The genomes contained 130 genes, 49 SSRs, and 49 long repetitive sequences. Comparative analysis with the chloroplast genomes of five Astilbe species indicated that A. uljinensis was closely related to A. chinensis and A. taquetii to A. koreana.

CONCLUSIONS

This study provides valuable references for the identification of two endemic Korean Astilbe species and contributes to a deeper understanding of the phylogeny and evolution of the genus Astilbe.

The complete chloroplast genome sequence of Primula medogensis (Primulaceae) and its phylogeny

Primula medogensis W.B Ju, B. Xu & X.F. Gao 2023, a new species categorized under P. sect. Cordifoliae, was officially described in 2023. Given its recent classification, the genetic resources for this species are currently very limited. Here, we sequenced and assembled the first complete chloroplast genome of P. medogensis using Illumina sequencing technology. The complete chloroplast genome of P. medogensis is 151,486 bp in length, exhibiting a typical quadripartite structure. It consists of a large single-copy region (LSC; 83,407 bp) and a small single-copy region (SSC;17675 bp), separated by a pair of inverted repeat regions (IRs; 25202 bp). A total of 131 genes were annotated, including 86 protein-coding, 37 tRNA, and eight rRNA genes. The overall GC content was 37.1%. Phylogenetic analysis of 59 Primula species revealed a close relationship between P. medogensis and P. calliantha subsp. bryophila.

Chloroplast genomes of Simarouba Aubl., molecular evolution and comparative analyses within Sapindales

Simarouba, a neotropical genus in the family Simaroubaceae, currently lacks comprehensive genomic data in existing databases. This study aims to fill this gap by providing genomic resources for three Simarouba species, S. amara, S. versicolor, and S. glauca. It also aims to perform comparative molecular evolutionary analyses in relation to other species within the order Sapindales. The analysis of these three Simarouba species revealed the presence of the typical quadripartite structure expected in plastomes. However, some pseudogenization events were identified in the psbC, infA, rpl22, and ycf1 genes. In particular, the CDS of the psbC gene in S. amara was reduced from 1422 bp to 584 bp due to a premature stop codon. Nucleotide diversity data pointed to gene and intergenic regions as promising candidates for species and family discrimination within the group, specifically matK, ycf1, ndhF, rpl32, petA-psbJ, and trnS-trnG. Selection signal analyses showed strong evidence for positive selection on the rpl23 gene. Phylogenetic analyses indicated that S. versicolor and S. glauca have a closer phylogenetic relationship than S. amara. We provide chloroplast genomes of three Simaruba species and use them to elucidate plastome evolution, highlight the presence of pseudogenization, and identify potential DNA barcode regions.

Chloroplast genome of Justicia procumbens: genomic features, comparative analysis, and phylogenetic relationships among Justicieae species

Justicia procumbens L. is a traditional medicinal plant that is widely distributed in China. However, little is known about the genetic diversity and evolution of this genus, and no genomic studies have been carried out on J. procumbens previously. In this study, we aimed to assemble and annotate the first complete chloroplast genome (cpDNA) of J. procumbens and compare it with all previously published cpDNAs within the tribe Justicieae. Genome structure and comparative and phylogenetic analyses were performed. The 150,454 bp-long J. procumbens cpDNA has a circular and quadripartite structure consisting of a large single copy, a small single copy, and two inverted repeat regions. It contains 133 genes, of which 88 are protein-coding genes, 37 are tRNA genes, and eight are rRNA genes. Twenty-four simple sequence repeats (SSRs) and 81 repeat sequences were identified. Comparative analyses with other Justicieae species revealed that the non-coding regions of J. procumbens cpDNA showed greater variation than did the coding regions. Moreover, phylogenetic analysis based on 14 cpDNA sequences from Justicieae species showed that J. procumbens and J. flava were most closely related. This study provides valuable genetic information to support further research on the genetic diversity and evolutionary development of the tribe Justicieae.

The complete chloroplast genome of Primula amethystina subsp. argutidens (Primulaceae)

Primula amethystina subsp. argutidens (Franchet) W. W. Smith & H. R. Fletcher (1942) is a blooming plant of the family Primulaceae. Here, we sequenced, assembled, and annotated the complete chloroplast (cp) genome of P. amethystina subsp. argutidens. The cp genome of P. amethystina subsp. argutidens is 151,560 bp in length with a GC content of 37%. The assembled genome has a typical quadripartite structure, containing a large single-copy (LSC) region of 83,516 bp, a small single-copy (SSC) region of 17,692 bp, and a pair of inverted repeat (IR) regions of 25,176 bp. The cp genome contains 115 unique genes, including 81 protein-coding genes, four rRNA genes, and 30 tRNA genes. Phylogenetic analysis showed that P. amethystina subsp. argutidens was closely related to P. amethystina.

The complete chloroplast genome of Primula vialii (Primulaceae), an ornamental plant

Primula vialii Delavay ex Franch. (1905) is an alpine species with an ornamental value. In this study, we sequenced, assembled, and annotated the chloroplast genome of P. vialii. The results showed that it was a double-stranded, closed circular DNA with 154,897 bp in length, comprising a small single-copy (SSC) region of 17,766 bp, a large single-copy (LSC) region of 85,379 bp and a pair of inverted repeat (IR) regions of 25,876 bp. A total of 113 unique genes were annotated, including 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. The phylogenetic analysis revealed that P. vialii is closely related to Primula flaccida. The cp genomic data will be useful for systematics and evolutionary studies of Primula.

ORTHOSKIM: in silico sequence capture from genomic and transcriptomic libraries for phylogenomic and barcoding applications

Low-coverage whole genome shotgun sequencing (or genome skimming) has emerged as a cost-effective method for acquiring genomic data in non-model organisms. This method provides sequence information on chloroplast genome (cpDNA), mitochondrial genome (mtDNA) and nuclear ribosomal regions (rDNA), which are over-represented within cells. However, numerous bioinformatic challenges remain to accurately and rapidly obtain such data in organisms with complex genomic structures and rearrangements, in particular for mtDNA in plants or for cpDNA in some plant families. Here we introduce the pipeline ORTHOSKIM, which performs in silico capture of targeted sequences from genomic and transcriptomic libraries without assembling whole organelle genomes. ORTHOSKIM proceeds in three steps: 1) global sequence assembly, 2) mapping against reference sequences, and 3) target sequence extraction; importantly it also includes a range of quality control tests. Different modes are implemented to capture both coding and non-coding regions of cpDNA, mtDNA and rDNA sequences, along with predefined nuclear sequences (e.g. ultra-conserved elements) or collections of single-copy ortholog genes. Moreover, aligned DNA matrices are produced for phylogenetic reconstructions, by performing multiple alignments of the captured sequences. While ORTHOSKIM is suitable for any eukaryote, a case study is presented here, using 114 genome-skimming libraries and 4 RNAseq libraries obtained for two plant families, Primulaceae and Ericaceae, the latter being a well-known problematic family for cpDNA assemblies. ORTHOSKIM recovered with high success rates cpDNA, mtDNA and rDNA sequences, well suited to accurately infer evolutionary relationships within these families. ORTHOSKIM is released under a GPL-3 license and is available at: https://github.com/cpouchon/ORTHOSKIM.

Mahonia vs. Berberis Unloaded: Generic Delimitation and Infrafamilial Classification of Berberidaceae Based on Plastid Phylogenomics

The early-diverging eudicot family Berberidaceae is composed of a morphologically diverse assemblage of disjunctly distributed genera long praised for their great horticultural and medicinal values. However, despite century-long studies, generic delimitation of Berberidaceae remains controversial and its tribal classification has never been formally proposed under a rigorous phylogenetic context. Currently, the number of accepted genera in Berberidaceae ranges consecutively from 13 to 19, depending on whether to define Berberis, Jeffersonia, and Podophyllum broadly, or to segregate these three genera further and recognize Alloberberis, Mahonia, and Moranothamnus, Plagiorhegma, and Dysosma, Diphylleia, and Sinopodophyllum, respectively. To resolve Berberidaceae's taxonomic disputes, we newly assembled 23 plastomes and, together with 85 plastomes from the GenBank, completed the generic sampling of the family. With 4 problematic and 14 redundant plastome sequences excluded, robust phylogenomic relationships were reconstructed based on 93 plastomes representing all 19 genera of Berberidaceae and three outgroups. Maximum likelihood phylogenomic relationships corroborated with divergence time estimation support the recognition of three subfamilies Berberidoideae, Nandinoideae, and Podophylloideae, with tribes Berberideae and Ranzanieae, Leonticeae and Nandineae, and Podophylleae, Achlydeae, Bongardieae tr. nov., Epimedieae, and Jeffersonieae tr. nov. in the former three subfamilies, respectively. By applying specifically stated criteria, our phylogenomic data also support the classification of 19 genera, recognizing Alloberberis, Mahonia, and Moranothamnus, Plagiorhegma, and Diphylleia, Dysosma, and Sinopodophyllum that are morphologically and evolutionarily distinct from Berberis, Jeffersonia, and Podophyllum, respectively. Comparison of plastome structures across Berberidaceae confirms inverted repeat expansion in the tribe Berberideae and reveals substantial length variation in accD gene caused by repeated sequences in Berberidoideae. Comparison of plastome tree with previous studies and nuclear ribosomal DNA (nrDNA) phylogeny also reveals considerable conflicts at different phylogenetic levels, suggesting that incomplete lineage sorting and/or hybridization had occurred throughout the evolutionary history of Berberidaceae and that Alloberberis and Moranothamnus could have resulted from reciprocal hybridization between Berberis and Mahonia in ancient times prior to the radiations of the latter two genera.

Highly variable chloroplast genome from two endangered Papaveraceae lithophytes Corydalis tomentella and Corydalis saxicola

The increasingly wide application of chloroplast (cp) genome super-barcode in taxonomy and the recent breakthrough in cp genetic engineering make the development of new cp gene resources urgent and significant. Corydalis is recognized as the most genotypes complicated and taxonomically challenging plant taxa in Papaveraceae. However, there currently are few reports about cp genomes of the genus Corydalis. In this study, we sequenced four complete cp genomes of two endangered lithophytes Corydalis saxicola and Corydalis tomentella in Corydalis, conducted a comparison of these cp genomes among each other as well as with others of Papaveraceae. The cp genomes have a large genome size of 189,029-190,247 bp, possessing a quadripartite structure and with two highly expanded inverted repeat (IR) regions (length: 41,955-42,350 bp). Comparison between the cp genomes of C. tomentella, C. saxicola, and Papaveraceae species, five NADH dehydrogenase-like genes (ndhF, ndhD, ndhL, ndhG, and ndhE) with psaC, rpl32, ccsA, and trnL-UAG normally located in the SSC region have migrated to IRs, resulting in IR expansion and gene duplication. An up to 9 kb inversion involving five genes (rpl23, ycf2, ycf15, trnI-CAU, and trnL-CAA) was found within IR regions. The accD gene was found to be absent and the ycf1 gene has shifted from the IR/SSC border to the SSC region as a single copy. Phylogenetic analysis based on the sequences of common CDS showed that the genus Corydalis is quite distantly related to the other genera of Papaveraceae, it provided a new clue for recent advocacy to establish a separate Fumariaceae family. Our results revealed one special cp genome structure in Papaveraceae, provided a useful resources for classification of the genus Corydalis, and will be valuable for understanding Papaveraceae evolutionary relationships.

The complete chloroplast genome of a distylous-homostylous species, Primula homogama (Primulaceae)

Primula homogama F. H. Chen & C. M. Hu (Primulaceae) is endemic to the Emei Mountain of China. In this study, we characterized the complete chloroplast genome of P. homogama based on next-generation sequencing (NGS). The complete chloroplast genome of P. homogama was 154,677 bp in size with a typical quadripartite structure, containing a large single-copy (LSC) region of 85,299 bp and a small single-copy (SSC) region of 17,816 bp. These two regions were separated by a pair of inverted repeat regions (IRs), each of 25,781 bp. A total of 130 functional genes were encoded, consisted of 86 protein-coding genes (PCG), 36 tRNA genes, and eight ribosomal RNA (rRNA) genes.

Determination of the optimum extraction regime of reducing compounds and flavonoids of Primula denticulata Smith leaves by a dispersion analysis

Herbal medicines are widely used in the complex treatment of various diseases. Therefore, theoretical and practical interest is the in-depth study of drumstick primrose (Primula denticulata Smith). The study aimed to determine the optimal extraction mode of flavonoids and reducing compounds of drumstick primrose leaves. The concentration of ethanol, the ratio of raw materials and extractant, and extraction method were studied by dispersion analysis. This allowed reducing the number of experiments from 64 to 16. To obtain the alcohol extract of drumstick primrose leaves with the highest content of reducing compounds and flavonoids, found that maceration is the optimal method of extraction, the ratio of raw materials to extractant should be 1 to 5 and 40% ethanol is the most appropriate extractant.

The complete chloroplast genome sequences of five pinnate-leaved Primula species and phylogenetic analyses

The six pinnate-leaved species are a very particular group in the genus Primula. In the present paper, we sequenced, assembled and annotated the chloroplast genomes of five of them (P. cicutarrifolia, P. hubeiensis, P. jiugongshanensis, P. merrilliana, P. ranunculoides). The five chloroplast genomes ranged from ~ 150 to 152 kb, containing 113 genes (four ribosomal RNA genes, 29 tRNA genes and 80 protein-coding genes). The six pinnate-leaved species exhibited synteny of gene order and possessed similar IR boundary regions in chloroplast genomes. The gene accD was pseudogenized in P. filchnerae. In the chloroplast genomes of the six pinnate-leaved Primula species, SSRs, repeating sequences and divergence hotspots were identified; ycf1 and trnH-psbA were the most variable markers among CDSs and noncoding sequences, respectively. Phylogenetic analyses showed that the six Primula species were separated into two distant clades: one was formed by P. filchnerae and P. sinensis and the other clade was consisting of two subclades, one formed by P. hubeiensis and P. ranunculoides, the other by P. merrilliana, P. cicutarrifolia and P. jiugongshanensis. P. hubeiensis was closely related with P. ranunculoides and therefore it should be placed into Sect. Ranunculoides. P. cicutarrifolia did not group first with P. ranunculoides but with P. merrilliana, although the former two were once united in one species, our results supported the separation of P. ranunculoides from P. cicutarrifolia as one distinct species.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae

Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

Complete chloroplast genome sequence of Barleria prionitis, comparative chloroplast genomics and phylogenetic relationships among Acanthoideae

BACKGROUND

The plastome of medicinal and endangered species in Kingdom of Saudi Arabia, Barleria prionitis was sequenced. The plastome was compared with that of seven Acanthoideae species in order to describe the plastome, spot the microsatellite, assess the dissimilarities within the sampled plastomes and to infer their phylogenetic relationships.

RESULTS

The plastome of B. prionitis was 152,217 bp in length with Guanine-Cytosine and Adenine-Thymine content of 38.3 and 61.7% respectively. It is circular and quadripartite in structure and constitute of a large single copy (LSC, 83, 772 bp), small single copy (SSC, 17, 803 bp) and a pair of inverted repeat (IRa and IRb 25, 321 bp each). 131 genes were identified in the plastome out of which 113 are unique and 18 were repeated in IR region. The genome consists of 4 rRNA, 30 tRNA and 80 protein-coding genes. The analysis of long repeat showed all types of repeats were present in the plastome and palindromic has the highest frequency. A total number of 98 SSR were also identified of which mostly were mononucleotide Adenine-Thymine and are located at the non coding regions. Comparative genomic analysis among the plastomes revealed that the pair of the inverted repeat is more conserved than the single copy region. In addition high variation is observed in the intergenic spacer region than the coding region. The genes, ycf1and ndhF and are located at the border junction of the small single copy region and IRb region of all the plastome. The analysis of sequence divergence in the protein coding genes indicates that the following genes undergo positive selection (atpF, petD, psbZ, rpl20, petB, rpl16, rps16, rpoC, rps7, rpl32 and ycf3). Phylogenetic analysis indicated sister relationship between Ruellieae and Justcieae. In addition, Barleria, Justicia and Ruellia are paraphyletic, suggesting that Justiceae, Ruellieae, Andrographideae and Barlerieae should be treated as tribes.

CONCLUSIONS

This study sequenced and assembled the first plastome of the taxon Barleria and reported the basics resources for evolutionary studies of B. prionitis and tools for phylogenetic relationship studies within the core Acanthaceae.

-Complete plastid genome sequences of two species of the Neotropical genus Brunellia (Brunelliaceae)

Here we present the first two complete plastid genomes for Brunelliaceae, a Neotropical family with a single genus, Brunellia. We surveyed the entire plastid genome in order to find variable cpDNA regions for further phylogenetic analyses across the family. We sampled morphologically different species, B. antioquensis and B. trianae, and found that the plastid genomes are 157,685 and 157,775 bp in length and display the typical quadripartite structure found in angiosperms. Despite the clear morphological distinction between both species, the molecular data show a very low level of divergence. The amount of nucleotide substitutions per site is one of the lowest reported to date among published congeneric studies (π = 0.00025). The plastid genomes have gene order and content coincident with other COM (Celastrales, Oxalidales, Malpighiales) relatives. Phylogenetic analyses of selected superrosid representatives show high bootstrap support for the ((C,M)O) topology. The N-fixing clade appears as the sister group of the COM clade and Zygophyllales as the sister to the rest of the fabids group.

Chloroplast Genomes and Comparative Analyses among Thirteen Taxa within Myrsinaceae s.str. Clade (Myrsinoideae, Primulaceae)

The Myrsinaceae s.str. clade is a tropical woody representative in Myrsinoideae of Primulaceae and has ca. 1300 species. The generic limits and alignments of this clade are unclear due to the limited number of genetic markers and/or taxon samplings in previous studies. Here, the chloroplast (cp) genomes of 13 taxa within the Myrsinaceae s.str. clade are sequenced and characterized. These cp genomes are typical quadripartite circle molecules and are highly conserved in size and gene content. Three pseudogenes are identified, of which ycf15 is totally absent from five taxa. Noncoding and large single copy region (LSC) exhibit higher levels of nucleotide diversity (Pi) than other regions. A total of ten hotspot fragments and 796 chloroplast simple sequence repeats (SSR) loci are found across all cp genomes. The results of phylogenetic analysis support the notion that the monophyletic Myrsinaceae s.str. clade has two subclades. Non-synonymous substitution rates (d_N) are higher in housekeeping (HK) genes than photosynthetic (PS) genes, but both groups have a nearly identical synonymous substitution rate (d_S). The results indicate that the PS genes are under stronger functional constraints compared with the HK genes. Overall, the study provides hypervariable molecular markers for phylogenetic reconstruction and contributes to a better understanding of plastid gene evolution in Myrsinaceae s.str. clade.

Rhopalocnemis phalloides has one of the most reduced and mutated plastid genomes known

Although most plant species are photosynthetic, several hundred species have lost the ability to photosynthesize and instead obtain nutrients via various types of heterotrophic feeding. Their plastid genomes markedly differ from the plastid genomes of photosynthetic plants. In this work, we describe the sequenced plastid genome of the heterotrophic plant Rhopalocnemis phalloides, which belongs to the family Balanophoraceae and feeds by parasitizing other plants. The genome is highly reduced (18,622 base pairs vs. approximately 150 kbp in autotrophic plants) and possesses an extraordinarily high AT content, 86.8%, which is inferior only to AT contents of plastid genomes of Balanophora, a genus from the same family. The gene content of this genome is quite typical of heterotrophic plants, with all of the genes related to photosynthesis having been lost. The remaining genes are notably distorted by a high mutation rate and the aforementioned AT content. The high AT content has led to sequence convergence between some of the remaining genes and their homologs from AT-rich plastid genomes of protists. Overall, the plastid genome of R. phalloides is one of the most unusual plastid genomes known.

The complete chloroplast genome of an endangered endemic herb species in China, Primula filchnerae (Primulaceae)

Primula filchnerae (Primulaceae) is an endangered endemic herb species in China. In this study, we characterized the complete chloroplast genome of P. filchnerae based on next generation sequencing (NGS). The chloroplast genome of P. filchnerae was 151,547 bp in size, containing a large single-copy (LSC) region of 82, 662 bp and a small single-copy (SSC) region of 17,749 bp. These two regions were separated by a pair of inverted repeat regions (IRs), each of 25,568 bp. A total of 130 functional genes were encoded, consisted of 86 protein-coding genes, 36 tRNA genes, and eight rRNA genes.

The plastid genome and its implications in barcoding specific-chemotypes of the medicinal herb Pogostemon cablin in China

Pogostemon cablin (Blanco) Benth. (Patchouli) is not only an important essential oil plant, but also a valuable medicinal plant in China. P. cablin in China can be divided into three cultivars (Shipai, Gaoyao, and Hainan) and two chemotypes (pogostone-type and patchoulol-type). The pogostone-type and patchoulol-type are, respectively, used for medicinals and perfumes. In this study, we sequenced and characterized the plastid genomes for all three Chinese cultivars and aimed to develop a chemotype-specific barcode for future quality control. The plastid genomes of P. cablin cultivars ranged from 152,461 to 152,462 bp in length and comprise 114 genes including 80 protein coding genes, 30 tRNA genes, and four rRNA genes. Phylogenetic analyses suggested that P. cablin cultivars clustered with the other two Pogostemon species with strong support. Although extremely conserved in P. cablin plastid genomes, 58 cpSSRs were filtered out among the three cultivars. One single variable locus, cpSSR, was discovered. The cpSSR genotypes successfully matched the chemotypes of Chinese patchouli, which was further supported by PCR-based Sanger sequences in more Chinese patchouli samples. The barcode developed in this study is thought to be a simple and reliable quality control method for Chinese P. cablin on the market.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

The Complete Chloroplast Genomes of Echinacanthus Species (Acanthaceae): Phylogenetic Relationships, Adaptive Evolution, and Screening of Molecular Markers

Among the four species of Echinacanthus (Acanthaceae), one distributed in the West Himalayan region and three restricted to the Sino-Vietnamese karst region. Because of its ecological significance, molecular markers are necessary for proper assessment of its genetic diversity and phylogenetic relationships. Herein, the complete chloroplast genomes of four Echinacanthus species were determined for the first time. The results indicated that all the chloroplast genomes were mapped as a circular structure and each genomes included 113 unique genes, of which 80 were protein-coding, 29 were tRNAs, and 4 were rRNAs. However, the four cp genomes ranged from 151,333 to 152,672 bp in length. Comparison of the four cp genomes showed that the divergence level was greater between geographic groups. We also analyzed IR expansion or contraction in the four cp genomes and the fifth type of the large single copy/inverted repeat region in Lamiales was suggested. Furthermore, based on the analyses of comparison and nucleotide variability, six most divergent sequences (rrn16, ycf1, ndhA, rps16-trnQ-UUG, trnS-GCU-trnG-UCC, and psaA-ycf3) were identified. A total of 37-45 simple sequence repeats were discovered in the four species and 22 SSRs were identified as candidate effective molecular markers for detecting interspecies polymorphisms. These SSRs and hotspot regions could be used as potential molecular markers for future study. Phylogenetic analysis based on Bayesian and parsimony methods did not support the monophyly of Echinacanthus. The phylogenetic relationships among the four species were clearly resolved and the results supported the recognition of the Sino-Vietnamese Echinacanthus species as a new genus. Based on the protein sequence evolution analysis, 12 genes (rpl14, rpl16, rps4, rps15, rps18, rps19, psbK, psbN, ndhC, ndhJ, rpoB, and infA) were detected under positive selection in branch of Sino-Vietnamese Echinacanthus species. These genes will lead to understanding the adaptation of Echinacanthus species to karst environment. The study will help to resolve the phylogenetic relationship and understand the adaptive evolution of Echinacanthus. It will also provide genomic resources and potential markers suitable for future species identification and speciation studies of the genus.

Comparative Plastid Genomes of Primula Species: Sequence Divergence and Phylogenetic Relationships

Compared to traditional DNA markers, genome-scale datasets can provide mass information to effectively address historically difficult phylogenies. Primula is the largest genus in the family Primulaceae, with members distributed mainly throughout temperate and arctic areas of the Northern Hemisphere. The phylogenetic relationships among Primula taxa still maintain unresolved, mainly due to intra- and interspecific morphological variation, which was caused by frequent hybridization and introgression. In this study, we sequenced and assembled four complete plastid genomes (Primula handeliana, Primula woodwardii, Primula knuthiana, and Androsace laxa) by Illumina paired-end sequencing. A total of 10 Primula species (including 7 published plastid genomes) were analyzed to investigate the plastid genome sequence divergence and their inferences for the phylogeny of Primula. The 10 Primula plastid genomes were similar in terms of their gene content and order, GC content, and codon usage, but slightly different in the number of the repeat. Moderate sequence divergence was observed among Primula plastid genomes. Phylogenetic analysis strongly supported that Primula was monophyletic and more closely related to Androsace in the Primulaceae family. The phylogenetic relationships among the 10 Primula species showed that the placement of P. knuthiana–P. veris clade was uncertain in the phylogenetic tree. This study indicated that plastid genome data were highly effective to investigate the phylogeny.

Plastid phylogenomics resolves infrafamilial relationships of the Styracaceae and sheds light on the backbone relationships of the Ericales

Relationships among the genera of the small, woody family Styracaceae and among families of the large, diverse order Ericales have resisted complete resolution with sequences from one or a few genes. We used plastome sequencing to attempt to resolve the backbone relationships of Styracaceae and Ericales and to explore plastome structural evolution. Complete plastomes for 23 species are newly reported here, including 18 taxa of Styracaceae and five of Ericales (including species of Sapotaceae, Clethraceae, Symplocaceae, and Diapensiaceae). Combined with publicly available complete plastome data, this resulted in a data set of 60 plastomes, including 11 of the 12 genera of Styracaceae and 12 of 22 families of Ericales. Styracaceae plastomes were found to possess the quadripartite structure typical of angiosperms, with sizes ranging from 155 to 159 kb. Most of the plastomes were found to possess the full complement of typical angiosperm plastome genes. Unusual structural features were detected in plastomes of Alniphyllum and Bruinsmia, including the presence of a large 20-kb inversion (14 genes) in the Large Single-Copy region, the loss or pseudogenization of the clpP and accD genes in Bruinsmia, and the loss of the first exon of rps16 in B. styracoides. Likewise, the second intron from clpP was found to be lost in Alniphyllum and Huodendron. Phylogenomic analyses including all 79 plastid protein-coding genes provided improved resolution for relationships among the genera of Styracaceae and families of Ericales. Styracaceae was strongly supported as monophyletic, with Styrax, Huodendron, and a clade of Alniphyllum + Bruinsmia successively sister to the remainder of the family, all with strong support. All genera of Styracaceae were recovered as monophyletic, except for Halesia and Pterostyrax, which were each recovered as polyphyletic with strong support. Within Ericales, all families were recovered as monophyletic with strong support, with Balsaminaceae sister to remaining Ericales. Most relationships recovered in plastome analyses are congruent with previous analyses based on smaller data sets. Our results demonstrate the power of plastid phylogenomics to improve phylogenetic hypotheses among genera and families, and provide new insight into plastome evolution across Ericales.

Complete Chloroplast Genome Sequence of Decaisnea insignis: Genome Organization, Genomic Resources and Comparative Analysis

Decaisnea insignis is a wild resource plant and is used as an ornamental, medicinal, and fruit plant. High-throughput sequencing of chloroplast genomes has provided insight into the overall evolutionary dynamics of chloroplast genomes and has enhanced our understanding of the evolutionary relationships within plant families. In the present study, we sequenced the complete chloroplast genome of D. insignis and used the data to assess its genomic resources. The D. insignis chloroplast genome is 158,683 bp in length and includes a pair of inverted repeats of 26,167 bp that are separated by small and large single copy regions of 19,162 bp and 87,187 bp, respectively. We identified 83 simple sequence repeats and 18 pairs of large repeats. Most simple-sequence repeats were located in the noncoding sections of the large single-copy/small single-copy region and exhibited a high A/T content. The D. insignis chloroplast genome bias was skewed towards A/T on the basis of codon usage. A phylogenetic tree based on 82 protein-coding genes of 33 angiosperms showed that D. insignis was clustered with Akebia in Lardizabalaceae. Overall, the results of this study will contribute to better understanding the evolution, molecular biology and genetic improvement of D. insignis.

The complete chloroplast genome sequences of Lychnis wilfordii and Silene capitata and comparative analyses with other Caryophyllaceae genomes

The complete chloroplast genomes of Lychnis wilfordii and Silene capitata were determined and compared with ten previously reported Caryophyllaceae chloroplast genomes. The chloroplast genome sequences of L. wilfordii and S. capitata contain 152,320 bp and 150,224 bp, respectively. The gene contents and orders among 12 Caryophyllaceae species are consistent, but several microstructural changes have occurred. Expansion of the inverted repeat (IR) regions at the large single copy (LSC)/IRb and small single copy (SSC)/IR boundaries led to partial or entire gene duplications. Additionally, rearrangements of the LSC region were caused by gene inversions and/or transpositions. The 18 kb inversions, which occurred three times in different lineages of tribe Sileneae, were thought to be facilitated by the intermolecular duplicated sequences. Sequence analyses of the L. wilfordii and S. capitata genomes revealed 39 and 43 repeats, respectively, including forward, palindromic, and reverse repeats. In addition, a total of 67 and 56 simple sequence repeats were discovered in the L. wilfordii and S. capitata chloroplast genomes, respectively. Finally, we constructed phylogenetic trees of the 12 Caryophyllaceae species and two Amaranthaceae species based on 73 protein-coding genes using both maximum parsimony and likelihood methods.

Plastome Rearrangements in the "Adenocalymma-Neojobertia" Clade (Bignonieae, Bignoniaceae) and Its Phylogenetic Implications

The chloroplast is one of the most important organelles of plants. This organelle has a circular DNA with approximately 130 genes. The use of plastid genomic data in phylogenetic and evolutionary studies became possible with high-throughput sequencing methods, which allowed us to rapidly obtain complete genomes at a reasonable cost. Here, we use high-throughput sequencing to study the "Adenocalymma-Neojobertia" clade (Bignonieae, Bignoniaceae). More specifically, we use Hi-Seq Illumina technology to sequence 10 complete plastid genomes. Plastomes were assembled using selected plastid reads and de novo approach with SPAdes. The 10 assembled genomes were analyzed in a phylogenetic context using five different partition schemes: (1) 91 protein-coding genes ("coding"); (2) 76 introns and spacers with alignment manually edited ("non-coding edited"); (3) 76 non-coding regions with poorly aligned regions removed using T-Coffee ("non-coding filtered"); (4) 91 coding regions plus 76 non-coding regions edited ("coding + non-coding edited"); and, (5) 91 protein-coding regions plus the 76 filtered non-coding regions ("coding + non-coding filtered"). Fragmented regions were aligned using Mafft. Phylogenetic analyses were conducted using Maximum Likelihood (ML) and Bayesian Criteria (BC). The analyses of the individual plastomes consistently recovered an expansion of the Inverted Repeated (IRs) regions and a compression of the Small Single Copy (SSC) region. Major genomic translocations were observed at the Large Single Copy (LSC) and IRs. ML phylogenetic analyses of the individual datasets led to the same topology, with the exception of the analysis of the "non-coding filtered" dataset. Overall, relationships were strongly supported, with the highest support values obtained through the analysis of the "coding + non-coding edited" dataset. Four regions at the LSC, SSC, and IR were selected for primer development. The "Adenocalymma-Neojobertia" clade shows an unusual pattern of plastid structure variation, including four major genomic translocations. These rearrangements challenge the current view of conserved plastid genome architecture in terms of gene order. It also complicates both genomic assemblies using reference genomes and sequence alignments using whole plastomes. Therefore, strategies that employ de novo assemblies and manual evaluation of sequence alignments are required to prevent assembly and alignment errors.