Complete chloroplast genome sequence and comparative analysis of loblolly pine (Pinus taeda L.) with related species

Comparative and phylogenetic analysis of the complete chloroplast genomes of 10 Artemisia selengensis resources based on high-throughput sequencing

BACKGROUND

Artemisia selengensis, classified within the genus Artemisia of the Asteraceae family, is a perennial herb recognized for its dual utility in culinary and medicinal domains. There are few studies on the chloroplast genome of A. selengensis, and the phylogeographic classification is vague, which makes phylogenetic analysis and evolutionary studies very difficult.

RESULTS

The chloroplast genomes of 10 A. selengensis in this study were highly conserved in terms of gene content, gene order, and gene intron number. The genome lengths ranged from 151,148 to 151,257 bp and were typical of a quadripartite structure with a total GC content of approximately 37.5%. The chloroplast genomes of all species encode 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Due to the contraction and expansion of the inverted repeats (IR), the overlap of ycf1 and ndhF genes occurred at the inverted repeats B (IRB) and short single copy sequence (SSC) boundaries. According to a codon use study, the frequent base in the chloroplast genome of A. selengensis' third codon position was A/T. The number of SSR repeats was 42-44, most of which were single nucleotide A/T repeats. Sequence alignment analysis of the chloroplast genome showed that variable regions were mainly distributed in single copy regions, nucleotide diversity values of 0 to 0.009 were calculated by sliding window analysis, 8 mutation hotspot regions were detected, and coding regions were more conserved than non-coding regions. Analysis of non-synonymous substitution (Ka) and synonymous substitution (Ks) revealed that accD, rps12, petB, and atpF genes were affected by positive selection and no genes were affected by neutral selection. Based on the findings of the phylogenetic analysis, Artemisia selengensis was sister to the genus Artemisia Chrysanthemum and formed a monophyletic group with other Artemisia genera.

CONCLUSIONS

In this research, the present study systematically compared the chloroplast genomic features of A. selengensis and provided important information for the study of the chloroplast genome of A. selengensis and the evolutionary relationships among Asteraceae species.

Plastid genomes provide insights into the phylogeny and chloroplast evolution of the paper daisy tribe Gnaphalieae (Asteraceae)

Chloroplast genomes, as an essential source of phylogenetic information, are increasingly utilized in the evolutionary study of angiosperms. Gnaphalieae is a medium-sized tribe of the sunflower family of Asteraceae, with about 2,100 species in 178 genera distributed in temperate habitats worldwide. There has been considerable progress in our understanding of their phylogenetic evolution using both nuclear and chloroplast sequences, but no focus on chloroplast genomic data. In this study, we performed sequencing, assembly, and annotation of 16 representative chloroplast genomes from all the major lineages of Gnaphalieae. Our results showed that the plastomes exhibited a typical circular tetrad structure with similar genomic structure gene content. But there were differences in genome size, SSRs, and codon usage within the tribe. Phylogenetic analysis revealed Relhania clade is the earliest diverged lineages with the Lasiopogon clade and the Gnaphalium s.s. clade diverged subsequently. The core group includes FLAG clade sister to the HAP and Australasian group. Compared with the outgroup species, chloroplast genome size of the FLAG clade is much reduced whereas those of Australasian, HAP, Gnaphalium s.s., Lasiopogon and Relhania clades are relatively expanded. Insertions and deletions in the intergenic regions associated with repetitive sequence variations are supposed to be the main factor leading to length variations in the chloroplast genomes of Gnaphalieae. The comparative analyses of chloroplast genomes would provide useful implications into understanding the taxonomic and evolutionary history of Gnaphalieae.

Complete chloroplast genome sequence of Pinus tabuliformis var. henryi (Mast.) C.T.Kuan 1983 (Pinaceae)

Pinus tabuliformis var. henryi (Mast.) C.T.Kuan 1983 is an endemic and rare subtropical pine, mainly distributed in central China. In this study, we sequenced the complete chloroplast (cp) genome of P. tabuliformis var. henryi and reported it for the first time. The cp genome was 119,634 base pairs (bp) in total length, including two inverted repeats (IRs, 495 bp), separated by a large single-copy region (LSC, 65,600 bp) and a small single-copy region (SSC, 53,044 bp). There are 114 different genes in the cp genome of P. tabuliformis var. henryi, including 74 protein-coding genes, 36 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of the cp genome was 38.5%. Our phylogenetic analysis of P. tabuliformis var. henryi demonstrated that it was closely related to P. tabuliformis and could be used to identify and analyze its genetic diversity, which was expected to provide new data for taxonomic and phylogenetic studies of Pinus.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

The complete chloroplast genome sequences of four Liparis species (Orchidaceae) and phylogenetic implications

Liparis Richard (Malaxideae, Epidendroideae) is a large and diverse genus of the family Orchidaceae, the taxonomy of which is complicated and controversial. In this study, we sequenced, assembled and analyzed four complete chloroplast genomes of Liparis species including L. kumokiri, L. makinoana, L. pauliana, and L. viridiflora, and evaluated their phylogenetic relationships with related species for the first time. These four chloroplast genomes (size range 153,095 to 158,239 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 83,533-86,752 bp), a small single copy (SSC, 17,938-18,156 bp) and a pair of inverted repeats (IRs, 26,421-26,933 bp). The genomes contain 133 genes, including 87 protein coding genes, 38 tRNAs and 8 rRNA genes. The genome arrangements, gene contents, gene order, long repeats and simple sequence repeats were similar with small differences observed among these four chloroplast genomes. Five highly variable regions including ycf1, ndhA, ndhF, trnQ and trnK were identified from the comparative analysis with other nine related Liparis species, which had the potential to be used as DNA markers for species identification and phylogenetic studies of Liparis species. Phylogenetic analysis based on the complete chloroplast genome sequences strongly supported the polyphyly of Liparis and its further division into three branches. These results provided valuable information to illustrate the complicated taxonomy, phylogeny and evolution process of the Liparis genus.

Assembly, annotation and analysis of the chloroplast genome of the Algarrobo tree Neltuma pallida (subfamily: Caesalpinioideae)

BACKGROUND

Neltuma pallida is a tree that grows in arid soils in northwestern Peru. As a predominant species of the Equatorial Dry Forest ecoregion, it holds significant economic and ecological value for both people and environment. Despite this, the species is severely threatened and there is a lack of genetic and genomic research, hindering the proposal of evidence-based conservation strategies.

RESULTS

In this work, we conducted the assembly, annotation, analysis and comparison of the chloroplast genome of a N. pallida specimen with those of related species. The assembled chloroplast genome has a length of 162,381 bp with a typical quadripartite structure (LSC-IRA-SSC-IRB). The calculated GC content was 35.97%. However, this is variable between regions, with a higher GC content observed in the IRs. A total of 132 genes were annotated, of which 19 were duplicates and 22 contained at least one intron in their sequence. A substantial number of repetitive sequences of different types were identified in the assembled genome, predominantly tandem repeats (> 300). In particular, 142 microsatellites (SSR) markers were identified. The phylogenetic reconstruction showed that N. pallida grouped with the other Neltuma species and with Prosopis cineraria. The analysis of sequence divergence between the chloroplast genome sequences of N. pallida, N. juliflora, P. farcta and Strombocarpa tamarugo revealed a high degree of similarity.

CONCLUSIONS

The N. pallida chloroplast genome was found to be similar to those of closely related species. With a size of 162,831 bp, it had the classical chloroplast quadripartite structure and GC content of 35.97%. Most of the 132 identified genes were protein-coding genes. Additionally, over 800 repetitive sequences were identified, including 142 SSR markers. In the phylogenetic analysis, N. pallida grouped with other Neltuma spp. and P. cineraria. Furthermore, N. pallida chloroplast was highly conserved when compared with genomes of closely related species. These findings can be of great potential for further diversity studies and genetic improvement of N. pallida.

DNA:RNA Hybrids Are Major Dinoflagellate Minicircle Molecular Types

Peridinin-containing dinoflagellate plastomes are predominantly encoded in nuclear genomes, with less than 20 essential chloroplast proteins carried on "minicircles". Each minicircle generally carries one gene and a short non-coding region (NCR) with a median length of approximately 400-1000 bp. We report here differential nuclease sensitivity and two-dimensional southern blot patterns, suggesting that dsDNA minicircles are in fact the minor forms, with substantial DNA:RNA hybrids (DRHs). Additionally, we observed large molecular weight intermediates, cell-lysate-dependent NCR secondary structures, multiple bidirectional predicted ssDNA structures, and different southern blot patterns when probed with different NCR fragments. In silico analysis suggested the existence of substantial secondary structures with inverted repeats (IR) and palindrome structures within the initial ~650 bp of the NCR sequences, in accordance with conversion event(s) outcomes with PCR. Based on these findings, we propose a new transcription-templating-translation model, which is associated with cross-hopping shift intermediates. Since dinoflagellate chloroplasts are cytosolic and lack nuclear envelope breakdown, the dynamic DRH minicircle transport could have contributed to the spatial-temporal dynamics required for photosystem repair. This represents a paradigm shift from the previous understanding of "minicircle DNAs" to a "working plastome", which will have significant implications for its molecular functionality and evolution.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

The chloroplast genome of black pepper (Piper nigrum L.) and its comparative analysis with related Piper species

Piper nigrum, also known as black pepper, is an economically and ecologically important crop of the genus Piper. It has been titled as the king of spices due to its wide consumption throughout the world. In the present investigation, the chloroplast genome of P. nigrum has been assembled from a whole genome sequence by integrating the short and long reads generated through Illumina and PacBio platforms, respectively. The chloroplast genome was observed to be 161,522 bp in size, having a quadripartite structure with a large single copy (LSC) region of 89,153 bp and a small single copy (SSC) region of 18,255 bp separated by a copy of inverted repeats (IRs), each 27,057 bp in length. Taking into consideration all the duplicated genes, a total of 131 genes were observed, which included 81 protein-coding genes, 37 tRNAs, 4 rRNAs, and 1 pseudogene. Individually, the LSC region consisted of 83 genes, the SSC region had 13 genes, and 18 genes were present in each IR region. Additionally, 216 SSRs were detected and 11 of these were validated through amplification in 12 species of Piper. The features of the chloroplast genome have been compared with those of the genus Piper. Our results provide useful insights into evolutionary and molecular studies of black pepper which will contribute to its further genetic improvement and breeding.

The complete chloroplast genome sequences of three Broussonetia species and comparative analysis within the Moraceae

Background

Species of Broussonetia (family Moraceae) are commonly used to make textiles and high-grade paper. The distribution of Broussonetia papyrifera L. is considered to be related to the spread and location of humans. The complete chloroplast (cp) genomes of B. papyrifera, Broussonetia kazinoki Sieb., and Broussonetia kaempferi Sieb. were analyzed to better understand the status and evolutionary biology of the genus Broussonetia.

Methods

The cp genomes were assembled and characterized using SOAPdenovo2 and DOGMA. Phylogenetic and molecular dating analysis were performed using the concatenated nucleotide sequences of 35 species in the Moraceae family and were based on 66 protein-coding genes (PCGs). An analysis of the sequence divergence (pi) of each PCG among the 35 cp genomes was conducted using DnaSP v6. Codon usage indices were calculated using the CodonW program.

Results

All three cp genomes had the typical land plant quadripartite structure, ranging in size from 160,239 bp to 160,841 bp. The ribosomal protein L22 gene (RPL22) was either incomplete or missing in all three Broussonetia species. Phylogenetic analysis revealed two clades. Clade 1 included Morus and Artocarpus, whereas clade 2 included the other seven genera. Malaisia scandens Lour. was clustered within the genus Broussonetia. The differentiation of Broussonetia was estimated to have taken place 26 million years ago. The PCGs' pi values ranged from 0.0005 to 0.0419, indicating small differences within the Moraceae family. The distribution of most of the genes in the effective number of codons plot (ENc-plot) fell on or near the trend line; the slopes of the trend line of neutrality plots were within the range of 0.0363-0.171. These results will facilitate the identification, taxonomy, and utilization of the Broussonetia species and further the evolutionary studies of the Moraceae family.

Comparative chloroplast genome and phylogenetic analyses of Chinese Polyspora

Polyspora Sweet (Theaceae) are winter ornamental landscape plants native to southern and southeastern Asia, some of which have medicinal value. The chloroplast (cp) genome data of Polyspora are scarce, and the gene evolution and interspecific relationship are still unclear. In this study, we sequenced and annotated Polyspora chrysandra cp genome and combined it with previously published genomes for other Chinese Polyspora species. The results showed that cp genomes of six Chinese Polyspora varied in length between 156,452 bp (P. chrysandra) and 157,066 bp (P. speciosa), but all contained 132 genes, with GC content of 37.3%, and highly similar genes distribution and codon usage. A total of eleven intergenic spacer regions were found having the highest levels of divergence, and eight divergence hotspots were identified as molecular markers for Phylogeography and genetic diversity studies in Polyspora. Gene selection pressure suggested that five genes were subjected to positive selection. Phylogenetic relationships among Polyspora species based on the complete cp genomes were supported strongly, indicating that the cp genomes have the potential to be used as super barcodes for further analysis of the phylogeny of the entire genus. The cp genomes of Chinese Polyspora species will provide valuable information for species identification, molecular breeding and evolutionary analysis of genus Polyspora.

Fortunella venosa (Champ. ex Benth.) C. C. Huang and F. hindsii (Champ. ex Benth.) Swingle as Independent Species: Evidence From Morphology and Molecular Systematics and Taxonomic Revision of Fortunella (Rutaceae)

Recently, the systematic status of Fortunella Swingle and its taxonomy has attracted much attention. Flora of China incorporates Fortunella into Citrus Linn. and treats all species of the traditional Fortunella as one species, namely Citrus japonica (Thunb.) Swingle. Furthermore, F. venosa (Champ. ex Benth.) C. C. Huang and F. hindsii (Champ. ex Benth.) Swingle are currently considered as synonyms of C. japonica. In this paper, morphological, palynological, and phylogenetic analyses were used to systematically explore the taxonomic status of traditional Fortunella. The key morphological features that differed among the Fortunella species were the leaf and the petiole hence could be key in its taxonomic classification of the species. Additionally, pollen morphological analysis based on the pollen size, germination grooves, polar, and equatorial axes also supported the separation of the species. The results of the phylogenetic analysis showed that each of the three species clustered separately, hence strongly supporting the conclusion of independent species. In addition, the phylogenetic analysis showed that the two genera clustered closely together hence our results support the incorporation of Fortunella into Citrus. Based on the above, this article has revised the classification of the traditional Fortunella and determined that this genus has three species, namely; F. venosa, F. hindsii, and F. japonica. F. venosa and F. hindsii are placed in the Citrus as separate species, and their species names still use the previous specific epithet. The revised scientific names of the new combinations of F. venosa and F. hindsii are as follows: Citrus venosa (Champ. ex Benth.) K. M. Liu, X. Z. Cai, and G. W. Hu, comb. nov. and Citrus hindsii (Champ. ex Benth.) K. M. Liu, G. W. Hu, and X. Z. Cai, comb. nov. F. venosa is the original species of Fortunella, F. venosa and F. hindsii are both listed as the second-class key protected wild plants in China. Therefore, the establishment of the taxonomic status of F. venosa and F. hindsii not only deepens our understanding, importance, and the complexity of the systematic classification of Fortunella, but is also significant for global biodiversity conservation, genetic resources for breeding purposes, and population genetics.

Comparative analysis of complete Ilex (Aquifoliaceae) chloroplast genomes: insights into evolutionary dynamics and phylogenetic relationships

BACKGROUND

Ilex (Aquifoliaceae) are of great horticultural importance throughout the world for their foliage and decorative berries, yet a dearth of genetic information has hampered our understanding of phylogenetic relationships and evolutionary history. Here, we compare chloroplast genomes from across Ilex and estimate phylogenetic relationships.

RESULTS

We sequenced the chloroplast genomes of seven Ilex species and compared them with 34 previously published Ilex plastomes. The length of the seven newly sequenced Ilex chloroplast genomes ranged from 157,182 bp to 158,009 bp, and contained a total of 118 genes, including 83 protein-coding, 31 rRNA, and four tRNA genes. GC content ranged from 37.6 to 37.69%. Comparative analysis showed shared genomic structures and gene rearrangements. Expansion and contraction of the inverted repeat regions at the LSC/IRa and IRa/SSC junctions were observed in 22 and 26 taxa, respectively; in contrast, the IRb boundary was largely invariant. A total of 2146 simple sequence repeats and 2843 large repeats were detected in the 41 Ilex plastomes. Additionally, six genes (psaC, rbcL, trnQ, trnR, trnT, and ycf1) and two intergenic spacer regions (ndhC-trnV and petN-psbM) were identified as hypervariable, and thus potentially useful for future phylogenetic studies and DNA barcoding. We recovered consistent phylogenetic relationships regardless of inference methodology or choice of loci. We recovered five distinct, major clades, which were inconsistent with traditional taxonomic systems.

CONCLUSION

Our findings challenge traditional circumscriptions of the genus Ilex and provide new insights into the evolutionary history of this important clade. Furthermore, we detail hypervariable and repetitive regions that will be useful for future phylogenetic and population genetic studies.

Comparative analyses of chloroplast genomes in 'Red Fuji' apples: low rate of chloroplast genome mutations

BACKGROUND

Fuji is a vital apple cultivar, and has been propagated clonally for nearly a century. The chloroplast genome variation of Fuji apples in China has not been investigated.

METHODS

This study used next-generation high-throughput sequencing and bioinformatics to compare and analyze the chloroplast genome of 24 Red Fuji varieties from nine regions in China.

RESULTS

The results showed that the 24 chloroplast genomes were highly conserved in genome size, structure, and organization. The length of the genomes ranged from 160,063 to 160,070 bp, and the GC content was 36.6%. Each of the 24 chloroplast genomes encoded 131 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The results of repeat sequence detection were consistent; the most common sequence was forward repeats (53.1%), and the least common sequence was complementary repeats (4.1%). The chloroplast genome sequence of Red Fuji was highly conserved. Two indels were detected, but the PI value was 0, and there were no SNP loci. The chloroplast genome variation rate of Red Fuji was low.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.‐Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single‐copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single‐copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132–133 genes, including 85–87 protein‐coding genes, 37–38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37–43 long repeats and 175–178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

The first complete chloroplast genome of Vicatia thibetica de Boiss.: genome features, comparative analysis, and phylogenetic relationships

Vicatia thibetica de Boiss.: a herb in the family Apiaceae, has been used for over a hundred years as an essential medicinal and edible plant in the Bai ethnic group of Dali City. However, due to the lack of study on plastid genomes of V. thibetica, studies of comparison and phylogeny with other related species remain scarce. In the current study, we assembled, annotated, and characterized the entire chloroplast (cp) genome of V. thibetica through high-throughput sequencing for the first time, compared with published whole chloroplast genomes from the same family. A phylogenetic analysis of the chloroplast genome has also been performed. The whole chloroplast genome of V. thibetica was 145,796 in size and consisted of a large single-copy region (LSC; 92,186 bp), a small single-copy region (SSC; 17,452 bp), and a pair of inverted repeat regions (IRs; 18,079 bp) forming a circular quadripartite structure. Annotation resulted in 128 genes, including 84 protein-coding genes (PCGs), 35 transfer RNA genes (tRNAs), eight ribosomal genes (rRNAs), and one pseudogene. Repeat sequence analysis displayed V. thibetica plastid genome contains 75 simple repeats, 37 long repeats, and 29 tandem repeats. Compared with the cp genome of other Apiaceae species, a common feature was that the IR regions of the genome were more conservative compared to the LSC and SSC regions. Highly variable hotspots included rps16, ndhC-trnV-UAC, clpP, ycf1, and ndhB in the genomes, which supply valuable molecular markers for phylogeny, identification, and classification in the Apiaceae family. The results of phylogenetic analysis strongly supported the genus Vicatia as an independent genus in the family Apiaceae, in which the closest affinities to the related species of Angelica, Peucedanum, and Ligusticum were observed. In conclusion, the first chloroplast genome of Vicatia reported in this study may  improve our understanding of phylogenetic relationship of different genera of Apiaceae. In addition, the current data will be valuable as chloroplast genomic resource for species identification and population genetics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01154-y.

cpSSR and High-Resolution Melting Analysis (HRM) for Pinus pseudostrobus Lindl. Variety Genotyping and Discrimination

The unambiguous identification of varieties within the Pseudostrobus complex is a key step to facilitate tree selection and monitoring in the wild as well as in plantations. Molecular tools provide a powerful approach for species delimitation; however, the use of DNA barcodes in this group has met limited success due to widespread haplotype sharing from lineage sorting, hybridization and introgression. Here, we evaluate the utility of real-time PCR coupled with high-resolution melting (HRM) to discriminate among Pinus pseudostrobus Lindl. var. pseudostrobus, apulcensis and oaxacana, from wild populations in central and southern Mexico, using chloroplast DNA sequence variants located within the clpP, ycf2, trnL(UAA)–trnT(UGU) and trnI(CAU)–trnF(GAA) loci. The markers ycf2/trnL(UAA)–trnT(UGU) produced clear melting patterns that separated the varieties pseudostrobus and oaxacana from type var. apulcensis, whereas clpP discriminated over 60% of var. oaxacana individuals. This assay underlines the usefulness of these less-used DNA regions as potential biological markers and exhibits the effect of geography on allele distribution and the likely presence of hybrids among the species and varieties.

Comparative chloroplast genome analysis of medicinally important Veratrum (Melanthiaceae) in China: Insights into genomic characterization and phylogenetic relationships

Members of Veratrum are perennial herbs widely used in traditional Chinese medicine to induce vomiting, resolve blood stasis and relieve pain. However, the intrageneric classification and phylogenetic relationships within Veratrum have long been controversial due to the complexity of morphological variations and lack of high-resolution molecular markers. In this study, we reevaluated the infrageneric relationships with the genus Veratrum using complete chloroplast genome sequence data. Herein, the complete cp genomes of ten species of Veratrum were newly sequenced and characterized. The complete cp genomes of ten species of Veratrum had the typical quadripartite structure, ranging from 151,597 bp to 153,711 bp in size and comprising a total of 135 genes. The structure of Veratrum cp genomes (i.e., gene order, content, and genome components) was highly similar across species. The number of simple sequence repeats (SSRs) ranged from 63 to 78, and of long repeats ranged from 31 to 35. Eight highly divergent regions (ndhF, psbC-psbZ, psbK-psbI, rpoB-trnC_GCA, trnK_UUU-trnQ_UUG, trnS_GCU-trnG_UCC, trnT_UGU-trnL_UAA and ycf1) were identified and are potentially useful for the DNA barcoding of Veratrum. Phylogenetic analysis among 29 taxa based on cp genomes, total genes, protein-coding genes and intergenic regions strongly supported the monophyly of Veratrum. The circumscription and relationships of the infrageneric taxa of Veratrum were well-presented with great resolution. These results will facilitate the identification, taxonomy, and utilization of Veratrum plants as well as the evolutionary studies of Melanthiaceae.

Comparative Genomics Analysis of Repetitive Elements in Ten Gymnosperm Species: "Dark Repeatome" and Its Abundance in Conifer and Gnetum Species

Repetitive elements (RE) and transposons (TE) can comprise up to 80% of some plant genomes and may be essential for regulating their evolution and adaptation. The “repeatome” information is often unavailable in assembled genomes because genomic areas of repeats are challenging to assemble and are often missing from final assembly. However, raw genomic sequencing data contain rich information about RE/TEs. Here, raw genomic NGS reads of 10 gymnosperm species were studied for the content and abundance patterns of their “repeatome”. We utilized a combination of alignment on databases of repetitive elements and de novo assembly of highly repetitive sequences from genomic sequencing reads to characterize and calculate the abundance of known and putative repetitive elements in the genomes of 10 conifer plants: Pinus taeda, Pinus sylvestris, Pinus sibirica, Picea glauca, Picea abies, Abies sibirica, Larix sibirica, Juniperus communis, Taxus baccata, and Gnetum gnemon. We found that genome abundances of known and newly discovered putative repeats are specific to phylogenetically close groups of species and match biological taxa. The grouping of species based on abundances of known repeats closely matches the grouping based on abundances of newly discovered putative repeats (kChains) and matches the known taxonomic relations.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Variation and Evolution of the Whole Chloroplast Genomes of Fragaria spp. (Rosaceae)

Species identification is vital for protecting species diversity and selecting high-quality germplasm resources. Wild Fragaria spp. comprise rich and excellent germplasm resources; however, the variation and evolution of the whole chloroplast (cp) genomes in the genus Fragaria have been ignored. In the present study, 27 complete chloroplast genomes of 11 wild Fragaria species were sequenced using the Illumina platform. Then, the variation among complete cp genomes of Fragaria was analyzed, and phylogenetic relationships were reconstructed from those genome sequences. There was an overall high similarity of sequences, with some divergence. According to analysis with mVISTA, non-coding regions were more variable than coding regions. Inverted repeats (IRs) were observed to contract or expand to different degrees, which resulted in different sizes of cp genomes. Additionally, five variable loci, trnS-trnG, trnR-atpA, trnC-petN, rbcL-accD, and psbE-petL, were identified that could be used to develop DNA barcoding for identification of Fragaria species. Phylogenetic analyses based on the whole cp genomes supported clustering all species into two groups (A and B). Group A species were mainly distributed in western China, while group B contained several species from Europe and Americas. These results support allopolyploid origins of the octoploid species F. chiloensis and F. virginiana and the tetraploid species F. moupinensis and F. tibetica. The complete cp genomes of these Fragaria spp. provide valuable information for selecting high-quality Fragaria germplasm resources in the future.

Complete Chloroplast Genomes of Fagus sylvatica L. Reveal Sequence Conservation in the Inverted Repeat and the Presence of Allelic Variation in NUPTs

Growing amounts of genomic data and more efficient assembly tools advance organelle genomics at an unprecedented scale. Genomic resources are increasingly used for phylogenetic analyses of many plant species, but are less frequently used to investigate within-species variability and phylogeography. In this study, we investigated genetic diversity of Fagus sylvatica, an important broadleaved tree species of European forests, based on complete chloroplast genomes of 18 individuals sampled widely across the species distribution. Our results confirm the hypothesis of a low cpDNA diversity in European beech. The chloroplast genome size was remarkably stable (158,428 ± 37 bp). The polymorphic markers, 12 microsatellites (SSR), four SNPs and one indel, were found only in the single copy regions, while inverted repeat regions were monomorphic both in terms of length and sequence, suggesting highly efficient suppression of mutation. The within-individual analysis of polymorphisms showed >9k of markers which were proportionally present in gene and non-gene areas. However, an investigation of the frequency of alternate alleles revealed that the source of this diversity originated likely from nuclear-encoded plastome remnants (NUPTs). Phylogeographic and Mantel correlation analysis based on the complete chloroplast genomes exhibited clustering of individuals according to geographic distance in the first distance class, suggesting that the novel markers and in particular the cpSSRs could provide a more detailed picture of beech population structure in Central Europe.

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.

Complete Chloroplast Genomes of Acanthochlamys bracteata (China) and Xerophyta (Africa) (Velloziaceae): Comparative Genomics and Phylogenomic Placement

Acanthochlamys P.C. Kao is a Chinese endemic monotypic genus, whereas Xerophyta Juss. is a genus endemic to Africa mainland, Arabian Peninsula and Madagascar with ca.70 species. In this recent study, the complete chloroplast genome of Acanthochlamys bracteata was sequenced and its genome structure compared with two African Xerophyta species (Xerophyta spekei and Xerophyta viscosa) present in the NCBI database. The genomes showed a quadripartite structure with their sizes ranging from 153,843 bp to 155,498 bp, having large single-copy (LSC) and small single-copy (SSC) regions divided by a pair of inverted repeats (IR regions). The total number of genes found in A. bracteata, X. spekei and X. viscosa cp genomes are 129, 130, and 132, respectively. About 50, 29, 28 palindromic, forward and reverse repeats and 90, 59, 53 simple sequence repeats (SSRs) were found in the A. bracteata, X. spekei, and X. viscosa cp genome, respectively. Nucleotide diversity analysis in all species was 0.03501, Ka/Ks ratio average score was calculated to be 0.26, and intergeneric K2P value within the Order Pandanales was averaged to be 0.0831. Genomic characterization was undertaken by comparing the genomes of the three species of Velloziaceae and it revealed that the coding regions were more conserved than the non-coding regions. However, key variations were noted mostly at the junctions of IRs/SSC regions. Phylogenetic analysis suggests that A. bracteata species has a closer genetic relationship to the genus Xerophyta. The present study reveals the complete chloroplast genome of A. bracteata and gives a genomic comparative analysis with the African species of Xerophyta. Thus, can be useful in developing DNA markers for use in the study of genetic variabilities and evolutionary studies in Velloziaceae.

Mangrove tree (Avicennia marina): insight into chloroplast genome evolutionary divergence and its comparison with related species from family Acanthaceae

Avicennia marina (family Acanthaceae) is a halotolerant woody shrub that grows wildly and cultivated in the coastal regions. Despite its importance, the species suffers from lack of genomic datasets to improve its taxonomy and phylogenetic placement across the related species. Here, we have aimed to sequence the plastid genome of A. marina and its comparison with related species in family Acanthaceae. Detailed next-generation sequencing and analysis showed a complete chloroplast genome of 150,279 bp, comprising 38.6% GC. Genome architecture is quadripartite revealing large single copy (82,522 bp), small single copy (17,523 bp), and pair of inverted repeats (25,117 bp). Furthermore, the genome contains 132 different genes, including 87 protein-coding genes, 8 rRNA, 37 tRNA genes, and 126 simple sequence repeats (122 mononucleotide, 2 dinucleotides, and 2 trinucleotides). Interestingly, about 25 forward, 15 reversed and 14 palindromic repeats were also found in the A. marina. High degree synteny was observed in the pairwise alignment with related genomes. The chloroplast genome comparative assessment showed a high degree of sequence similarity in coding regions and varying divergence in the intergenic spacers among ten Acanthaceae species. The pairwise distance showed that A. marina exhibited the highest divergence (0.084) with Justicia flava and showed lowest divergence with Aphelandra knappiae (0.059). Current genomic datasets are a valuable resource for investigating the population and evolutionary genetics of family Acanthaceae members' specifically A. marina and related species.

Complete Chloroplast Genome of Argania spinosa: Structural Organization and Phylogenetic Relationships in Sapotaceae

Argania spinosa (Sapotaceae), an important endemic Moroccan oil tree, is a primary source of argan oil, which has numerous dietary and medicinal proprieties. The plant species occupies the mid-western part of Morocco and provides great environmental and socioeconomic benefits. The complete chloroplast (cp) genome of A. spinosa was sequenced, assembled, and analyzed in comparison with those of two Sapotaceae members. The A. spinosa cp genome is 158,848 bp long, with an average GC content of 36.8%. The cp genome exhibits a typical quadripartite and circular structure consisting of a pair of inverted regions (IR) of 25,945 bp in length separating small single-copy (SSC) and large single-copy (LSC) regions of 18,591 and 88,367 bp, respectively. The annotation of A. spinosa cp genome predicted 130 genes, including 85 protein-coding genes (CDS), 8 ribosomal RNA (rRNA) genes, and 37 transfer RNA (tRNA) genes. A total of 44 long repeats and 88 simple sequence repeats (SSR) divided into mononucleotides (76), dinucleotides (7), trinucleotides (3), tetranucleotides (1), and hexanucleotides (1) were identified in the A. spinosa cp genome. Phylogenetic analyses using the maximum likelihood (ML) method were performed based on 69 protein-coding genes from 11 species of Ericales. The results confirmed the close position of A. spinosa to the Sideroxylon genus, supporting the revisiting of its taxonomic status. The complete chloroplast genome sequence will be valuable for further studies on the conservation and breeding of this medicinally and culinary important species and also contribute to clarifying the phylogenetic position of the species within Sapotaceae.

Molecular structure and phylogenetic analyses of the complete chloroplast genomes of three original species of Pyrrosiae Folium

Pyrrosia petiolosa, Pyrrosia lingua and Pyrrosia sheareri are recorded as original plants of Pyrrosiae Folium (PF) and commonly used as Chinese herbal medicines. Due to the similar morphological features of PF and its adulterants, common DNA barcodes cannot accurately distinguish PF species. Knowledge of the chloroplast (cp) genome is widely used in species identification, molecular marker and phylogenetic analyses. Herein, we determined the complete cp genomes of three original species of PF via high-throughput sequencing technologies. The three cp genomes exhibited a typical quadripartite structure with sizes ranging from 158 165 to 163 026 bp. The cp genomes of P. petiolosa and P. lingua encoded 130 genes, whilst that of P. sheareri encoded 131 genes. The complete cp genomes were compared, and five highly divergent regions of petA-psbJ, matK-rps16, ndhC-trnM, psbM-petN and psaC-ndhE were screened as potential DNA barcodes for identification of Pyrrosia genus species. The phylogenetic tree we obtained indicated that P. petiolosa and P. lingua are clustered in a single clade and, thus, share a close relationship. This study provides invaluable information for further studies on the species identification, taxonomy and phylogeny of Pyrrosia genus species.

A comparative analysis of the complete chloroplast genomes of three Chrysanthemum boreale strains

BACKGROUND

Chrysanthemum boreale Makino (Anthemideae, Asteraceae) is a plant of economic, ornamental and medicinal importance. We characterized and compared the chloroplast genomes of three C. boreale strains. These were collected from different geographic regions of Korea and varied in floral morphology.

METHODS

The chloroplast genomes were obtained by next-generation sequencing techniques, assembled de novo, annotated, and compared with one another. Phylogenetic analysis placed them within the Anthemideae tribe.

RESULTS

The sizes of the complete chloroplast genomes of the C. boreale strains were 151,012 bp (strain 121002), 151,098 bp (strain IT232531) and 151,010 bp (strain IT301358). Each genome contained 80 unique protein-coding genes, 4 rRNA genes and 29 tRNA genes. Comparative analyses revealed a high degree of conservation in the overall sequence, gene content, gene order and GC content among the strains. We identified 298 single nucleotide polymorphisms (SNPs) and 106 insertions/deletions (indels) in the chloroplast genomes. These variations were more abundant in non-coding regions than in coding regions. Long dispersed repeats and simple sequence repeats were present in both coding and noncoding regions, with greater frequency in the latter. Regardless of their location, these repeats can be used for molecular marker development. Phylogenetic analysis revealed the evolutionary relationship of the species in the Anthemideae tribe. The three complete chloroplast genomes will be valuable genetic resources for studying the population genetics and evolutionary relationships of Asteraceae species.

Phylogenetic analysis of Fritillaria cirrhosa D. Don and its closely related species based on complete chloroplast genomes

Fritillaria cirrhosa D. Don, whose bulb is used in a well-known traditional Chinese medicine to relieve cough and eliminate phlegm, is one of the most important medicinal plants of Fritillaria L. The species is widely distributed among the alpine regions in southwestern China and possesses complex morphological variations in different distributions. A series of newly related species were reported, based on obscure morphological differences. As a result, F. cirrhosa and its closely related species constitute a taxonomically complex group. However, it is difficult to accurately identify these species and reveal their phylogenetic relationships using traditional taxonomy. Molecular markers and gene fragments have been adopted but they are not able to afford sufficient phylogenetic resolution in the genus. Here, we report the complete chloroplast genome sequences of F. cirrhosa and its closely related species using next generation sequencing (NGS) technology. Eight plastid genomes ranged from 151,058 bp to 152,064 bp in length and consisted of 115 genes. Gene content, gene order, GC content, and IR/SC boundary structures were highly similar among these genomes. SSRs and five large repeat sequences were identified and the total number of them ranged from 73 to 79 and 63 to 75, respectively. Six highly divergent regions were successfully identified that could be used as potential genetic markers of Fritillaria. Phylogenetic analyses revealed that eight Fritillaria species were clustered into three clades with strong supports and F. cirrhosa was closely related to F. przewalskii and F. sinica. Overall, this study indicated that the complete chloroplast genome sequence was an efficient tool for identifying species in taxonomically complex groups and exploring their phylogenetic relationships.

Complete Chloroplast Genome of Pinus densiflora Siebold & Zucc. and Comparative Analysis with Five Pine Trees

Pinus densiflora (Korean red pine) is widely distributed in East Asia and considered one of the most important species in Korea. In this study, the complete chloroplast genome of P. densiflora was sequenced by combining the advantages of Oxford Nanopore MinION and Illumina MiSeq. The sequenced genome was then compared with that of a previously published conifer plastome. The chloroplast genome was found to be circular and comprised of a quadripartite structure, including 113 genes encoding 73 proteins, 36 tRNAs and 4 rRNAs. It had short inverted repeat regions and lacked ndh gene family genes, which is consistent with other Pinaceae species. The gene content of P. densiflora was found to be most similar to that of P. sylvestris. The newly attempted sequencing method could be considered an alternative method for obtaining accurate genetic information, and the chloroplast genome sequence of P. densiflora revealed in this study can be used in the phylogenetic analysis of Pinus species.

Araucaria angustifolia chloroplast genome sequence and its relation to other Araucariaceae

Araucaria angustifolia is endemic to southern Brazil. Known as Brazilian pine, A. angustifolia is the only native conifer species with economic and social relevance in this country. Due to massive exploitation, it has suffered a significant population decline and currently is classified as critically endangered. This encouraged the scientific community to investigate genetic features in Brazilian pine to increase resources for management and preservation. In this work, RNA-Seq data was used to determine the complete nucleotide sequence of the A. angustifolia chloroplast genome (cpDNA). The cpDNA is 146,203 bp in length and contains 122 genes, including 80 protein-coding genes, 5 ribosomal RNA genes, and 37 tRNA genes. Coding regions comprise 45.02%, 4.96% correspond to rRNAs and tRNAs, and 50.02% of the genome encompasses non-coding regions. Genes found in the inverted repeat (IR) are present as single copy, with exception of the rrn5 and trnI-CAU loci. The typical LSC, SSC, IRa and IRb organization reported in several land-plant groups is not present in A. angustifolia cpDNA. Phylogenetic analyses using Bayesian and Maximum Likelihood methods clustered A. angustifolia in the Araucariaceae family, with A. heterophylla and A. columnaris as congeneric species. The screening of A. angustifolia cpDNA reveled 100 SSRs, 14 of them corresponding to tetrapolymer loci.

The complete chloroplast genome sequences of three Spondias species reveal close relationship among the species

This study reports the complete chloroplast sequences of three Spondias species. The genome sequences were obtained for Spondias tuberosa, Spondias bahienses, and Spondias mombin using the Illumina sequencing technology by a combination of de novo methods and a reference-guided assembly using Sapindus mukorossi as reference. The genomes of S. tuberosa, S. bahiensis, and S. mombin had 162,036, 162,218, and 162,302 bp, respectively. The coding regions exhibited 130 genes, including 34-35 tRNAs and 4 rRNAs. The results revealed synteny among the genomes, with high conservation in the gene order and content and CG content. The inverted repeat regions (IRA and IRB) and the large and small single copies were very similar among the three genomes. The phylogenomic analysis reported similar topologies as that of previous studies, which used partial chloroplast, wherein S. mombin was the first diverging lineage, while S. tuberosa and S. bahiensis were derived, indicating that the phylogenetic analysis using partial or complete genome produces similar results. In summary, (1) we presented the first complete chloroplast genome for the genus Spondias, (2) phylogenies analyzed using the complete chloroplast genomes revealed a robust phylogenetic topology for Spondias, and (3) gene order, content, and orientation in Spondias are highly conserved.