Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla.

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20-30 bp and 10-20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species.

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

Assembly and comparative analysis of the complete mitochondrial and chloroplast genome of Cyperus stoloniferus (Cyperaceae), a coastal plant possessing saline-alkali tolerance.

BACKGROUND: Cyperus stoloniferus is an important species in coastal ecosystems and possesses economic and ecological value. To elucidate the structural characteristics, variation, and evolution of the organelle genome of C. stoloniferus, we sequenced, assembled, and compared its mitochondrial and chloroplast genomes. RESULTS: We assembled the mitochondrial and chloroplast genomes of C. stoloniferus. The total length of the mitochondrial genome (mtDNA) was 927,413 bp, with a GC content of 40.59%. It consists of two circular DNAs, including 37 protein-coding genes (PCGs), 22 tRNAs, and five rRNAs. The length of the chloroplast genome (cpDNA) was 186,204 bp, containing 93 PCGs, 40 tRNAs, and 8 rRNAs. The mtDNA and cpDNA contained 81 and 129 tandem repeats, respectively, and 346 and 1,170 dispersed repeats, respectively, both of which have 270 simple sequence repeats. The third high-frequency codon (RSCU > 1) in the organellar genome tended to end at A or U, whereas the low-frequency codon (RSCU < 1) tended to end at G or C. The RNA editing sites of the PCGs were relatively few, with only 9 and 23 sites in the mtDNA and cpDNA, respectively. A total of 28 mitochondrial plastid DNAs (MTPTs) in the mtDNA were derived from cpDNA, including three complete trnT-GGU, trnH-GUG, and trnS-GCU. Phylogeny and collinearity indicated that the relationship between C. stoloniferus and C. rotundus are closest. The mitochondrial rns gene exhibited the greatest nucleotide variability, whereas the chloroplast gene with the greatest nucleotide variability was infA. Most PCGs in the organellar genome are negatively selected and highly evolutionarily conserved. Only six mitochondrial genes and two chloroplast genes exhibited Ka/Ks > 1; in particular, atp9, atp6, and rps7 may have undergone potential positive selection. CONCLUSION: We assembled and validated the mtDNA of C. stoloniferus, which contains a 15,034 bp reverse complementary sequence. The organelle genome sequence of C. stoloniferus provides valuable genomic resources for species identification, evolution, and comparative genomic research in Cyperaceae.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination.

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention.

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

Exploring phylogeny and genomic vulnerability of Melastoma (Melastomataceae) endemic to a World Natural Heritage site, the Bonin Islands.

Melastoma (Melastomataceae) includes ca. 100 species across tropical to subtropical regions of Asia and Oceania. The Bonin Islands harbor three endemic taxa: M. tetramerum, M. tetramerum var. pentapetalum, and M. candidum var. alessandrense. Of these, M. tetramerum is critically endangered and faces near extinction in the wild. This study investigates the phylogenetic relationships among these endemic Melastoma species in the Bonin Islands based on the whole chloroplast genome and nuclear SNPs. The results revealed that M. candidum var. alessandrense was placed in the clade of the widespread East Asian M. candidum and has a distinct evolutionary origin from the other two taxa. The population genomics analyses (heterozygosity, rates of deleterious mutations, and numbers and lengths of runs of homozygosity) indicated lower genetic diversity and more vulnerable genomes of endemic Melastoma, especially M. tetramerum var. pentapetalum. M. tetramerum var. pentapetalum is not a target of any protection programs, however, conservation plans might be required for this variety because M. tetramerum var. pentapetalum would have a more vulnerable genome than M. tetramerum, which faces near extinction in the wild. This information can facilitate the development of effective conservation strategies in a precautionary way that anticipates imminent threats to the survival of the species.

Chloroplast Genomes Evolution and Phylogenetic Relationships of Caragana species.

Caragana sensu lato (s.l.) includes approximately 100 species that are mainly distributed in arid and semi-arid regions. Caragana species are ecologically valuable for their roles in windbreaking and sand fixation. However, the taxonomy and phylogenetic relationships of the genus Caragana are still unclear. In this study, we sequenced and assembled the chloroplast genomes of representative species of Caragana and reconstructed robust phylogenetic relationships at the section level. The Caragana chloroplast genome has lost the inverted repeat region and wascategorized in the inverted repeat loss clade (IRLC). The chloroplast genomes of the eight species ranged from 128,458 bp to 135,401 bp and contained 110 unique genes. All the Caragana chloroplast genomes have a highly conserved structure and gene order. The number of long repeats and simple sequence repeats (SSRs) showed significant variation among the eight species, indicating heterogeneous evolution in Caragana. Selective pressure analysis of the genes revealed that most of the protein-coding genes evolved under purifying selection. The phylogenetic analyses indicated that each section forms a clade, except the section Spinosae, which was divided into two clades. This study elucidated the evolution of the chloroplast genome within the widely distributed genus Caragana. The detailed information obtained from this study can serve as a valuable resource for understanding the molecular dynamics and phylogenetic relationships within Caragana.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage.

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

Complete chloroplast genomes of eight Delphinium taxa (Ranunculaceae) endemic to Xinjiang, China: insights into genome structure, comparative analysis, and phylogenetic relationships.

BACKGROUND: Delphinium L. represents a taxonomically intricate genus of significant phylogenetic and economic importance in Ranunculaceae. Despite the existence of few chloroplast genome datasets, a comprehensive understanding of genome structures and selective pressures within the genus remains unknown. Furthermore, several taxa in this genus are exclusively found in Xinjiang, China, a region renowned for its distribution and diversity of Chinese and Central Asian Delphinium species. Therefore, investigating the features of chloroplast genomes in this area will provide valuable insights into the evolutionary processes and phylogenetic relationships of the genus. RESULTS: In this study, the eight newly completed chloroplast genomes are examined, ranging in length from 153,979 bp to 154,284 bp. Alongside these, analysing six previously reported taxa re-annotated in Delphinium, 111 unique genes are identified across all samples. Genome structure, distributions of simple sequence repeats and short dispersed repeats, as well as gene content are similar among these Delphinium taxa. Nine hypervariable intergenic spacers and protein coding regions, including ndhF-trnL(TAG), rpl16-intron, rpl33, rps15, rps18, trnK(TTT)-trnQ(TTG), trnP(TGG)-psaJ, trnT(GGT)-psbD and ycf1, are identified among 13 perennial Delphinium. Selective pressure and codon usage bias of all the plastid genes are performed within 14 Delphinium taxa. Phylogenetic analysis based on 14 Delphinium plastomes, alongside two Aconitum (Ranunculaceae) species serving as outgroup taxa, reveals the monophyletic nature of Delphinium. Our findings further discern Delphinium into two distinct clades: perennial species (clade I) and annual species (clade II). In addition, compared with the nrDNA ITS topology, cytological data and morphological characters, D. mollifolium and D. maackianum showed potential involvement in hybridization or polyploidization processes. Excluding these two species, the perennial Delphinium (clade I) exhibits a stronger consistency with the morphology-based system that utilized seed morphology. CONCLUSION: This study represents the first comprehensive analysis of plastomic variations among Delphinium taxa, based on the examination of 14 complete plastomes. The chloroplast genome structure of Delphinium is similar to other angiosperms and possesses the typical quadripartite structure with the conserved genome arrangement and gene features. In addition, the variation of non-coding regions is larger than coding regions of the chloroplast genome. Through DNA sequence divergence across Delphinium plastomes and subsequent phylogenomic analyses ndhF-trnL(TAG) and ycf1 are identified as promising molecular markers. These highly variable loci held significant potential for future phylogenetic and phylogeographic studies on Delphinium. Our phylogenomic analyses based on the whole plastomes, concatenation of 132 unique intergenic spacer regions, concatenation of 77 unique protein-coding genes and nrDNA ITS, all support the monophyly of Delphinium and perennial taxa clusters together into one clade within this genus. These findings provide crucial data for systematic, phylogenomic and evolutionary research in the genus for future studies.

Comparative Analysis of Complete Chloroplast Genomes of Rubus in China: Hypervariable Regions and Phylogenetic Relationships.

With more than 200 species of native Rubus, China is considered a center of diversity for this genus. Due to a paucity of molecular markers, the phylogenetic relationships for this genus are poorly understood. In this study, we sequenced and assembled the plastomes of 22 out of 204 Chinese Rubus species (including varieties) from three of the eight sections reported in China, i.e., the sections Chamaebatus, Idaeobatus, and Malachobatus. Plastomes were annotated and comparatively analyzed with the inclusion of two published plastomes. The plastomes of all 24 Rubus species were composed of a large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs), and ranged in length from 155,464 to 156,506 bp. We identified 112 unique genes, including 79 protein-coding genes, 29 transfer RNAs, and four ribosomal RNAs. With highly consistent gene order, these Rubus plastomes showed strong collinearity, and no significant changes in IR boundaries were noted. Nine divergent hotspots were identified based on nucleotide polymorphism analysis: trnH-psbA, trnK-rps16, rps16-trnQ-psbK, petN-psbM, trnT-trnL, petA-psbJ, rpl16 intron, ndhF-trnL, and ycf1. Based on whole plastome sequences, we obtained a clearer phylogenetic understanding of these Rubus species. All sampled Rubus species formed a monophyletic group; however, sections Idaeobatus and Malachobatus were polyphyletic. These data and analyses demonstrate the phylogenetic utility of plastomes for systematic research within Rubus.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data.

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.

Chloroplast Genomes of Vitis flexuosa and Vitis amurensis: Molecular Structure, Phylogenetic, and Comparative Analyses for Wild Plant Conservation.

The chloroplast genome plays a crucial role in elucidating genetic diversity and phylogenetic relationships. Vitis vinifera L. (grapevine) is an economically important species, prompting exploration of wild genetic resources to enhance stress resilience. We meticulously assembled the chloroplast genomes of two Korean Vitis L. species, V. flexuosa Thunb. and V. amurensis Rupr., contributing valuable data to the Korea Crop Wild Relatives inventory. Through exhaustive specimen collection spanning diverse ecological niches across South Korea, we ensured comprehensive representation of genetic diversity. Our analysis, which included rigorous codon usage bias assessment and repeat analysis, provides valuable insights into amino acid preferences and facilitates the identification of potential molecular markers. The assembled chloroplast genomes were subjected to meticulous annotation, revealing divergence hotspots enriched with nucleotide diversity, thereby presenting promising candidates for DNA barcodes. Additionally, phylogenetic analysis reaffirmed intra-genus relationships and identified related crops, shedding light on evolutionary patterns within the genus. Comparative examination with chloroplast genomes of other crops uncovered conserved sequences and variable regions, offering critical insights into genetic evolution and adaptation. Our study advances the understanding of chloroplast genomes, genetic diversity, and phylogenetic relationships within Vitis species, thereby laying a foundation for enhancing grapevine genetic diversity and resilience to environmental challenges.

Complete Organelle Genome of the Desiccation-Tolerant (DT) Moss Tortula atrovirens and Comparative Analysis of the Pottiaceae Family.

Tortula atrovirens (Sm.) Lindb. is an important component of biological soil crusts and possesses an extraordinary tolerance against desiccation in dryland habitats. However, knowledge of the organelle genome of this desiccation-tolerant (DT) moss is still lacking. Here, we assembled the first reported Tortula organelle genome and conducted a comprehensive analysis within the Pottiaceae family. T. atrovirens exhibited the second largest chloroplast genome (129,646 bp) within the Pottiaceae, whereas its mitogenome (105,877 bp) and those of other mosses were smaller in size compared to other land plants. The chloroplast and mitochondrial genomes of T. atrovirens were characterized by the expansion of IR boundaries and the absence of homologous recombination-mediated by large repeats. A total of 57 RNA editing sites were detected through mapping RNA-seq data. Moreover, the gene content and order were highly conserved among the Pottiaceae organelle genomes. Phylogenetic analysis showed that bryophytes are paraphyletic, with their three lineages (hornworts, mosses, and liverworts) and vascular plants forming successive sister clades. Timmiella anomala is clearly separated from the monophyletic Pottiaceae, and T. atrovirens is closely related to Syntrichia filaris within the Pottioideae. In addition, we detected four hypervariable regions for candidate-molecular markers. Our findings provide valuable insights into the organelle genomes of T. atrovirens and the evolutionary relationships within the Pottiaceae family, facilitating future discovery of DT genetic resources from bryophytes.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Four Bromus (Poaceae, Bromeae) Species.

Bromus (Poaceae Bromeae) is a forage grass with high adaptability and ecological and economic value. Here, we sequenced Bromus ciliatus, Bromus benekenii, Bromus riparius, and Bromus rubens chloroplast genomes and compared them with four previously described species. The genome sizes of Bromus species ranged from 136,934 bp (Bromus vulgaris) to 137,189 bp (Bromus ciliates, Bromus biebersteinii), with a typical quadripartite structure. The studied species had 129 genes, consisting of 83 protein-coding, 38 tRNA-coding, and 8 rRNA-coding genes. The highest GC content was found in the inverted repeat (IR) region (43.85-44.15%), followed by the large single-copy (LSC) region (36.25-36.65%) and the small single-copy (SSC) region (32.21-32.46%). There were 33 high-frequency codons, with those ending in A/U accounting for 90.91%. A total of 350 simple sequence repeats (SSRs) were identified, with single-nucleotide repeats being the most common (61.43%). A total of 228 forward and 141 palindromic repeats were identified. No reverse or complementary repeats were detected. The sequence identities of all sequences were very similar, especially with respect to the protein-coding and inverted repeat regions. Seven highly variable regions were detected, which could be used for molecular marker development. The constructed phylogenetic tree indicates that Bromus is a monophyletic taxon closely related to Triticum. This comparative analysis of the chloroplast genome of Bromus provides a scientific basis for species identification and phylogenetic studies.

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora.

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

The chloroplast genomes of two medicinal species (Veronica anagallis-aquatica L. and Veronica undulata Wall.) and its comparative analysis with related Veronica species.

Veronica anagallis-aquatica L. and Veronica undulata Wall. are widely used ethnomedicinal plants in China. The two species have different clinical efficacies, while their extremely similar morphology and unclear interspecific relationship make it difficult to accurately identify them, leading to increased instances of mixed usage. This article reports on the complete chloroplast genomes sequence of these two species and their related Veronica species to conduct a comparative genomics analysis and phylogenetic construction. The results showed that the chloroplast (cp) genomes of Veronica exhibited typical circular quadripartite structures, with total lengths of 149,386 to 152,319 base pairs (bp), and GC content of 37.9 to 38.1%, and the number of genes was between 129-134. The total number of simple sequence repeats (SSRs) in V. anagallis-aquatica and V. undulata is 37 and 36, while V. arvensis had the highest total number of 56, predominantly characterized by A/T single bases. The vast majority of long repeat sequence types are forward repeats and palindromic repeats. Selective Ka/Ks values showed that three genes were under positive selection. Sequence differences often occur in the non-coding regions of the large single-copy region (LSC) and small single-copy region (SSC), with the lowest sequence variation in the inverted repeat regions (IR). Seven highly variable regions (trnT-GGU-psbD, rps8-rpl16, trnQ-UUG, trnN-GUU-ndhF, petL, ycf3, and ycf1) were detected, which may be potential molecular markers for identifying V. anagallis-aquatica and V. undulata. The phylogenetic tree indicates that there is a close genetic relationship between the genera Veronica and Neopicrorhiza, and V. anagallis-aquatica and V. undulata are sister groups. The molecular clock analysis results indicate that the divergence time of Veronica may occur at ∼ 9.09 Ma, and the divergence time of these two species occurs at ∼ 0.48 Ma. It is speculated that climate change may be the cause of Veronica species diversity and promote the radiation of the genus. The chloroplast genome data of nine Veronica specie provides important insights into the characteristics and evolution of the chloroplast genome of this genus, as well as the phylogenetic relationships of the genus Veronica.

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis.

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

Complete chloroplast genomes and phylogenetic relationships of Pedicularis chinensis and Pedicularis kansuensis.

The complete cp genomes of Pedicularis chinensis (GenBank accession number: OQ587614) and Pedicularis kansuensis (GenBank accession number: OQ587613) were sequenced, assembled, and annotated. Their chloroplast (cp) genome lengths were 146,452 bp, and 146,852 bp, respectively; 120 and 116 genes were identified, comprising 75 and 72 protein-coding genes (PCGs), 37 and 36 transfer RNA (tRNA) genes, and 8 and 8 ribosomal RNA (rRNA) genes, for P. chinensis and P. kansuensis, respectively. A simple sequence repeat (SSR) analysis revealed that the repetitive sequences were mainly composed of mononucleotide repeats (A/T motif) and dinucleotide repeats (AT/TA motif). Comparative genomics identified several variant genes (rpl22, rps19, rpl12, ycf1, trnH, psbA, and ndhH) and variant regions (trnS-GGA, trnV-UAC, ndhJ-trnV, ycf4-cemA, ndhE-nhdG, and rpl32-trnL) with a high Pi, indicating the potential to serve as deoxyribo nucleic acid (DNA) barcodes for Pedicularis species identification. The results show that the cp genomes of P. chinensis and P. kansuensis were the same as those of other plants in Pedicularis, with different degrees of AT preference for codons. Large differences in the number of SSRs and the expansion of the inverted repeat (IR) region showed strong variability and interspecific differentiation between these two species and other species represented in the genus Pedicularis. A phylogenetic analysis showed that P. kansuensis had the closest relationship with P. oliveriana, and P. chinensis had the closest relationship with P. aschistorhyncha. These results will facilitate the study of the phylogenetic classification and interspecific evolution of Pedicularis plants.

Characterization of the DNA accessibility of chloroplast genomes in grasses.

Although the chloroplast genome (cpDNA) of higher plants is known to exist as a large protein-DNA complex called 'plastid nucleoid', researches on its DNA state and regulatory elements are limited. In this study, we performed the assay for transposase-accessible chromatin sequencing (ATAC-seq) on five common tissues across five grasses, and found that the accessibility of different regions in cpDNA varied widely, with the transcribed regions being highly accessible and accessibility patterns around gene start and end sites varying depending on the level of gene expression. Further analysis identified a total of 3970 putative protein binding footprints on cpDNAs of five grasses. These footprints were enriched in intergenic regions and co-localized with known functional elements. Footprints and their flanking accessibility varied dynamically among tissues. Cross-species analysis showed that footprints in coding regions tended to overlap non-degenerate sites and contain a high proportion of highly conserved sites, indicating that they are subject to evolutionary constraints. Taken together, our results suggest that the accessibility of cpDNA has biological implications and provide new insights into the transcriptional regulation of chloroplasts.

Comparative analysis the chloroplast genomes of Celastrus (Celastraceae) species: Provide insights into molecular evolution, species identification and phylogenetic relationships.

BACKGROUND: The genus Celastrus is an important medicinal plant resource. The similarity of morphology and the lack of complete chloroplast genome analysis have significantly impeded the exploration of species identification, molecular evolution and phylogeny of Celastrus. PURPOSE: In order to resolve the phylogenic controversy of Celastrus species, the chloroplast genome comparative analysis was performed to provide genetic evidence. METHODS: In this study, we collected and sequenced ten chloroplast genomes of Celastrus species from China and downloaded three chloroplast genomes from the databases. The chloroplast genomes were compared and analyzed to explore their characteristics and evolution. Furthermore, the phylogenetic relationships of Celastrus species were inferred based on the whole chloroplast genomes and protein-coding genes. RESULTS: All the 13 Celastrus species chloroplast genomes showed a typical quadripartite structure with genome sizes ranging from 155,113 to 157,366 bp. The intron loss of the rps16 gene occurred in all the 13 Celastrus species. The GC content, gene sequence, repeat types and codon bias pattern were highly conserved. Ten highly variation regions were identified, which can be used as potential DNA markers in molecular identification of Celastrus species. Eight genes, including accD, atp4, ndhB, rpoC1, rbcL, rpl2, rpl20 and ycf1, were detected to experience positive selection. Phylogenetic analysis showed that Celastrus was a monophyletic group and Tripterygium was the closest sister-group. Noteworthy, C. gemmatus Loes. and C. orbiculatus Thunb. can be discriminated using the chloroplast genome as a super barcode. The comparative and phylogenetic analysis results proposed that C. tonkinensis Pitard. was the synonym of C. hindsii Benth. CONCLUSION: The comparative analysis of the Celastrus chloroplast genomes can provide comprehensive genetic evidence for molecular evolution, species identification and phylogenetic relationships.