CPGAVAS2, an integrated plastome sequence annotator and analyzer

The complete chloroplast genome sequence of Paris polyphylla var. alba H.Li & R.J.Mitchell and its phylogenetic analysis

Paris polyphylla var. alba is a medicinal plant commonly used in the southwest of China. This study characterized the complete chloroplast (cp) genome sequence of P. polyphylla var. alba to investigate its phylogenetic relationship in Melanthiaceae. The cp genome of P. polyphylla var. alba is 165,079 base pairs (bp) in length with 36.96% G + C content. The cp genome is divided into (a) large single copy (LSC) (84,393 bp), (b) small single copy (SSC) (16,066 bp), and (c) two inverted regions (32,310 bp). The cp genome contains 135 genes, including 89 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis indicated that P. polyphylla var. alba is closest to P. polyphylla var. emeiensis, and Paris had a close relationship with Trillium in Melanthiaceae.

The genome of a thorny species: comparative genomic analysis among South and North American Cactaceae

The first South American cactus nuclear genome assembly associated with comparative genomic analyses provides insights into nuclear and plastid genomic features, such as size, transposable elements, and metabolic processes related to cactus development. Here, we assembled the partial genome, plastome, and transcriptome of Cereus fernambucensis (Cereeae, Cactaceae), a representative species of the South American core Cactoideae. We accessed other genomes and transcriptomes available for cactus species to compare the heterozygosity level, genome size, transposable elements, orthologous genes, and plastome structure. These estimates were obtained from the literature or using the same pipeline adopted for C. fermabucensis. In addition to the C. fernambucensis plastome, we also performed de novo plastome assembly of Pachycereus pringlei, Stenocereus thurberi, and Pereskia humboldtii based on the sequences available in public databases. We estimated a genome size of ~ 1.58 Gb for C. fernambucensis, the largest genome among the compared species. The genome heterozygosity was 0.88% in C. fernambucensis but ranged from 0.36 (Carnegiea gigantea) to 17.4% (Lophocereus schottii) in the other taxa. The genome lengths of the studied cacti are constituted by a high amount of transposable elements, ranging from ~ 57 to ~ 67%. Putative satellite DNAs are present in all species, excepting C. gigantea. The plastome of C. fernambucensis was ~ 104 kb, showing events of translocation, inversion, and gene loss. We observed a low number of shared unique orthologs, which may suggest gene duplication events and the simultaneous expression of paralogous genes. We recovered 37 genes that have undergone positive selection along the Cereus branch that are associated with different metabolic processes, such as improving photosynthesis during drought stress and nutrient absorption, which may be related to the adaptation to xeric areas of the Neotropics.

Chloroplast genome analysis of Chrysotila dentata

Chrysotila dentata is an ecologically important marine alga contributing to the coccolith formation. In this study, a complete chloroplast (cp DNA) genome of Chrysotila dentata was sequenced by using Illumina Hiseq and was analyzed with the help of a bioinformatics tool CPGAVAS2. The circular chloroplast genome of Chrysotila dentata has a size of 109,017 bp with two inverted repeats (IRs) regions (4513 bp each) which is a common feature in most land plants and algal species. The Chrysotila dentata cp genome consists of 61 identified protein-coding genes, 30 tRNA genes and 6 rRNAs with 21 microsatellites. The phylogenetic relationship with other select algal species revealed a close phylogeny of Chrysotila dentata with Phaeocystis antarctica. This is the first report of the cp genome analysis of genus Chrysotila and the results from this study will be helpful for understanding the genetic structure and function of chloroplast in other species of Chrysotila.

The complete chloroplast genome sequence of Prunus simonii ‘Weiwang’

Prunus simonii ‘Weiwang’ is an important economic fruit crops. In this study, we reported the complete chloroplast genome sequence of P. simonii ‘Weiwang’. The genome has a circular structure of 157,924 bp containing a large single-copy region (LSC) of 86,187 bp, a small copy region (SSC) of 19,031 bp, and two inverted repeats (IR) of 26,353 bp by each. It harbors 110 unique genes, including 78 protein-coding genes, 4 ribosomal RNA genes, and 28 transfer RNA genes. The phylogenomic analysis shows that Prunus simonii ‘Weiwang’ is clustered with Prunus salicina.

The complete chloroplast genome of Rhododendron kawakamii (Ericaceae)

Rhododendron kawakamii is endemic in Taiwan island and is a unique and epiphytic species. Here, we report its complete chloroplast genome. The length of the R. kawakamii chloroplast genome is 230,777 bp, with a large single-copy region of 146,155 bp, a small single-copy region of 72,082 bp, and a pair of inverted repeat regions (IRA) of 6,270 bp each. The genome contains 77 protein-coding genes, 29 transfer RNA genes, and four ribosomal RNA genes. In addition, the genome contains 81 simple sequence repeats. Phylogenetic analysis revealed that R. kawakamii is genetically related to R. datiandingense.

The complete plastid genome of Cotinus coggygria and phylogenetic analysis of the Anacardiaceae

Cotinus coggygria Scop. (Anacardiaceae) is an important ornamental tree with beautiful characteristics that is grown in China. In this study, the complete plastid genome of C. coggygria was sequenced and assembled. This genome was 158,843 bp in size and presented a typical tetrad structure, consisting of a large single-copy region (87,121 bp), a pair of inverted repeat regions (26,829 bp), and a small single-copy region (18,064 bp). A total of 134 genes were annotated, including 88 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. We observed a deletion that caused the loss of the rpl32 gene, and a small expansion of IR regions resulted in the trnH gene accessing IR regions; two copies were obtained. Phylogenetic analysis showed that C. coggygria was most closely related to Pistacia, with 100% bootstrap support within Anacardiaceae. In this study, we report the plastid genome of Cotinus species for the first time, which provides insight into the evolution of the plastid genome in Anacardiaceae and promotes the understanding of Cotinus plants.

The complete chloroplast genome of Agropyron pectinatum (M. Bieb.) P. Beauv

Agropyron pectinatum is a perennial forage widely cultivated in China, and it belongs to the Gramineous family. In this study, we assembled the complete chloroplast genome of A. pectinatum. The whole chloroplast genome of A. pectinatum is 135,041 bp in length, comprising a pair of inverted repeat (IR) regions (20,821 bp) that are separated by a large single copy (LSC) region (80,632 bp) and a small single copy (SSC) region (12,767 bp). The chloroplast genome of A. pectinatum contains 133 genes, and 87 of them are protein-coding genes, 38 are tRNA, and eight are rRNA genes. The chloroplast genome of A. pectinatum could provide valuable information for varieties identification and evolution of the Agropyron Gaertn.

Chloroplast Genome Sequencing, Comparative Analysis, and Discovery of Unique Cytoplasmic Variants in Pomegranate (Punica granatum L.)

Here we report on comprehensive chloroplast (cp) genome analysis of 16 pomegranate (Punica granatum L.) genotypes representing commercial cultivars, ornamental and wild types, through large-scale sequencing and assembling using next-generation sequencing (NGS) technology. Comparative genome analysis revealed that the size of cp genomes varied from 158,593 bp (in wild, “1201” and “1181”) to 158,662 bp (cultivar, “Gul-e-Shah Red”) among the genotypes, with characteristic quadripartite structures separated by a pair of inverted repeats (IRs). The higher conservation for the total number of coding and non-coding genes (rRNA and tRNA) and their sizes, and IRs (IR-A and IR-B) were observed across all the cp genomes. Interestingly, high variations were observed in sizes of large single copy (LSC, 88,976 to 89,044 bp) and small single copy (SSC, 18,682 to 18,684 bp) regions. Although, the structural organization of newly assembled cp genomes were comparable to that of previously reported cp genomes of pomegranate (“Helow,” “Tunisia,” and “Bhagawa”), the striking differences were observed with the Lagerstroemia lines, viz., Lagerstroemia intermedia (NC_0346620) and Lagerstroemia speciosa (NC_031414), which clearly confirmed previous findings. Furthermore, phylogenetic analysis also revealed that members outside the genus Punica were clubbed into a separate clade. The contraction and expansion analysis revealed that the structural variations in IRs, LSC, and SSC have significantly accounted for the evolution of cp genomes of Punica and L. intermedia over the periods. Microsatellite survey across cp genomes resulted in the identification of a total of 233 to 234 SSRs, with majority of them being mono- (A/T or C/G, 164–165 numbers), followed by di- (AT/AT or AG/CT, 54), tri- (6), tetra- (8), and pentanucleotides (1). Furthermore, the comparative structural variant analyses across cp genomes resulted in the identification of many varietal specific SNP/indel markers. In summary, our study has offered a successful development of large-scale cp genomics resources to leverage future genetic, taxonomical, and phylogenetic studies in pomegranate.

Complete chloroplast genome sequence of Oxytropis glabra (Leguminosae)

Oxytropis glabra DC. is a perennial poisonous plant to livestock belonging to the genus Oxytropis, Leguminosae, mainly distributed in Northwestern China. As a poisonous grass, this species protects plant diversity in degraded grasslands by sheltering adjacent plants. In this study, the complete chloroplast genome with a total size of 122,094 bp was reported. Our annotations showed that the chloroplast genome contains 109 genes, including 76 protein-coding genes, 29 tRNA genes, and four rRNA genes. This work presents complete chloroplast genome information, which will be valuable for studying the evolution and genetic diversity of O. glabra.

Complete plastid genome sequence of Halimodendron halodendron (Leguminosae)

Halimodendron halodendron (Pall.) Voss. is a deciduous shrub belonging to the genus Halimodendron, Leguminosae, and is mainly distributed in dry areas. This species can be used for saline-alkali soil improvement and sand fixation. The complete plastid genome of H. halodendron first reported here is 129,342 bp in length, and contains 110 genes, including 76 protein coding genes, 30 tRNA genes, and 4 rRNA genes. A total of 105 simple sequence repeats (SSRs) were identified in the chloroplast genome. This information will be useful for study on the evolution and genetic diversity of Halimodendron halodendron in the future.

The complete chloroplast genome of the subtropical species Camellia japonica 'Huaheling'

Camellia japonica ‘Huaheling’ is a rare subtropical Camellia species in China with high ornamental and medicinal value. The complete chloroplast genome of C. japonica ‘Huaheling’ is a 157,001-bp circular DNA molecule containing a large single-copy region (LSC, 86,704 bp), a small single-copy region (SSC, 18,393 bp), and two inverted repeat sequences (IR). Of the 131 genes identified, 86 are protein-coding genes, 8 are rRNA genes, and 37 are tRNA genes. A total of 54 simple sequence repeats (SSRs) were identified in the chloroplast genome. The phylogenetic analysis showed that C. japonica ‘Huaheling’ is clustered with C. japonica. This work provides valuable information for future study of the evolution and genetic diversity of C. japonica ‘Huaheling.’

Characterization of the complete chloroplast genome of Alopecurus pratensis L. (Poaceae)

Alopecurus pratensis L. is one of the most important fodder grasses distributed in sub-frigid regions of the world. In this study, the complete chloroplast genome of A. pratensis was deciphered and is 136,157 bp in length. The genome includes a large single-copy region of 80,275 bp, small single-copy region of 12,830 bp, and a pair of inverted repeat regions of 21,526 bp. The GC content of the complete chloroplast genome of A. pratensis is 38.30%. Among the 134 unique genes in the circular genome, 38 tRNA, eight rRNA, and 88 protein-coding genes were annotated. We constructed the Maximum likelihood (ML) tree with 13 species from the Poaceae and found that A. pratensis was phylogenetically related to A. arundinaceus. The published A. pratensis chloroplast genome will provide useful information for phylogenetic and evolutionary study of the genus Alopecurus in the Poaceae.

Chromosomal-level assembly of the Leptodermis oblonga (Rubiaceae) genome and its phylogenetic implications

Rubiaceae is the fourth largest and a taxonomically complex family of angiosperms. Many species in this family harbor low reproductive isolation and frequently exhibit inconsistent phenotypic characteristics. Therefore, taxonomic classification and their phylogenetic relationships in the Rubiaceae family is challenging, especially in the genus Leptodermis. Considering the low taxonomic confusion and wide distribution, Leptodermis oblonga is selected as a representative Leptodermis for genome sequencing. The assemblies resulted in 497 Mbp nuclear and 155,100 bp chloroplast genomes, respectively. Using the nuclear genome as a reference, SNPs were called from 37 Leptodermis species or varieties. The phylogenetic tree based on SNPs exhibited high resolution for species delimitation of the complex and well-resolved phylogenetic relationships in the genus. Moreover, 28,987 genes were predicted in the nuclear genome and used for comparative genomics study. As the first chromosomal-level genome of the subfamily Rubioideae in Rubiaceae, it will provide fruitfully evolutionary understanding in the family.

Complete chloroplast genome of medicinal plant Sabia parviflora Wall. ex Roxb. (Sabiaceae)

Sabia parviflora Wall. ex Roxb., an evergreen climbing woody vine, is a Chinese herbal medicine commonly used by ethnic minorities in some areas of China. In this study, the chloroplast genome of S. parviflora was sequenced for the first time. Its genome is 162,054 bp in length with 38.6% of GC content. The genome consists of a large single copy (LSC) region of 90,001 bp, a small single copy (SSC) region of 18,887 bp, and two inverted repeat (IRa and IRb) regions of 26,583 bp each. A total of 130 genes were annotated, including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis was conducted by nine species from order Proteales, which demonstrated a close relationship between the family Sabiaceae and Nelumbonaceae.

The complete chloroplast genome of Atractylodes japonica Koidz. ex Kitam. and its phylogenetic inference

Atractylodes japonica Koidz. ex Kitam. is a perennial herbal plant, and its dried rhizomes have been widely used as traditional medicine in China and Japan. In this study, we assembled and annotated the complete chloroplast (cp) genome sequence of A. japonica using the high-throughput sequencing approach. The cp genome of A. japonica is 153,208 bp in length with the overall GC content of 37.7%, including two inverted repeat (IR) regions of 25,147 bp, which was separated by a large single-copy (LSC) region of 84,255 bp and a small single-copy (SSC) region of 18,659 bp. 113 unique genes were annotated in the genome, including 80 protein-coding genes, 29 represented tRNA genes, and four denoted rRNA genes. A maximum-likelihood phylogenetic analysis with 38 complete cp sequences showed that Atractylodes formed a monophyletic clade, and A. japonica and A. koreana formed a subclade in Atractylodes. This study provides the chloroplast genome structure features and phylogenetic relationship of A. japonica.

The complete chloroplast genome of Bromus biebersteinii

Bromus biebersteinii is a perennial gramineous grass, which is mainly distributed in Southwest Asia. In our study, we obtained the complete chloroplast genome of B. biebersteinii and found it is 137,189 bp in length. The GC content of its whole chloroplast genome is 38.37%. Among the 134 unique genes in the circular genome, 38 tRNA, 8 rRNA and 88 protein-coding genes were successfully annotated. We constructed the Maximum likelihood (ML) tree with 12 species, and found that B. biebersteinii was phylogenetically close to Bromus vulgaris of the genus Bromus.

The complete chloroplast genome of Androsace erecta (Primulaceae) and its phylogenetic implication

With about 153 species, the genus Androsace (Primulaceae) is known for its horticultural and economic importance. In this study, we report the complete chloroplast genome of Androsace erecta Maximowicz, a morphologically distinct species of Sect. Orthocaulon native to the Western China. The plastome of A. erecta is highly conserved in genome size, structure, and content when compared to all previously published plastomes of the genus. The phylogenomic analysis strongly supported A. erecta as sister to a clade comprising species of Sections Aizoideia and Chamaejasme. Lastly, we selected the four most variable regions across the Androsace species plastomes (trnK^UUU-rps16, trnS^GCU-trnG^UCC, psbE-petL, and infA-rps8), which were considered to be suitable candidate DNA barcodes for Androsace.

Chromosome-level reference genome of the soursop (Annona muricata): A new resource for Magnoliid research and tropical pomology

The flowering plant family Annonaceae includes important commercially grown tropical crops, but development of promising species is hindered by a lack of genomic resources to build breeding programs. Annonaceae are part of the magnoliids, an ancient lineage of angiosperms for which evolutionary relationships with other major clades remain unclear. To provide resources to breeders and evolutionary researchers, we report a chromosome-level genome assembly of the soursop (Annona muricata). We assembled the genome using 444.32 Gb of DNA sequences (676× sequencing depth) from PacBio and Illumina short-reads, in combination with 10× Genomics and Bionano data (v1). A total of 949 scaffolds were assembled to a final size of 656.77 Mb, with a scaffold N50 of 3.43 Mb (v1), and then further improved to seven pseudo-chromosomes using Hi-C sequencing data (v2; scaffold N50: 93.2 Mb, total size in chromosomes: 639.6 Mb). Heterozygosity was very low (0.06%), while repeat sequences accounted for 54.87% of the genome, and 23,375 protein-coding genes with an average of 4.79 exons per gene were annotated using de novo, RNA-seq and homology-based approaches. Reconstruction of the historical population size showed a slow continuous contraction, probably related to Cenozoic climate changes. The soursop is the first genome assembled in Annonaceae, supporting further studies of floral evolution in magnoliids, providing an essential resource for delineating relationships of ancient angiosperm lineages. Both genome-assisted improvement and conservation efforts will be strengthened by the availability of the soursop genome. As a community resource, this assembly will further strengthen the role of Annonaceae as model species for research on the ecology, evolution and domestication potential of tropical species in pomology and agroforestry.

The complete chloroplast genome of Alhagi sparsifolia Shap. (Leguminosae)

Alhagi sparsifolia Shap. is a perennial herbaceous plant belonging to the genus Alhagi, Leguminosae. This species is of high nutritional, medicinal and ecological values. The complete chloroplast genome was 128,418 bp and lost an IR (inverted repeat) region. Further annotation revealed the chloroplast genome contains 108 genes, including 75 protein coding genes, 29 tRNA genes, and 4 rRNA genes. A total of 103 simple sequence repeats (SSRs) were identified in the chloroplast genome. This chloroplast genome resource will be useful for study on the evolution and genetic diversity of A. sparsifolia in the future.

New Insight into Taxonomy of European Mountain Pines, Pinus mugo Complex, Based on Complete Chloroplast Genomes Sequencing

The Pinus mugo complex is a large group of closely related mountain pines, which are an important component of the ecosystems of the most important mountain ranges, such as the Alps, Carpathians and Pyrenees. The phylogenetic relationships between taxa in this complex have been under discussion for many years. Despite the use of many different approaches, they still need to be clarified and supplemented with new data, especially those obtained with high-throughput methods. Therefore, in this study, the complete sequences of the chloroplast genomes of the three most recognized members of the Pinus mugo complex, i.e., Pinus mugo, Pinus rotundata and Pinus uncinata, were sequenced and analyzed to gain new insight into their phylogenetic relationships. Comparative analysis of their complete chloroplast genome sequences revealed several mutational hotspots potentially useful for the genetic identification of taxa from the Pinus mugo complex. Phylogenetic inference based on sixteen complete chloroplast genomes of different coniferous representatives showed that pines from the Pinus mugo complex form one distinct monophyletic group. The results obtained in this study provide new and valuable omics data for further research within the European mountain pine complex. They also indicate which regions may be useful in the search for diagnostic DNA markers for the members of Pinus mugo complex and set the baseline in the conservation of genetic resources of its endangered taxa.

The complete chloroplast genome of Euphorbia hirta (Euphorbiaceae), a commonly used medicinal plant in China

Plants in the genus Euphorbia have been widely used as herbal medicine, and for ornamental horticulture and biofuel production. In this study, we report the complete chloroplast genome of Euphorbia hirta which is known as the ‘asthma-plant’ due to its medicinal use. The chloroplast genome of this species is 164,340 bp in length, including a pair of inverted repeat regions (IRs) (27,354 bp) that are divided by a large single-copy region (LSC) (91,373 bp) and a small single-copy region (SSC) (18,259 bp). The chloroplast genome of E. hirta contains 111 unique genes (77 protein-coding, 30 tRNA, and four rRNA), 19 of which are duplicated in the IR regions. The overall GC content is 35.4%. Phylogenetic analysis fully resolved E. hirta groups with other species of Euphorbia. The complete chloroplast genome of E. hirta provides useful information that can be used to distinguish related species and reconstruct evolutionary relationships.

Characterization of the complete chloroplast genome of mangrove Rhizophora apiculata Blume (Rhizophoraceae)

The chloroplast (cp) genome sequence of Rhizophora apiculata was characterized. The cp genome length was 164,343 bp in length, containing a typical structure of a large single copy (LSC) of 93,155 bp, a small single copy (SSC) of 19,376 bp, and two inverted repeats (IRs) of 25,906 bp, with a GC content of 34.9%. There were 131 genes were annotated in the cp genome, including 85 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. A phylogenetic analysis using cp genomes of mangroves and ecologically associated species resolved R. apiculata in Rhizophora with R. stylosa and R. x lamarckii. This complete chloroplast sequence offers a promising tool for further species identification and evolutionary studies of Rhizophora, as well as for mangroves.

Characterization of the complete chloroplast genome of mangrove Bruguiera gymnorrhiza (L.) Lam. ex Savigny

The chloroplast (cp) genome sequence of Bruguiera gymnorrhiza was characterized. The cp genome length was 163,795 bp in length, with a GC content of 35.3%, containing a large single copy (LSC) of 90,830 bp, a small single copy (SSC) of 20,207 bp, and a pair of inverted repeats (IRs) of 26,379 bp. The genome contained 121 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. A phylogenetic analysis using cp genomes of mangroves and ecologically associated species resolved B. gymnorrhiza in Bruguiera with B. sexangula var. rhynchopetala. This complete chloroplast sequence offers a promising tool for further species identification and evolutionary studies of Bruguiera, as well as for mangroves.

The complete chloroplast genome of Atractylodes koreana (Nakai) Kitam and its phylogenetic analysis

Atractylodes koreana (Nakai) Kitam is a perennial herb of Asteraceae, mainly distributed in China and Korea, which is the main adulterant of traditional herbal medicine ‘Cangzhu’. In the present study, we reported the complete chloroplast (cp) genome of A. koreana with the total length of 153,232 bp, which is consisted of four regions, including one large single copy (LSC) region of 84,250 bp, one small single copy (SSC) region of 18,690 bp, and two inverted repeat regions (IRa and IRb) of 25,146 bp. The GC content of the complete cp genome is 37.7%. A total of 110 unique genes were annotated, comprising 79 protein-coding genes, 27 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Moreover, nine protein-coding genes contained one intron and three protein-coding genes (clpP, ycf3, and rps12) contained two introns. The phylogenetic analysis indicated that A. koreana is a sister group of A. chinensis and A. lancea.

Comparative genome analysis revealed gene inversions, boundary expansions and contractions, and gene loss in the Stemona sessilifolia (Miq.) Miq. chloroplast genome

Stemona sessilifolia (Miq.) Miq., commonly known as Baibu, is one of the most popular herbal medicines in Asia. In the Chinese Pharmacopoeia, Baibu has multiple authentic sources and there are many similar herbs sold as Baibu in herbal medicine markets. The existence of counterfeits of Baibu brings challenges to its identification. To assist in its accurate identification, we sequenced and analyzed the complete chloroplast genome of S. sessilifolia using next-generation sequencing technology. The genome was found to be 154,037 bp in length, possessing a typical quadripartite structure consisting of a pair of inverted repeats (IRs: 27,090 bp) separated by a large single copy (LSC: 81,949 bp) and a small single copy (SSC: 17,908 bp). A total of 112 unique genes were identified, including 80 protein-coding, 28 transfer RNA and four ribosomal RNA genes. In addition, 45 tandem, 27 forward, 23 palindromic and 104 simple sequence repeats were detected in the genome by repeated analysis. Compared with its counterfeits (Asparagus officinalis and Carludovica palmata) we found that IR expansion and SSC contraction events of S. sessilifolia resulted in two copies of the rpl22 gene in the IR regions and a partial duplication of the ndhF gene in the SSC region. An approximately 3-kb-long inversion was also identified in the LSC region, leading to the petA and cemA genes being presented in the complementary strand of the chloroplast DNA molecule. Comparative analysis revealed some highly variable regions, including trnF-GAA_ndhJ, atpB_rbcL, rps15_ycf1, trnG-UCC_trnR-UCU, ndhF_rpl32, accD_psaI, rps2_rpoC2, trnS-GCU_trnG-UCC, trnT-UGU_trnL-UAA and rps16_trnQ-UUG. Finally, gene loss events were investigated in the context of phylogenetic relationships. In summary, the complete plastome of S. sessilifolia will provide valuable information for the distinction between Baibu and its counterfeits and assist in elucidating the evolution of S. sessilifolia.

The complete chloroplast genome of Agropyron desertorum (Fisch. ex Link) Schult

Agropyron desertorum is one of the most important fodder grasses which distributes in the temperate regions of the world. In this study, the complete chloroplast genome of A. desertorum was sequenced. The genome was 135,459 bp in length, with a large single-copy region of 79,617 bp, a small single-copy region of 12,760 bp, and two inverted repeat regions of 21,541 bp. The genome consisted of 132 genes, including 86 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. The GC contents was 38.32%. We constructed the Maximum likelihood (ML) tree with 13 species from the Hordeinae genus, and concluded that A. desertorum was closely related to plant species of the Elymus genus in the phylogenetic tree.

Characterization of the complete chloroplast genome of Clerodendrum bungei Steud. (Lamiaceae)

Clerodendrum bungei Steud. is one kind of traditional medical herb which can be used for airway hyperreactivity treatment. In this study, the complete chloroplast genome sequence of C. bungei was assembled. Its complete circular chloroplast DNA length was 151,680 bp. The genome was made up of a large single-copy region of 83,189 bp, a small single-copy region of 17,311 bp, and a pair of inverted repeat regions of 25,590 bp. The genome totally encoded 130 genes, containing 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The phylogenetic analysis indicates that C. bungei belongs to the Lamiaceae family.

The complete chloroplast genome of Verbena officinalis, an herbal species of Verbenaceae family

Verbena officinalis is one kind of traditional medical herb which has potential for multiple diseases’ treatment. In this study, the complete chloroplast genome sequence of V. officinalis was assembled. Its complete circular chloroplast DNA length was 153,491 bp. The genome was made up of a large single-copy region of 84,518 bp, a small single-copy region of 17,357 bp, and a pair of inverted repeat regions of 25,808 bp. The genome totally encoded 128 genes, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis indicates that V. officinalis belongs to the verbenaceae family.

Genome skimming reveals novel plastid markers for the molecular identification of illegally logged African timber species

Tropical forests represent vast carbon stocks and continue to be key carbon sinks and buffer climate changes. The international policy constructed several mechanisms aiming at conservation and sustainable use of these forests. Illegal logging is an important threat of forests, especially in the tropics. Several laws and regulations have been set up to combat illegal timber trade. Despite significant enforcement efforts of these regulations, illegal logging continues to be a serious problem and impacts for the functioning of the forest ecosystem and global biodiversity in the tropics. Microscopic analysis of wood samples and the use of conventional plant DNA barcodes often do not allow to distinguish closely-related species. The use of novel molecular technologies could make an important contribution for the identification of tree species. In this study, we used high-throughput sequencing technologies and bioinformatics tools to obtain the complete de-novo chloroplast genome of 62 commercial African timber species using the genome skimming method. Then, we performed a comparative genomic analysis that revealed new candidate genetic regions for the discrimination of closely-related species. We concluded that genome skimming is a promising method for the development of plant genetic markers to combat illegal logging activities supporting CITES, FLEGT and the EU Timber Regulation.

The complete chloroplast genome of Cymbidium dayanum (Orchidaceae)

Cymbidium dayanum, a wild orchid species in the Orchid family (Orchidaceae), is considered highly valuable because of its long flowering period and beautiful plant shape. We sequenced the complete chloroplast genome of C. dayanum using the Illumina Hiseq platform (Illumina, San Diego, CA). The size of the C. dayanum chloroplast genome is 155,408 bp, with an average GC content of 36.76%. This chloroplast genome has containing a large single-copy (LSC) region of 84,189 bp, a small single-copy (SSC) region of 17,991 bp, and two inverted (IRa and IRb) repeat regions of two 26,614 bp. A total of 118 unique genes were annotated, including 76 protein-coding genes, 38 tRNA genes, and 4 rRNA genes. A maximum-likelihood phylogenetic tree indicated that C. dayanum is closely related to C. tracyanum in the genus Cymbidium based on 9 whole chloroplast genome sequences.

The complete chloroplast genome of Pittosporum brevicalyx

Here, we report the complete chloroplast genome of Pittosporum brevicalyx. The genome is 153,388 bp in size, which is comprised of a large single-copy (LSC) region of 84,724 bp, a small single-copy (SSC) region of 18,716 bp, and two inverted repeat (IR) regions of 24,974 bp. The overall GC content of the plastome was 38.3%. The new sequence comprised 127 unique genes, including 82 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. Phylogenetic analysis showed that P.brevicalyx was close to Pittosporum kerrii and Pittosporum eugenioides. These data may providing useful information for phyletic evolution of P.brevicalyx within the Pittosporaceae family.

Characterization of the complete chloroplast genome of Gymnocladus chinensis Baill

Gymnocladus chinensis Baill (Fam.: Leguminosae; Trib.: Caesalpinieae) are widely distributed in China. In this study, we assembled the complete chloroplast (cp) genome of G. chinensis. The total cp genome size was 165,315 bp in length, containing a large single-copy region of 92,356 bp, a small single-copy region of 20,449bp, and a pair of inverted repeat regions of 26,255 bp. The all GC content of G. chinensis cp was 34.95%. It encodes a total of 105 unique genes, including 75 protein-coding genes, 26 tRNA genes, and four rRNA genes. Seventeen genes contain a single intron, and two genes (ycf3 and clpP) have two introns. Phylogenetic analysis results strongly supported that G. chinensis was closely related to Angylocalyx braunii.

NOVOWrap: An automated solution for plastid genome assembly and structure standardization

Plastid genomes play an important role in genomics and evolutionary biology. Next-generation sequencing has revolutionized plastid genomic data acquisition to the point that genome assembly has become a bottleneck for widespread utilization of plastid genome data. To solve this problem, we developed an open-source, cross-platform tool known as, NOVOWrap, which includes both command-line and graphical interfaces for automatically assembling plastid genomes on personal computers. With minimal inputs, settings, and user intervention, NOVOWrap can automatically assemble plastid genomes, validate results and standardize the structure using affordable computer resources. The performance of this software has been successfully benchmarked against the plastid genomes of 11 species belonging to lycopods, gymnosperms, and angiosperms. By liberating researchers from laborious and cumbersome computer manipulations and create reliable and standardized genomic data, NOVOWrap is expected to accelerate plastid genome assembly, ease the process of data exchange, and contribute to downstream analysis.

The complete chloroplast genome of Rhododendron datiandingense (Ericaceae)

Rhododendron datiandingense is newly reported and endemic to China. The genome of R. datiandingense is 207,311 bp in length, including a large single-copy region of 190,689 bp and a small single-copy region of 2582 bp, a pair of inverted repeat regions (IRA) of 7020 bp each. The genome encodes 110 genes, comprising 77 protein-coding genes, four ribosomal RNA genes, and 29 transfer RNA genes. Repeat analysis revealed 62 simple sequence repeats (SSRs) in the genome. Phylogenetic analysis revealed that R. datiandingense is clearly separated from the other Rhododendron species and shown in the basal position.

The complete plastome of Andreaea rupestris Hedw. (Andreaeaceae, Bryophyta)

Andreaea rupestris Hedw., one of the lantern mosses, is the lectotype of the genus Andreaea Hedw. (Andreaeaceae). Here we present its complete plastome. The plastome of A. rupestris is successfully assembled from raw reads sequenced by HiSeq X ten system. Its total length is 135,214 bp consisting of four regions: large single copy (LSC) region (92,780 bp), small single copy (SSC) region (21,102 bp), and two inverted repeat regions (IRs; 10,666 bp per each). It contains 134 genes (88 coding genes, 8 rRNAs, and 38 tRNAs). The overall GC content is 30.3% and in the LSC, SSC, and IR regions are 27.5%, 26.5%, and 46.2%, respectively. The present data will be an important sequence resource for further studies on the important early diverging lineage of mosses.

The complete plastome of Polytrichum commune Hedw. (Polytrichaceae, Bryophyta)

Polytrichum commune, one of hair-cap mosses, is the type species of the genus Polytrichum Hedw. (Polytrichaceae). Here we present its complete plastome. The plastome of P. commune is successfully assembled from raw reads sequenced by HiSeq X ten system. Its total length is 126,323 bp consisting of four regions: large single copy (LSC) region (88,070 bp), small single copy (SSC) region (16,717 bp), and inverted repeats (IRs; 9,680 bp per each). It contains 128 genes (84 coding genes, eight rRNAs, and 36 tRNAs); nine genes (four rRNAs and five tRNAs) are duplicated in IR regions. The overall GC content is 28.9% and in the LSC, SSC and IR regions is 26.1%, 25.1%, and 45.5%, respectively. This plastome is an important sequence resource for further studies on the class Polytrichopsida.

Complete chloroplast genomes of Achnatherum inebrians and comparative analyses with related species from Poaceae

This article reports the complete chloroplast genome of Achnatherum inebrians, a poisonous herb that is widely distributed in the rangelands of Northern China. The genome is 137 714 bp in total and consists of a large single‐copy (81 758 bp) region and small single‐copy (12 682 bp) region separated by a pair of inverted repeats (21 637 bp). The genome contains 130 genes, including 84 protein‐coding genes, 38 tRNA genes and 8 ribosomal RNA genes, and the guanine + cytosine content is 36.17%. We subsequently performed comparative analysis of complete genomes from A. inebrians and other Poaceae‐related species from GenBank. Thirty‐eight simple sequence repeats were identified, further demonstrating rapid evolution in Poaceae. Finally, the phylogenetic trees of 37 species of Poaceae and 2 species of Amaranthaceae were constructed by using maximum likelihood and Bayesian inference methods, based on the genes of the complete chloroplast genome. We identified hotspots that can be used as molecular markers and barcodes for phylogenetic analysis, as well as for species identification. Phylogenetic analysis indicated that A. inebrians is a member of the genus Stipa rather than Achnatherum.

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.

Decoding first complete chloroplast genome of toothbrush tree (Salvadora persica L.): insight into genome evolution, sequence divergence and phylogenetic relationship within Brassicales

BACKGROUND

Salvadora persica L. (Toothbrush tree - Miswak; family-Salvadoraceae) grows in the arid-land ecosystem and possesses economic and medicinal importance. The species, genus and the family have no genomic datasets available specifically on chloroplast (cp) genomics and taxonomic evolution. Herein, we have sequenced the complete chloroplast genome of S. persica for the first time and compared it with 11 related specie's cp genomes from the order Brassicales.

RESULTS

The S. persica cp genome was 153,379 bp in length containing a sizeable single-copy region (LSC) of 83,818 bp which separated from the small single-copy region (SSC) of 17,683 bp by two inverted repeats (IRs) each 25,939 bp. Among these genomes, the largest cp genome size (160,600 bp) was found in M. oleifera, while in S. persica it was the smallest (153,379 bp). The cp genome of S. persica encoded 131 genes, including 37 tRNA genes, eight rRNA genes and 86 protein-coding genes. Besides, S. persica contains 27 forward, 36 tandem and 19 palindromic repeats. The S. persica cp genome had 154 SSRs with the highest number in the LSC region. Complete cp genome comparisons showed an overall high degree of sequence resemblance between S. persica and related cp genomes. Some divergence was observed in the intergenic spaces of other species. Phylogenomic analyses of 60 shared genes indicated that S. persica formed a single clade with A. tetracantha with high bootstrap values. The family Salvadoraceae is closely related to Capparaceae and Petadiplandraceae rather than to Bataceae and Koberliniacaea.

CONCLUSION

The current genomic datasets provide pivotal genetic resources to determine the phylogenetic relationships, genome evolution and future genetic diversity-related studies of S. persica in complex angiosperm families.

Precise Species Detection in Traditional Herbal Patent Medicine, Qingguo Wan, Using Shotgun Metabarcoding

As one of the high-incidence diseases in the world, pharyngitis seriously affects the lives of those with the condition. Qingguo Wan is a herbal medicine used for treating pharyngitis, and its quality evaluation is currently only accomplished via traditional identification. However, precise identification becomes challenging with fake products on the market or fungal contamination during the production process. This study used the Illumina NovaSeq platform for targeting the ITS2, psbA-trnH, matK, and rbcL sequences to survey the species composition of lab-made and commercial samples. The results showed that a total of 34.56 Gb of raw data that was obtained represented more than 0.23 billion reads. After assembly, annotation, and operational taxonomic unit clustering, 103, 12, 10, and 12 OTUs were obtained, which belonged to the ITS2, psbA-trnH, matK, and rbcL sequences of the mock lab-made and commercial samples. The analytical results indicated that the sequences of all the prescription ingredients were successfully obtained in the two lab-made samples. The positive control medicinal Panax quinquefolius L. sequence was obtained in HSZY175, while Scutellaria baicalensis Georgi, Lonicera japonica Thunb. Menispermum dauricum DC. and Paeonia lactiflora Pall. were detected in the three commercial samples. The detection results of the other four herbs in different fragments were not all the same. In addition, a total of 28 fungi OTUs, representing 19 families and 20 genera, were obtained from both the commercial and mock lab-made samples. Aspergillus, Cladosporium, and Penicillium dominated among the 20 genera. This study demonstrated that the shotgun metabarcoding method is a powerful tool for the molecular identification of the biological ingredients in Qingguo Wan. It can be used to effectively supplement traditional methods while providing a new technique for the quality evaluation of Qingguo Wan.

Comparative analysis of the complete chloroplast genomes of Cirsium japonicum from China and Korea

Cirsium japonicum (C. japonicum) is a traditional Chinese medicine belonging to the family Asteraceae. The previous studies have indicated that the chemical compound content of C. japonicum from different places was different. To distinguish C. japonicum from different geographies, the chloroplast genome of C. japonicum from China was sequenced and compared with that from Korea. The total length of this genome is 152,602 bp, similar to that of Korea (152,606 bp). It has a conservative quartile structure which is composed of a large single-copy (LSC) region, a small single-copy (SSC) region and a pair of inverted repeats (IRs) regions, with lengths of 83,487 bp, 18,721 bp, and 25,197 bp, respectively. It encodes 79 protein-coding, 27 transfer RNAs, and 4 ribosomal RNA genes. The overall GC content of the genome is 37.70%. A total of 20 single nucleotide polymorphisms and 6 insertions and deletions were identified between the chloroplast genome of C. japonicum from China and Korea. These results can be applied to develop molecular markers to distinguish C. japonicum from different geographical origins.

Comparative and phylogenetic analysis of the complete chloroplast genome sequences of Lactuca raddeana and Lactuca sativa

Lactuca raddeana is a biennial plant of the Lactuca genus belonging to the Asteraceae family. The classification of Lactuca is controversial, and no consistent conclusions have been reached based on the analysis of morphological characters and different molecular markers. Here, we sequenced and assembled the complete chloroplast genome of L. raddeana. This genome has a total length of 152,339 bp, a conservative quartile structure that is composed of a large single-copy (LSC) region of 83,976 bp, a small-copy (SSC) region of 18,521 bp, and a pair of inverted repeats (IRs) of 24,921 bp. The genome contains 112 unique genes, including 79 protein-coding, four rRNA, and 29 tRNA genes. Repeat analysis obtained 17 microsatellite, 16 tandem, and 17 interspersed repeats. Comparison of sequence divergence between L. raddeana and L. sativa found the intergenic spacer trnC-petN exhibited the highest degree of variation. Three phylogenetic trees based on the 72 shared protein, matK gene, and rbcL gene sequences showed that L. raddeana and L. sativa are clustered together. The acquisition and comparative analysis of the chloroplast genome provide a valuable resource for the taxonomic and phylogenetic studies of Lactuca.

The complete chloroplast genome of Photinia davidsoniae: molecular structures and comparative analysis

Photinia davidsoniae is a common ornamental arbor in the genus Photinia (family Rosaceae). Here, we sequenced and assembled the complete plastome of P. davidsoniae using the next-generation DNA sequencing technology. And we then compared it with nine Photinia species using a range of bioinformatics software tools. The ten plastomes had sizes ranging from 159,230 bp for P. beckii to 160,346 bp for P. davidsoniae. They all had a conservative quartile structure. It contained two single-copy regions: a large single-copy (LSC) region, a small single-copy (SSC) region, and a pair of inverted repeat (IR) regions. Each of the plastomes encoded 113 unique genes, including 79 protein-coding genes, four rRNA genes, and 30 tRNA genes. Furthermore, we detected six hypervariable regions (matK-rps16, rpoB-trnC, trnT-psbD, ndhC-trnV, psbE-petL, ndhF-rpl32-trnL), which could be used as potential molecular markers. We constructed two phylogenetic trees with plastomes or concatenated protein sequences of 25 species of 8 genera of Rosaceae. The tree constructed with complete plastomes has much stronger support. The results placed P. davidsoniae in the upper part of the phylogenetic tree. It shows that P. davidsoniae and P. lanuginosa are closely related. In summary, the plastomes of Photinia are conserved overall but carry significant minor variations, as expected. The results will be indispensable for distinguishing species, understanding the interspecific diversity, and elucidating the evolutionary processes of Photinia species.

Elatostema qinzhouense (Urticaceae), a new species from limestone karst in Guangxi, China

Elatostema qinzhouense L.F. Fu, A.K. Monro & Y.G. Wei, a new species from Guangxi, China is described and illustrated. Morphologically, E. qinzhouense is most similar to E. hezhouense from which it differs by having smaller size of leaf laminae, fewer and smaller staminate peduncle bracts, longer pistillate peduncle bracts and a larger achene. This result is supported by the molecular evidence. The phylogenetic position of the new species within Elatostema is evaluated using three DNA regions, ITS, trnH-psbA and psbM-trnD, for 107 taxa of Elatostema s.l. (including E. qinzhouense). Bayesian inference (BI) and maximum likelihood (ML) analyses each recovered the same strongly supported tree topologies, indicating that E. qinzhouense is a member of the core Elatostema clade and sister to E. hezhouense. Along with the phylogenetic studies, plastid genome and ribosomal DNA (rDNA) sequences of the new species are assembled and annotated. The plastid genome is 150,398 bp in length and comprises two inverted repeats (IRs) of 24,688 bp separated by a large single-copy of 83,919 bp and a small single-copy of 17,103 bp. A total of 113 functional genes are recovered, comprising 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The rDNA is 5,804 bp in length and comprised the 18S ribosomal RNA partial sequence (1,809 bp), internal transcribed spacer 1 (213 bp), 5.8S ribosomal RNA (164 bp), internal transcribed spacer 2 (248 bp) and 26S ribosomal RNA partial sequence (3,370 bp). In addition, the chromosome number of E. qinzhouense is observed to be 2n = 26, suggesting that the species is diploid. Given a consistent relationship between ploidy level and reproductive system in Elatostema, the new species is also considered to be sexually reproducing. Our assessment of the extinction threat for E. qinzhouense is that it is Endangered (EN) according to the criteria of the International Union for Conservation of Nature.

Complete chloroplast genome sequence of Peganum harmala (Zygophyllaceae)

Peganum harmala L. is a perennial herbaceous plant belonging to the family of Zygophyllaceae, and is grows in semi-arid climates, such as Xinjiang, Gansu, Ningxia, Qinghai, and Inner Mongolia in China, and also Middle East and North Africa. This species is of high medicinal value. The complete chloroplast genome was reported in this study. The chloroplast genome with a total size of 159,957 bp consists of two inverted repeats (IR, 26,550 bp) separated by a large single-copy region (LSC, 88,098 bp) and a small single-copy region (SSC, 18,759 bp). Further annotation revealed the chloroplast genome contains 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. A total of 90 simple sequence repeats (SSRs) were identified in the chloroplast genome. This information will be useful for study on the evolution and genetic diversity of Peganum harmala in the future.

The complete chloroplast genome of Bupleurum euphorbioides, a traditional medicinal plant

Bupleurum euphorbioides is a rare native plant attributed with analgesic, gallbladder-supportive, and other functions in China and the Republic of Korea. However, the complete chloroplast genome sequence of the native plant B. euphorbioides has not been determined. In this study, we sequenced the complete chloroplast genome sequence, and examined the molecular phylogeny and genetic information of B. euphorbioides. The total chloroplast genome of B. euphorbioides was 154,871 bp in length with a large single-copy region (85,089 bp), small single-copy region (17,714 bp), and pair of inverted repeats regions (26,034 bp). The chloroplast genome encoded a total of 176 genes, including 131 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The phylogenetic tree indicated that B. euphorbioides was the most closely related to B. latissimum.

The complete chloroplast genome of Melicope pteleifolia (Champ. ex Benth.) T. G. Hartley (Rutaceae)

Melicope pteleifolia, an important medicinal and horticultural plant, has antipyretic, anti-inflammatory, and analgesic effects. Here, the complete chloroplast genome of M. pteleifolia was sequenced and its phylogenetic relationship was investigated. The complete chloroplast genome of M. pteleifolia was 159,014 bp in size, including a pair of inverted repeat regions (IR, 27,640 bp), a large single copy region (LSC, 85,124 bp), and a small single copy region (SSC, 18,610 bp). The GC content of the chloroplast genome was 38.3%. A total of 133 genes were annotated, including 88 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. Phylogenetic analysis revealed that M. pteleifolia clusters together with species of Toddalia, Zanthoxylum, Tetradium, Phellodendron, and Casimiroa.

Peculiarities of the inverted repeats in the complete chloroplast genome of Strobilanthes bantonensis Lindau

Strobilanthes bantonensis Lindau belongs to the family Acanthaceae. It is an antiviral herb that can be used to prevent Influenza virus infections in the border areas between China and Vietnam. Local people call it ‘Purple Ban-lan-gen’ because its root is very similar to that of Strobilanthes cusia (Nees) Kuntze, which is called ‘Southern Ban-lan-gen’ and is listed in Chinese Pharmacopeia. The two species have been used interchangeably locally. However, their pharmacological equivalence has caused concern for years. We have sequenced the chloroplast genome of S. cusia previously. In this study, we sequenced the complete chloroplast genome sequence of S. bantonensis to preform in-depth comparative genetic analysis of the two Strobilanthes species. The chloroplast genome of S. bantonensis is a circular DNA molecule with a total length of 144,591 bp and encodes 84 protein-coding, 8 ribosomes, and 37 transfer RNA genes. The chloroplast genome has a conservative quadripartite structure, including a large single-copy (LSC) region, a small single-copy (SSC) region, and a pair of inverted repeat (IR) regions, with lengths of 92,068 bp, 17,767 bp, and 17,378 bp, respectively. Phylogenetic analysis confirmed that S. bantonensis is closely related to the S. cusia. Compared with other species from Acanthaceae, S. bantonensis has a significantly shortened IR region, suggesting the occurrence of IR contraction events. This study will help future taxonomic, evolutionary, phylogenetic, and bioprospecting studies of the sizeable Strobilanthes genus, which contains over 400 species.

Complete chloroplast genome of Myracrodruon urundeuva and its phylogenetics relationships in Anacardiaceae family

Continuous exploratory use of tree species is threatening the existence of several plants in South America. One of these threatened species is Myracroduron urundeuva, highly exploited due to the high quality and durability of its wood. The chloroplast (cp) has been used for several evolutionary studies as well traceability of timber origin, based on its gene sequences and simple sequence repeats (SSR) variability. Cp genome organization is usually consisting of a large single copy and a small single copy region separated by two inverted repeats regions. We sequenced the complete cp genome from M. urundeuva based on Illumina next-generation sequencing. Our results show that the cp genome is 159,883 bp in size. The 36 SSR identified ranging from mono- to hexanucleotides. Positive selection analysis revealed nine genes related to photosystem, protein synthesis, and DNA replication, and protease are under positive selection. Genome comparison a other Anacardiaceae chloroplast genomes showed great variability in the family. The phylogenetic analysis using complete chloroplast genome sequences of other Anacardiaceae family members showed a close relationship with two other economically important genera, Pistacia and Rhus. These results will help future investigations of timber monitoring and population and evolutionary studies.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-00989-1.

A The complete chloroplast genome and phylogenetic analysis of Bupleurum yinchowense Shan & Yin Li

Bupleurum yinchowense Shan & Yin Li was first described as a new Bupleurum species in 1974, but its classification status has always been disputed. Here, its complete chloroplast genome was provided to resolve this issue. The length of the B. yinchowense chloroplast genome is 155,851 bp and composed of two inverted repeats (IR: 26,307 bp), a large single-copy region (LSC: 85,625 bp), and a small single-copy region (SSC: 17,612 bp). The overall GC content is 37.6%. The chloroplast genome consists of 113 genes, including 79 protein-coding genes, four rRNA genes, and 30 tRNA genes. Phylogenetic analysis suggested that Bupleurum yinchowense holds a distinct phylogenetic position and can be considered as an accepted species.