Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Plastome comparison and phylogenomics of Chinese endemic Schnabelia (Lamiaceae): insights into plastome evolution and species divergence

BACKGROUND

Schnabelia species, herbaceous perennial plants within the Lamiaceae family, possess medicinal value and are endemic to China. While previous studies have focused on morphological classification, molecular systematics, and medicinal components, there has been limited research on phylogenomics. To reveal their plastid genome characteristics and phylogenetic relationships, we sequenced and assembled the plastomes of all five Schnabelia species (S. oligophylla, S. tetrodonta, S. nepetifolia, S. terniflora, S. aureoglandulosa), conducted comparative genomic analyses, and constructed a phylogenetic tree incorporating closely related taxa in subfamily Ajugoideae, as well as conducting divergence time estimation.

RESULTS

Plastome size of the five species ranged from 155,733 bp to 156,944 bp, encompassing 115 unique genes, with a GC content of 37.8% same for all species. Five intergenic spacer regions (trnH-GUG-psbA, trnK-UUU-matK, petB-petD, ndhD-psaC, ndhA-ndhH) were identified as divergence hotspots. Gene selection pressure analysis demonstrated that all genes were under negative selection. Phylogenetic relationship of Ajugoideae species based on plastomes confirmed the monophyly of Schnabelia. Two clades within Schnabelia were supported, one containing two original species and the other comprising three species transferred from Caryopteris. The stem age of the Schnabelia is estimated to be approximately 30.24 Ma, with the split of two Sections occurring around 12.60 Ma.

CONCLUSIONS

We revealed plastid genome evolutionary features for five species within the genus Schnabelia. The identified highly variable regions can provide a tool for future identification of these medicinal plants. The diversification of Schnabelia during middle Miocene and the Quaternary suggests that historical geological and climatic shifts facilitated species differentiation. These findings enhance our understanding of Schnabelia's evolution and support future research on chloroplast diversity, aiding conservation and sustainable use.

Gibberellic acid and light effects on seed germination in the seagrass Zostera marina

Seagrass meadows have been heavily affected by human activities, with Zostera marina L. (Zosteraceae) being one of the most impacted species. Seed-based methods are currently the preferred approach for their restoration, yet low germination rates and poor seedling establishment remain significant challenges. This study explored the combined effects of light spectra (white, red, and darkness), photoperiod, and gibberellic acid (GA3-0, 50, 500, and 1000 mg L-1) on Z. marina seed germination using a fully crossed incubation experiment. Penalised logistic regression and Cox proportional hazards analysis were chosen to account for low germination events and to analyse the temporal dynamics of germination. We found that light conditions, particularly red light and darkness, when combined with GA3, significantly enhanced germination probability. Furthermore, mid (50 mg L-1) and high (500 mg L-1) GA3 concentrations reduced time-to-germination. Morphometric analysis of the cotyledonary and leaf tissue development indicates no adverse effects of the treatments on seedling development. Our findings suggest that light and GA3 treatments effectively improve germination success and reduce dormancy in Z. marina seeds. Seed treatments can mitigate stress- or manipulation-induced dormancy and can represent a viable strategy for on-demand germination, such as in the context of seed-based restoration efforts.

The organelle genomes of the endangered seagrass Zostera caespitosa reveal sequence divergences, massive gene transfer, and uncommon RNA editing types

Introduction

Zostera caespitosa, a rare submerged angiosperm, is considered endemic to the northwestern Pacific.

Methods

This study assembled and compared the mitochondrial (mt) and chloroplast (cp) genomes of Z. caespitosa to understand the organelle evolutionary patterns.

Results and discussion

The cp genome (143,972 bp) was the second smallest within the seagrasses, whereas the mt genomes (192,246 bp) of Z. caespitosa and other seagrasses were smaller compared to those of other monocotyledons. The protein-coding genes (PCGs) in the organelle genome exhibit a strong A/U bias at codon endings, a selection-driven codon bias. The rates of nonsynonymous (Ka) and synonymous (Ks) substitutions in the mt genes of Zostera were two times higher than those in the cp genes. Additionally, 50 mitochondrial plastid DNA (MTPT) segments, totaling 44,662 bp, were identified, constituting 23.23% of the mt genome, which is significantly higher than those in most land plants. Phylogenetic analysis of 13 seagrass core cp-PCGs supported previous studies showing two genera in family Zosteraceae: Phyllospadix and Zostera, the latter comprising Zostera and Zosterella as subgenera. RNA editing was remarkably abundant in the 167 mt-PCGs and 172 in cp-PCGs, particularly in the cp genome. There are 11 different RNA editing types in the cp and 3 in the mt, most of which are C to U. Unexpectedly rare editing events, such as A to C, A to U, U to A, G to C, and U to G, have also been found in the cp.

The first complete chloroplast genome of Halodule uninervis (Forssk.) Boiss. 1882 (Cymodoceaceae)

Halodule uninervis (Forssk.) Boiss. 1882 is a seagrass species. The complete chloroplast genome is 164,005 bp in length, with a GC content of 36.3%. This genome consists of a large single-copy region of 87,522 bp, a small single-copy region of 11,481 bp and a pair of inverse repeat regions of 32,501 bp each. A total of 131 genes were annotated in the genome, including 87 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. Phylogenetic analysis revealed that H. uninervis shares the closest taxonomic relationship with Cymodocea rotundata and Syringodium isoetifolium, and all three species belong to the family Cymodoceaceae.

Complete chloroplast genome sequence of Artemisia littoricola (Asteraceae) from Dokdo Island Korea: genome structure, phylogenetic analysis, and biogeography study

The Asteraceae family, particularly the Artemisia genus, presents taxonomic challenges due to limited morphological characteristics and frequent natural hybridization. Molecular tools, such as chloroplast genome analysis, offer solutions for accurate species identification. In this study, we sequenced and annotated the chloroplast genome of Artemisia littoricola sourced from Dokdo Island, employing comparative analyses across six diverse Artemisia species. Our findings reveal conserved genome structures with variations in repeat sequences and junction boundaries. Notably, the chloroplast genome of A. littoricola spans 150,985 bp, consistent with other Artemisia species, and comprises 131 genes, including 86 protein-coding, 37 tRNA, and 8 rRNA genes. Among these genes, 16 possess a single intron, while clpP and ycf3 exhibit two introns each. Furthermore, 18 genes display duplicated copies within the IR regions. Moreover, the genome possesses 42 Simple Sequence Repeats (SSRs), predominantly abundant in A/T content and located within intergenic spacer regions. The analysis of codon usage revealed that the codons for leucine were the most frequent, with a preference for ending with A/U. While the chloroplast genome exhibited conservation overall, non-coding regions showed lower conservation compared to coding regions, with the Inverted Repeat (IR) region displaying higher conservation than single-copy regions. Phylogenetic analyses position A. littoricola within subgenus Dracunculus, indicating a close relationship with A. scoparia and A. desertorum. Additionally, biogeographic reconstructions suggest ancestral origins in East Asia, emphasizing Mongolia, China (North East and North Central and South Central China), and Korea. This study underscores the importance of chloroplast genomics in understanding Artemisia diversity and evolution, offering valuable insights into taxonomy, evolutionary patterns, and biogeographic history. These findings not only enhance our understanding of Artemisia's intricate biology but also contribute to conservation efforts and facilitate the development of molecular markers for further research and applications in medicine and agriculture.

Combined transcriptomics and metabolomics analysis reveals salinity stress specific signaling and tolerance responses in the seagrass Zostera japonica

KEY MESSAGE

Multiple regulatory pathways of Zostera japonica to salt stress were identified through growth, physiological, transcriptomic and metabolomic analyses. Seagrasses are marine higher submerged plants that evolved from terrestrial monocotyledons and have fully adapted to the high saline seawater environment during the long evolutionary process. As one of the seagrasses growing in the intertidal zone, Zostera japonica not only has the ability to quickly adapt to short-term salt stress but can also survive at salinities ranging from the lower salinity of the Yellow River estuary to the higher salinity of the bay, making it a good natural model for studying the mechanism underlying the adaptation of plants to salt stress. In this work, we screened the growth, physiological, metabolomic, and transcriptomic changes of Z. japonica after a 5-day exposure to different salinities. We found that high salinity treatment impeded the growth of Z. japonica, hindered its photosynthesis, and elicited oxidative damage, while Z. japonica increased antioxidant enzyme activity. At the transcriptomic level, hypersaline stress greatly reduced the expression levels of photosynthesis-related genes while increasing the expression of genes associated with flavonoid biosynthesis. Meanwhile, the expression of candidate genes involved in ion transport and cell wall remodeling was dramatically changed under hypersaline stress. Moreover, transcription factors signaling pathways such as mitogen-activated protein kinase (MAPK) were also significantly influenced by salt stress. At the metabolomic level, Z. japonica displayed an accumulation of osmolytes and TCA mediators under hypersaline stress. In conclusion, our results revealed a complex regulatory mechanism in Z. japonica under salt stress, and the findings will provide important guidance for improving salt resistance in crops.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

The complete chloroplast genomes of three Alismataceae species, including the medicinally important Alisma orientale

Alismataceae is one of the early diverged families of monocotyledonous plants. We report the complete chloroplast genomes of three Alisma species, including Alisma orientale (Sam.) Juzep. 1934, A. subcordatum Raf. 1908, and A. triviale Pursh 1813, of which A. orientale is a traditional Chinese medical plant used widely to treat diuretics, diabetes, hepatitis, and inflammation. We sequenced the complete chloroplast genomes with the Illumina Nova-Seq 6000 platform using herbarium collections. The chloroplast genomes of A. orientale, A. subcordatum and A. triviale are 159,861 bp, 160,180 bp, and 159,727 bp in length, respectively. The three chloroplast genomes each contain 113 genes, including four rRNAs, 30 tRNAs genes, and 79 protein-coding genes, and the average GC content is 36.0%. Based on the whole chloroplast genomes of 19 species of Alismataceae and the close allies, the medicinally important A. orientale was found to be closely related to another medicinal plant Alisma plantago-aquatica L. 1753 in the phylogenetic analysis. The genus Alisma was supported to be monophyletic.

Plastid genome and its phylogenetic implications of Asiatic Spiraea (Rosaceae)

BACKGROUND

Spiraea L. is a genus comprising approximately 90 species that are distributed throughout the northern temperate regions. China is recognized as the center of species diversity for this genus, hosting more than 70 species, including 47 endemic species. While Spiraea is well-known for its ornamental value, its taxonomic and phylogenetic studies have been insufficient.

RESULTS

In this study, we conducted sequencing and assembly of the plastid genomes (plastomes) of 34 Asiatic Spiraea accessions (representing 27 Asiatic Spiraea species) from China and neighboring regions. The Spiraea plastid genome exhibits typical quadripartite structures and encodes 113-114 genes, including 78-79 protein-coding genes (PCGs), 30 tRNA genes, and 4 rRNA genes. Linear regression analysis revealed a significant correlation between genome size and the length of the SC region. By the sliding windows method, we identified several hypervariable hotspots within the Spiraea plastome, all of which were localized in the SC regions. Our phylogenomic analysis successfully established a robust phylogenetic framework for Spiraea, but it did not support the current defined section boundaries. Additionally, we discovered that the genus underwent diversification after the Early Oligocene (~ 30 Ma), followed by a rapid speciation process during the Pliocene and Pleistocene periods.

CONCLUSIONS

The plastomes of Spiraea provided us invaluable insights into its phylogenetic relationships and evolutionary history. In conjunction with plastome data, further investigations utilizing other genomes, such as the nuclear genome, are urgently needed to enhance our understanding of the evolutionary history of this genus.

Comparative plastome genomics and phylogenetic relationships of the genus Trollius

Trollius, a genus in the Ranunculaceae family, has significant medicinal and ornamental value. It is widely distributed in China with 16 different species accepted. However, due to the lack of enough samples and information sites, the molecular phylogenetic relationships of Trollius have been unresolved till now. Here we sequenced, assembled and annotated the plastomes of 16 Trollius species to investigate their genomic characteristics, inverted repeat (IR) boundaries, sequence repeats, and hypervariable loci. In addition, the phylogenetic relationships of this genus was reconstructed based on the whole plastomes and the protein-coding sequences data-sets. The plastomes of Trollius ranged between 159,597 bp and 160,202 bp in length, and contained 113 unique genes, including 79 protein coding, 30 tRNA, and 4 rRNA. The IR boundaries were relatively conserved within the genus Trollius. 959 simple sequence repeats and 657 long sequence repeats were detected in the Trollius plastomes. We identified 12 highly polymorphic loci (Pi > 0.0115) that can be used as plastid markers in molecular identification and phylogenetic investigation of the genus. Besides, Trollius was a monophyletic group with the earliest divergence clade being Trollius lilacinus Bunge, and the remaining species were divided into two strongly-supported clades. The phylogeny in our study supported the traditional classification systems based on the color of sepal, but not the previous classification system based on the types and relative lengths of the nectaries, and distribution. The genomic resources provided in our study can be used in the taxonomy of the genus Trollius, promoting the development and utilization of this genus.

Comparative chloroplast genomics reveals the phylogeny and the adaptive evolution of Begonia in China

BACKGROUND

The Begonia species are common shade plants that are mostly found in southwest China. They have not been well studied despite their medicinal and decorative uses because gene penetration, decreased adaptability, and restricted availability are all caused by frequent interspecific hybridization.

RESULT

To understand the patterns of mutation in the chloroplast genomes of different species of Begonia, as well as their evolutionary relationships, we collected seven Begonia species in China and sequenced their chloroplast genomes. Begonia species exhibit a quadripartite structure of chloroplast genomes (157,634 - 169,694 bp), consisting of two pairs of inverted repeats (IR: 26,529 - 37,674 bp), a large single copy (LSC: 75,477 - 86,500 bp), and a small single copy (SSC: 17,861 - 18,367 bp). 128-143 genes (comprising 82-93 protein-coding genes, 8 ribosomal RNAs, and 36-43 transfer RNAs) are found in the chloroplast genomes. Based on comparative analyses, this taxon has a relatively similar genome structure. A total of six substantially divergent DNA regions (trnT-UGU-trnL-UAA, atpF-atpH, ycf4-cemA, psbC-trnS-UGA, rpl32-trnL-UAG, and ccsA-ndhD) are found in the seventeen chloroplast genomes. These regions are suitable for species identification and phylogeographic analysis. Phylogenetic analysis shows that Begonia species that were suited to comparable environments grouped in a small clade and that all Begonia species formed one big clade in the phylogenetic tree, supporting the genus' monophyly. In addition, positive selection sites were discovered in eight genes (rpoC1, rpoB, psbE, psbK, petA, rps12, rpl2, and rpl22), the majority of which are involved in protein production and photosynthesis.

CONCLUSION

Using these genome resources, we can resolve deep-level phylogenetic relationships between Begonia species and their families, leading to a better understanding of evolutionary processes. In addition to enhancing species identification and phylogenetic resolution, these results demonstrate the utility of complete chloroplast genomes in phylogenetically and taxonomically challenging plant groupings.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Chloroplast genomic comparison provides insights into the evolution of seagrasses

BACKGROUND

Seagrasses are a polyphyletic group of monocotyledonous angiosperms that have evolved to live entirely submerged in marine waters. Thus, these species are ideal for studying plant adaptation to marine environments. Herein, we sequenced the chloroplast (cp) genomes of two seagrass species (Zostera muelleri and Halophila ovalis) and performed a comparative analysis of them with 10 previously published seagrasses, resulting in various novel findings.

RESULTS

The cp genomes of the seagrasses ranged in size from 143,877 bp (Zostera marina) to 178,261 bp (Thalassia hemprichii), and also varied in size among different families in the following order: Hydrocharitaceae > Cymodoceaceae > Ruppiaceae > Zosteraceae. The length differences between families were mainly related to the expansion and contraction of the IR region. In addition, we screened out 2,751 simple sequence repeats and 1,757 long repeat sequence types in the cp genome sequences of the 12 seagrass species, ultimately finding seven hot spots in coding regions. Interestingly, we found nine genes with positive selection sites, including two ATP subunit genes (atpA and atpF), three ribosome subunit genes (rps4, rps7, and rpl20), one photosystem subunit gene (psbH), and the ycf2, accD, and rbcL genes. These gene regions may have played critical roles in the adaptation of seagrasses to diverse environments. In addition, phylogenetic analysis strongly supported the division of the 12 seagrass species into four previously recognized major clades. Finally, the divergence time of the seagrasses inferred from the cp genome sequences was generally consistent with previous studies.

CONCLUSIONS

In this study, we compared chloroplast genomes from 12 seagrass species, covering the main phylogenetic clades. Our findings will provide valuable genetic data for research into the taxonomy, phylogeny, and species evolution of seagrasses.

CPJSdraw: analysis and visualization of junction sites of chloroplast genomes

Background

Chloroplast genomes are usually circular molecules, and most of them are tetrad structures with two inverted repeat (IR) regions, a large single-copy region, and a small single-copy region. IR contraction and expansion are among the genetic diversities during the evolution of plant chloroplast genomes. The only previously released tool for the visualization of junction sites of the regions does not consider the diversity of the starting point of genomes, which leads to incorrect results or even no results for the examination of IR contraction and expansion.

Results

In this work, a new tool named CPJSdraw was developed for visualizing the junction sites of chloroplast genomes. CPJSdraw can format the starting point of the irregular linearized genome, correct the junction sites of IR and single-copy regions, display the tetrad structure, visualize the junction sites of any number (≥1) of chloroplast genomes, show the transcription direction of genes adjacent to junction sites, and indicate the IR expansion or contraction of chloroplast genomes.

Conclusions

CPJSdraw is a software that is universal and reliable in analysis and visualization of IR expansion or contraction of chloroplast genomes. CPJSdraw has more accurate analysis and more complete functions when compared with previously released tool. CPJSdraw as a perl package and tested data are available at http://dx.doi.org/10.5281/zenodo.7669480 for English users. In addition, an online version with a Chinese interface is available at http://cloud.genepioneer.com:9929/#/tool/alltool/detail/335.

Comparative analysis of mitochondrial genomes reveals marine adaptation in seagrasses

BACKGROUND

Seagrasses are higher marine flowering plants that evolved from terrestrial plants, but returned to the sea during the early evolution of monocotyledons through several separate lineages. Thus, they become a good model for studying the adaptation of plants to the marine environment. Sequencing of the mitochondrial (mt) genome of seagrasses is essential for understanding their evolutionary characteristics.

RESULTS

In this study, we sequenced the mt genome of two endangered seagrasses (Zostera japonica and Phyllospadix iwatensis). These data and data on previously sequenced mt genomes from monocotyledons provide new evolutionary evidence of genome size reduction, gene loss, and adaptive evolution in seagrasses. The mt genomes of Z. japonica and P. iwatensis are circular. The sizes of the three seagrasses (including Zostera marine) that have been sequenced to date are smaller than that of other monocotyledons. Additionally, we found a large number of repeat sequences in seagrasses. The most abundant long repeat sequences were 31-40 bp repeats. Our study also found that seagrass species lost extensive ribosomal protein genes during evolution. The rps7 gene and the rpl16 gene of P. iwatensis are exceptions to this trend. The phylogenetic analysis based on the mt genome strongly supports the previous results. Furthermore, we identified five positive selection genes (atp8, nad3, nad6, ccmFn, and matR) in seagrasses that may be associated with their adaptation to the marine environment.

CONCLUSIONS

In this study, we sequenced and annotated the mt genomes of Z. japonica and P. iwatensis and compared them with the genome of other monocotyledons. The results of this study will enhance our understanding of seagrass adaptation to the marine environment and can inform further investigations of the seagrass mt genome.

Comparative and phylogenetic analyses of the chloroplast genome reveal the taxonomy of the Morus genus

Mulberry (genus Morus) is an economically important woody plant with an altered ploidy level. The variable number of Morus species recognized by different studies indicates that the genus is in need of revision. In this study, the chloroplast (CP) genomes of 123 Morus varieties were de novo assembled and systematically analyzed. The 123 varieties represented six Morus species, namely, Morus alba, Morus nigra, Morus notabilis, Morus rubra, Morus celtidifolia, and Morus serrata. The Morus CP genome was found to be 158,969~159,548 bp in size with 125 genes, including 81 protein coding, 36 tRNA, and 8 rRNA genes. The 87 out of 123 mulberry accessions were assigned to 14 diverse groups with identical CP genome, which indicated that they are maternally inherited and share 14 common ancestors. Then 50 diverse CP genomes occurred in 123 mulberry accessions for further study. The CP genomes of the Morus genus with a quadripartite structure have two inverted repeat (IR) regions (25,654~25,702 bp) dividing the circular genome into a large single-copy (LSC) region (87,873~88,243 bp) and small single-copy (SSC) region (19,740~19,994 bp). Analysis of the phylogenetic tree constructed using the complete CP genome sequences of Morus revealed a monophyletic genus and that M. alba consisted of two clades, M. alba var. alba and M. alba var. multicaulis. The Japanese cultivated germplasms were derived from M. alba var. multicaulis. We propose that the Morus genus be classified into six species, M. nigra, M. notabilis, M. serrata, M. celtidifolia, M. rubra, and M. alba with two subspecies, M. alba var. alba and M. alba var. multicaulis. Our findings provide a valuable resource for the classification, domestication, and breeding improvement of mulberry.

Comparative Analyses of Plastomes of Four Anubias (Araceae) Taxa, Tropical Aquatic Plants Endemic to Africa

Anubias Schott (Araceae) have high ornamental properties as aquarium plants. However, the genus has difficulties in species identification, and the mechanism of its adaptation to the aquatic environment is unknown. To better identify species and understand the evolutionary history of Anubias, the plastomes of Anubias barteri Schott, A. barteri var. nana (Engl.) Crusio, and A. hastifolia Engl., were sequenced. The sizes of the plastomes of Anubias ranged from 169,841 bp to 170,037 bp. These plastomes were composed of conserved quadripartite circular structures and comprised 112 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. The comparative analysis of genome structure, repeat sequences, codon usage and RNA editing sites revealed high similarities among the Anubias plastomes, indicating the conservation of plastomes of Anubias. Three spacer regions with relatively high nucleotide diversity, trnL-CAA-ndhB, ycf1-ndhF, and rps15-ycf1, were found within the plastomes of Anubias. Phylogenetic analysis, based on 75 protein-coding genes, showed that Anubias was sister to Montrichardia arborescens (L.) Schott (BS = 99). In addition, four genes (ccsA, matK, ndhF, and ycf4) that contain sites undergoing positive selection were identified within the Anubias plastomes. These genes may play an important role in the adaptation of Anubias to the aquatic environment. The present study provides a valuable resource for further studies on species identification and the evolutionary history of Anubias.

Sequence Characteristics and Phylogenetic Analysis of the Artemisia argyi Chloroplast Genome

Artemisia argyi Levl. et Van is an important Asteraceae species with a high medicinal value. There are abundant A. argyi germplasm resources in Asia, especially in China, but the evolutionary relationships of these varieties and the systematic localization of A. argyi in the family Asteraceae are still unclear. In this study, the chloroplast (cp) genomes of 72 A. argyi varieties were systematically analyzed. The 72 varieties originated from 47 regions in China at different longitudes, latitudes and altitudes, and included both wild and cultivated varieties. The A. argyi cp genome was found to be ∼151 kb in size and to contain 114 genes, including 82 protein-coding, 28 tRNA, and 4 rRNA genes. The number of short sequence repeats (SSRs) in A. argyi cp genomes ranged from 35 to 42, and most of them were mononucleotide A/T repeats. A total of 196 polymorphic sites were detected in the cp genomes of the 72 varieties. Phylogenetic analysis demonstrated that the genetic relationship between A. argyi varieties had a weak relationship with their geographical distribution. Furthermore, inverted repeat (IR) boundaries of 10 Artemisia species were found to be significantly different. A sequence divergence analysis of Asteraceae cp genomes showed that the variable regions were mostly located in single-copy (SC) regions and that the coding regions were more conserved than the non-coding regions. A phylogenetic tree was constructed using 43 protein-coding genes common to 67 Asteraceae species. The resulting tree was consistent with the traditional classification system; Artemisia species were clustered into one group, and A. argyi was shown to be closely related to Artemisia lactiflora and Artemisia montana. In summary, this study systematically analyzed the cp genome characteristics of A. argyi and compared cp genomes of Asteraceae species. The results provide valuable information for the definitive identification of A. argyi varieties and for the understanding of the evolutionary relationships between Asteraceae species.

Genome skimming approach reveals the gene arrangements in the chloroplast genomes of the highly endangered Crocus L. species: Crocus istanbulensis (B.Mathew) Rukšāns

Crocus istanbulensis (B.Mathew) Rukšāns is one of the most endangered Crocus species in the world and has an extremely limited distribution range in Istanbul. Our recent field work indicates that no more than one hundred individuals remain in the wild. In the present study, we used genome skimming to determine the complete chloroplast (cp) genome sequences of six C. istanbulensis individuals collected from the locus classicus. The cp genome of C. istanbulensis has 151,199 base pairs (bp), with a large single-copy (LSC) (81,197 bp), small single copy (SSC) (17,524 bp) and two inverted repeat (IR) regions of 26,236 bp each. The cp genome contains 132 genes, of which 86 are protein-coding (PCGs), 8 are rRNA and 38 are tRNA genes. Most of the repeats are found in intergenic spacers of Crocus species. Mononucleotide repeats were most abundant, accounting for over 80% of total repeats. The cp genome contained four palindrome repeats and one forward repeat. Comparative analyses among other Iridaceae species identified one inversion in the terminal positions of LSC region and three different gene (psbA, rps3 and rpl22) arrangements in C. istanbulensis that were not reported previously. To measure selective pressure in the exons of chloroplast coding sequences, we performed a sequence analysis of plastome-encoded genes. A total of seven genes (accD, rpoC2, psbK, rps12, ccsA, clpP and ycf2) were detected under positive selection in the cp genome. Alignment-free sequence comparison showed an extremely low sequence diversity across naturally occurring C. istanbulensis specimens. All six sequenced individuals shared the same cp haplotype. In summary, this study will aid further research on the molecular evolution and development of ex situ conservation strategies of C. istanbulensis.

Complete Chloroplast Genomes Provide Insights Into Evolution and Phylogeny of Campylotropis (Fabaceae)

The genus Campylotropis Bunge (Desmodieae, Papilionoideae) comprises about 37 species distributed in temperate and tropical Asia. Despite the great potential in soil conservation, horticulture, and medicine usage, little is known about the evolutionary history and phylogenetic relationships of Campylotropis due to insufficient genetic resources. Here, we sequenced and assembled 21 complete chloroplast genomes of Campylotropis species. In combination with the previously published chloroplast genomes of C. macrocarpa and closely related species, we conducted comparative genomics and phylogenomic analysis on these data. Comparative analysis of the genome size, structure, expansion and contraction of inverted repeat (IR) boundaries, number of genes, GC content, and pattern of simple sequence repeats (SSRs) revealed high similarities among the Campylotropis chloroplast genomes. The activities of long sequence repeats contributed to the variation in genome size and gene content in Campylotropis chloroplast genomes. The Campylotropis chloroplast genomes showed moderate sequence variation, and 13 highly variable regions were identified for species identification and further phylogenetic studies. We also reported one more case of matK pseudogene in the legume family. The phylogenetic analysis confirmed the monophyly of Campylotropis and the sister relationship between Lespedeza and Kummerowia, the latter two genera were then sister to Campylotropis. The intrageneric relationships of Campylotropis based on genomic scale data were firstly reported in this study. The two positively selected genes (atpF and rps19) and eight fast-evolving genes identified in this study may help us to understand the adaptation of Campylotropis species. Overall, this study enhances our understanding of the chloroplast genome evolution and phylogenetic relationships of Campylotropis.

Comparative chloroplast genome analyses of Paraboea (Gesneriaceae): Insights into adaptive evolution and phylogenetic analysis

Paraboea (Gesneriaceae) distributed in the karst areas of South and Southwest China and Southeast Asia, is an ideal genus to study the phylogeny and adaptive evolution of karst plants. In this study, the complete chloroplast genomes of twelve Paraboea species were sequenced and analyzed. Twelve chloroplast genomes ranged in size from 153166 to 154245 bp. Each chloroplast genome had a typical quartile structure, and relatively conserved type and number of gene components, including 131 genes which are composed of 87 protein coding genes, 36 transfer RNAs and 8 ribosomal RNAs. A total of 600 simple sequence repeats and 389 non-overlapped sequence repeats were obtained from the twelve Paraboea chloroplast genomes. We found ten divergent regions (trnH-GUG-psbA, trnM-CAU, trnC-GCA, atpF-atpH, ycf1, trnK-UUU-rps16, rps15, petL, trnS-GCU-trnR-UCU and psaJ-rpl33) among the 12 Paraboea species to be potential molecular markers. In the phylogenetic tree of 31 Gesneriaceae plants including twelve Paraboea species, all Paraboea species clustered in a clade and confirmed the monophyly of Paraboea. Nine genes with positive selection sites were detected, most of which were related to photosynthesis and protein synthesis, and might played crucial roles in the adaptability of Paraboea to diverse karst environments. These findings are valuable for further study of the phylogeny and karst adaptability of Gesneriaceae plants.