Intraspecific Differentiation of Styrax japonicus (Styracaceae) as Revealed by Comparative Chloroplast and Evolutionary Analyses

Styrax japonicus is a medicinal and ornamental shrub belonging to the Styracaceae family. To explore the diversity and characteristics of the chloroplast genome of S. japonicus, we conducted sequencing and comparison of the chloroplast genomes of four naturally distributed S. japonicus. The results demonstrated that the four chloroplast genomes (157,914–157,962 bp) exhibited a typical quadripartite structure consisting of a large single copy (LSC) region, a small single copy (SSC) region, and a pair of reverse repeats (IRa and IRb), and the structure was highly conserved. DNA polymorphism analysis revealed that three coding genes (infA, psbK, and rpl33) and five intergene regions (petA-psbJ, trnC-petN, trnD-trnY, trnE-trnT, and trnY-trnE) were identified as mutation hotspots. These genetic fragments have the potential to be utilized as DNA barcodes for future identification purposes. When comparing the boundary genes, a small contraction was observed in the IR region of four S. japonicus. Selection pressure analysis indicated positive selection for ycf1 and ndhD. These findings collectively suggest the adaptive evolution of S. japonicus. The phylogenetic structure revealed conflicting relationships among several S. japonicus, indicating divergent evolutionary paths within this species. Our study concludes by uncovering the genetic traits of the chloroplast genome in the differentiation of S. japonicus variety, offering fresh perspectives on the evolutionary lineage of this species.


Introduction
Chloroplasts in higher plants play a crucial role in the biosynthesis of amino acids, fatty acids, vitamins, and pigments [1], as well as being actively involved in photosynthesis [2].The chloroplast genome structure is characterized by a circular quadripartite organization consisting of two inverted repeat regions (IRs), a small single copy (SSC) region, and a large single copy (LSC) region that are highly conserved [3].Furthermore, chloroplast genes demonstrate a relatively high level of conservation and play a crucial role in fundamental biological processes such as photosynthesis, transcription, and translation [4].The size of the chloroplast genome varies among different species, typically ranging from 120 to 160 kb, and contains two sets of four ribosomal RNA genes and 30 tRNA genes [5,6].The majority of angiosperm species rely on matrilineal inheritance to ensure the stability of species evolution.Simultaneously, genetic mutation events provide valuable information for evolutionary studies, population classification, and species identification [7][8][9][10].Consequently, chloroplast genes present an ideal research subject for investigating species evolution [11].
The diversity of chloroplast genomes offers valuable resources for uncovering phylogenetic relationships across various dimensions [12][13][14].Furthermore, due to the maternal inheritance of chloroplasts in most angiosperms [15], variations in their association with nuclear phylogeny could offer valuable insights into the processes of speciation [16,17].The modification of gene content within the chloroplast genome may contribute to the adaptation of species to specific habitats and living strategies [18,19].Changes in the environment and habitat can impose selective genetic pressure, influencing the genes involved.Genes undergoing positive selection typically confer advantages in terms of individual fitness and reproductive capability [20].Therefore, the investigation of gene selection pressure and adaptive evolution has emerged as a prominent area of focus in molecular research, providing a solid groundwork for the study of germplasm resources.
S. japonicus Siebold & Zucc is a member of the Styracaceae family, which is a lowbranched landscape tree with aromatic flowers possessing significant ornamental and medicinal properties [21][22][23].S. japonicus has been utilized in the treatment of various conditions, including cough, bronchitis, and laryngitis [24,25] and is also employed in traditional Chinese herbal medicine for relieving toothaches and sore throats [26].Moreover, some scholars suggest that the distinctive properties of its nectar make it a potential additive in honey processing [27].It exhibits wide distribution across eastern Asia, including Korea, Japan, and southern China [28].Due to its wide distribution and morphological characteristics, there has been some confusion in the taxonomy of S. japonicus, making it one of the most heterogeneous species in Styrax [29][30][31].However, there is limited research on S. japonicus variety; Ju et al. [32] conducted an analysis of the mechanism of calyx color variation in two S. japonicus varieties using integrated metabolomics and transcriptomics.The extent of genomic diversity in chloroplasts remains entirely unexplored.
Although the nuclear genome predominantly dictates the morphological and physiological traits of plants, the chloroplast genome indirectly influences these attributes to a notable extent.For instance, Lee and colleagues discovered variations in the chloroplast genome of Cynanchum wilfordii, proposing a potential link between these genetic differences and the species' calyx variations [33].Similarly, Tian et al., in their investigation of Lindera obtusiloba, highlighted variations in chloroplast genome size and intraspecific diversity, which might be connected to the plant's adaptive responses to environmental pressures.Specifically, genes under positive selection could drive morphological adaptations in plants facing environmental stress [34].Sarwar et al. conducted a comprehensive analysis of Olea europaea cultivars, revealing significant disparities in plant height, blade size and shape, phyllotaxy, venation, and chlorophyll content among varieties from different origins.These variations might originate from the diversity within the chloroplast genome [35].In essence, the chloroplast genome's diversity significantly shapes plant morphology, with these morphological traits not only mirroring genetic variability but also potentially indicating adaptability to specific environmental conditions.Moreover, the interplay between the chloroplast and nuclear genomes could be crucial in plant speciation and evolutionary processes.Understanding these interactions offers valuable insights into the complex genetic architecture underlying plant diversity and adaptation.
To explore the diversity and evolution of the chloroplast genome of S. japonicus, we collected three samples with typical leaf variation and one wild-type sample with standard morphology from the natural population of S. japonicus in Shandong, China, and conducted chloroplast genome sequencing.The objective is to analyze the evolutionary mechanisms, intraspecific differentiation, and chloroplast genome diversity of S. japonicus from the perspective of the chloroplast genome.Specifically, (1) A1 and Table A1).Each of the three chloroplasts contains a total of 133 genes, including 87 protein-coding genes, 37 transfer RNA (tRNA) genes, 8 ribosomal RNA (rRNA) genes, and a pseudogene (ycf1).All genes can be categorized into three groups: (1) photosynthesis, (2) self-replication, and (3) others (Appendix B, Table A2).

Comparative Analysis
Based on the previously published S. japonicus chloroplast genome from NCBI (Accession number: MT178456; MZ285743; NC047429), the results indicate that four samples exhibit consistency with other S. japonicus boundary genes.However, there were varying degrees of expansion or contraction in the IR region (Figure 1), although such variations were minimal, with only a few bp.The chloroplast genome sequences of three varieties were analyzed using the reference chloroplast genome of wild type in mVISTA (Appendix B Figure A2), revealing a significant level of similarity and conservation among the chloroplast genomes within the three varieties.Importantly, coding regions exhibited higher levels of conservation compared to non-coding regions, while the IR regions demonstrated lower variability than the LSC and SSC regions.Through the analysis of the CDS sequence and complete chloroplast genome nucleotide polymorphisms, it was found that infA, psbK, and rpl33 exhibited high Pi values (Figure 2A).Furthermore, the intergenic regions petA-psbJ, trnC-petN, trnD-trnY, trnE-trnT, and trnY-trnE also exhibited elevated Pi values (Figure 2B).

Repetitive and Codon Usage Bias Analysis
Through SSR analysis, the four chloroplast genomes were found to contain a total of 58-61 SSRs, with the majority (77.59-79.66%)located in the LSC region.The IR regions contained between 3.39 and 6.90% of SSR loci, while the SSC region included between 14.75 and 16.95% (Figure 3A).In contrast, only WT exhibited dinucleotide repeat units (Figure 3B).Furthermore, the majority of repeats in all four species consist of mononucleotide repeats (Figure 3C).
The examination of chloroplast genomes demonstrated that the GC and GC3s composition within the codons was consistently lower than 0.5, suggesting a predilection towards A/T bases and A/T-ending codons in three chloroplast genomes.

RSCU Analysis and Selection in Evolution
The relative synonymous codon usage (RSCU) values showed similarity across the four genomes (Figure 4A).A total of 33 codons exhibited an RSCU value greater than 1 (Figure 4B), with only two of these codons ending in G (AUG and UUG).Among the codons with an RSCU value less than 1, except for AUA, CUA, and UGA, which terminate in A, the remaining codons conclude with either C or G. Due to the majority of genes having Ks values of 0, this leads to an ineffective Ka/Ks ratio (Figure 4C).Out of the 79 common genes, only 8 had valid values, with ycf1 and ndhD exhibiting values greater than 1.

Phylogenetic Analysis
The best-fit model GTR+G determined by ModelFinder was utilized for maximum likelihood analysis.Based on the Styrax topology, the interbranches within the genus are notably reduced, suggesting a lower level of divergence between Styrax species (Figure 5A).Phylogenetic analysis of the complete chloroplast genome and CDS sequences revealed that S. japonicus formed three clades within Styrax.The four samples sequenced in this study grouped together, with YZ and WT forming one branch, and LY and JR forming another.Three additional reference sequences reveal that S. japonicus (MZ285743) forms a branch with S. calvescens and as a sister clade to two other S. japonicus in the complete chloroplast genome tree (Figure 5B).However, in the CDS tree, S. japonicus (MZ285743) appears as a single branch, while the remaining two S.japonicus (NC047429 and MT178456) consistently form a branch in both trees.In general, the phylogenetic analysis revealed that seven samples of S. japonicus formed three distinct clusters, indicating a high level of diversity in the chloroplast genome and potential phylogenetic conflict within this species.

Discussion
In this study, we performed sequencing of three S. japonicus varieties and conducted a comparative analysis with the wild type and other S. japonicus to explore the correlation between leaf morphological variation and evolution in the chloroplast genome dimension.The results revealed a high degree of similarity among their chloroplast genomes; however, sliding window analysis identified potential mutation hotspots that could serve as DNA barcodes for the identification and germplasm resource of S. japonicus taxa.Additionally, (1) the boundary genes detection revealed varying degrees of expansion and contraction in the IR region of the three varieties; (2) ndhD and ycf1 were found to be under positive selection; (3) phylogenetic analysis revealed the formation of three clades by seven S. japonicus.In the following, we will discuss and explore these three points.

IR Region Contraction within S. japonicus
The findings of this study demonstrate that the boundary genes in all S. japonicus are consistent; however, the IR regions of the three varieties exhibit varying degrees of contraction compared to different S. japonicus samples.Despite being minor, only a few bp, it is noteworthy that such changes occur within the same species.
The IR region of the chloroplast genome is generally considered to be highly conserved [36], but variations in its boundary genes can lead to changes in gene loss or the formation of pseudogenes [37].This phenomenon is common in chloroplast genome evolution and is a major contributor to its length variation [38].However, the expansion or contraction of IR regions can be attributed to various factors, including genome rearrangement [39], random mutations [40], evolutionary pressure [41], or gene transfer [42].
In conjunction with the current study, we hypothesize that the contraction of the IR region of S. japonicus in this study is primarily propelled by adaptive evolution due to specific environmental stress adaptation.This variation directly contributes to changes in genetic diversity, leading to the generation of new genetic variants.Moreover, it has been demonstrated that the expansion of the IR region can enhance the stability of the chloroplast genome [43].Conversely, the contraction of the IR region in this study indicates a reduction in the structural stability of the chloroplast genome, which may be a primary factor contributing to phenotypic change among varieties.

Positive Selection of ndhD and ycf1
The results of the selection pressure analysis in this study revealed a positive selection of the ndhD and ycf1 genes, directly indicating the adaptive evolution of the four S. japonicus.Positive selection may represent the optimization of gene function in response to environmental stress [44][45][46].
The ndhD-encoded protein is a constituent subunit of the chloroplast NDH complex [47], which plays a crucial role in cyclic electron transport and respiration around photosystem I (PSI), generating a transthylakoid membrane proton gradient and ATP for CO 2 assimilation.The expression of the ndhD gene is subject to regulation through RNA editing.RNA editing can create start codons to control the efficiency of translation and promote gene expression.For instance, in the chloroplast genome of Arabidopsis thaliana, the ACG sequence is altered to AUG, thereby generating the start codon for ndhD [48].Hirose and Sugiura also observed a correlation between the editing efficiency of tobacco ndhD and the stage of leaf development, with the highest efficiency occurring in the early stages and gradually decreasing [49].Additionally, previous research has indicated that environmental factors such as temperature and light can impact the expression of the ndhD gene [50].In this study, we speculate that four S. japonicus ndhD genes were positively selected, leading to alterations in the photosynthetic process and subsequent changes in the morphology of S. japonicus.Although research on this process is limited, these findings directly indicate the ongoing evolution of S. japonicus.
The ycf1 gene is typically the largest in chloroplasts and is a component of the chloroplast inner membrane translocon [51].However, this study found that ycf1 was a pseudogene, a phenomenon also observed in rice, maize, palm, and some Geraniaceae species [52][53][54].Additionally, the ycf1 gene has been actively selected in other lineages, such as Oenothera, Citrus, and Cardamine [55][56][57].If the genes encoded by chloroplasts become pseudogenes, it typically indicates their inability to encode functional proteins.In this study, ycf1 was identified as a pseudogene due to the presence of a premature stop codon.However, we hypothesize that it may still possess certain functionalities, such as (1) influencing the expression of nearby genes due to the presence of enhancers or promoters in the sequence [58]; (2) transcribing into non-coding RNA with regulatory functions [59]; (3) becoming part of the transposable element, affecting the dynamic changes of the genome [60].Of course, pseudogenes may also acquire new functions through evolutionary processes.The positive selection of ycf1 as a pseudogene in this study provides robust support for the ongoing evolutionary event in the S. japonicus chloroplast genome, a phenomenon not previously observed in studies of this species.
In this section, we endeavor to speculate on the underlying factors contributing to chloroplast diversity in S. japonicus by conducting an analysis of ndhD and ycf1.Despite our limited understanding, these findings provide novel insights and establish a groundwork for future genetic manipulation and breeding initiatives.

Phylogenetic Analysis
The phylogenetic structure reconstructed based on the chloroplast genome in this study reveals conflicting relationships within seven S. japonicus samples.This discrepancy, not previously observed in Styrax studies, prompts a discussion of the reasons for these conflicts in the species' phylogeny.
In previous studies, phylogenetic analysis based on chloroplast genomes revealed a close clustering of S. japonicus with S. grandiflorus, S. confuses, and S. calvescens [61], which is similar to the findings of the phylogenetic study of Styrax by Song et al. [62].However, these analyses did not capture the full extent of chloroplast genome diversity in individual species.Our study, which combined analysis of published S. japonicus chloroplast genomes with four newly sequenced samples, revealed that S. japonicus forms at least three distinct clusters within Styrax; it shares sister branches with S. calvescens, S. dasyanthus, and S. faberi, respectively.
Generally, samples of the same species may diverge into different branches for several reasons.Firstly, genetic diversity can lead to significant variations within the species [63].Secondly, geographic isolation or differences in flora distribution between populations may restrict gene flow [64].Thirdly, niche differentiation can promote the development of various adaptive evolutionary paths within a species [65,66].Fourth, mutations and gene recombination in sexual reproduction can introduce new genetic variations [67,68].Additionally, population bottleneck events or founder effects of small populations can result in a reduction in genetic diversity, reflected in the phylogenetic tree where some samples may be markedly different from others [69][70][71].Furthermore, hybridization and incomplete lineage sorting with closely related species may also contribute to the formation of distinct clades in phylogenetic analysis [16,17,72].
In combination with this study, we inferred that the broad distribution of S. japonicus and its overlapping habit with other Styrax species may increase the likelihood of hybridization and gene introgression, resulting in direct phylogenetic conflict and diversity of chloroplast genomes [73][74][75].For example, Wang et al. conducted a comparison between domesticated soybean and wild-type soybean, revealing that gene introgression events contributed to the genetic diversity of domesticated soybean, as well as the asymmetric evolution of the nuclear genome and chloroplast genome [76].
Hybridization and gene introgression can result in recombination and mutation of the chloroplast genome, leading to the emergence of new chloroplast genotypes with potentially distinct biological functions and adaptations [77,78].For instance, variations in chloroplast genes may impact photosynthesis efficiency, subsequently influencing plant growth and survival [79,80].Gene introgression could also alter the selection pressure on certain genes within the chloroplast genome, thereby facilitating positive or negative selection processes and impacting the adaptive evolution of the chloroplast genome [81].
The complexity of relationships among the seven S. japonicus samples in this study indicates that the chloroplast genome of S. japonicus exhibits diversity and may be undergoing evolutionary divergence.Although the evolutionary mechanism has not been fully elucidated, the findings of this study provide compelling evidence for the evolution of the chloroplast genome of S. japonicus and establish a solid foundation for future germplasm resources.

Plant Materials
In this study, fresh blades of four samples were collected from the natural population of S. japonicus in Weide Mountain, Rongcheng, Shandong Province, China (122.422• E, 37.273 • N).
We designated them as Jian Ren (JR), Lian Yi (LY), Yan Zi (YZ), and the wild type (WT) for differentiation.To be specific, three varieties exhibited morphological variations in leaf color, shape, and texture when compared to the WT (Figure 6A).The leaves of JR are thick and leathery (Figure 6B), those of LY exhibit leaf folds (Figure 6C), and the leaves of YZ are purplish red in color (Figure 6D).The voucher specimens were deposited in the herbarium of Nanjing Forestry University.The methods are exclusively conducted on Styrax plants for experimental purposes, strictly adhering to relevant institutional, national, and international guidelines and regulations.

DNA Extraction and Sequencing
The fresh and healthy leaf tissues of four samples were frozen in liquid nitrogen and stored in an Ultra-low temperature freezer for DNA extraction.The DNA extraction was carried out using the Plant Genomic DNA Kit (Nanjing Genepioneer Biotechnologies Inc., Nanjing, China) according to the manufacturer's instructions.The concentration and purity of the DNA were evaluated using a Nandrop 2000 instrument (Thermo Fisher Scientific, Waltham, MA, USA).The nucleic acids utilized for sequencing demonstrated a purity range of 1.8-2.0,as indicated by A260/A280 ratios and sample concentrations ≥10 ng/µL.After the quality inspection of the genomic DNA was performed, the specific process included DNA sample detection, fragmentation, end repair, 3 ′ end plus A, adapter ligation, fragment selection by agarose gel electrophoresis, polymerase chain reaction (PCR) amplification for sequencing library formation, and quality verification.Subsequently, it underwent sequencing on the Illumina NovaSeq 6000 (Illumina, Cambridge, MA, USA) platform with a read length of 150 bp.The total DNA extracting and sequencing were conducted by Nanjing Genepioneer Biotechnologies Inc. (Nanjing, China).

Chloroplast Genome Comparative Analysis
The boundaries of the plastome sequences collected and generated in this study were visualized in IRscope [91] (http://genocat.tools/tools/irscope.html, accessed on 1 June 2024).The Shufe-LAGAN model in the mVISTA program was utilized to compare the chloroplast genome structure of three varieties, with WT serving as a reference.DnaSP version 6 [92] was utilized to compute the nucleotide diversity (Pi) of coding sequences (CDS) and complete chloroplast genomes.The analysis employed a step size of 50 bp with a window length of 600 bp.

Repeat and Analysis of Codon Usage Bias
REPuter [93] (https://bibiserv.cebitec.uni-bielefeld.de/reputer,accessed on 1 June 2024) was employed to identify the dimensions and positions of forward, reverse, palindromic, and complementary repeats.Simple sequence repeats (SSRs) were determined using a Perl script MISA [94], encompassing mono-, di-, tri-, tetra-, penta-, and hexanucleotides with minimum thresholds of 10, 6, 5, 5, 5, and 5, respectively.For the analysis of codon bias, full-length CDS sequences were extracted for the three variants using ATG as the start codon and TGA, TAG, or TAA as the stop codon to determine codon usage preferences.CodonW was utilized to calculate the nucleotide composition at the third position (A3s, U3s, and G3s), codon adaptation index (CAI), codon bias index (CBI), and effective number of codons (ENC).

RSCU Analysis and Selective Pressures in the Evolution
Relative synonymous codon usage (RSCU) was analyzed using BioPython 1.84 [95].The synonymous (Ks) and non-synonymous (Ka) substitution, as well as their ratios (Ka/Ks), were calculated by extracting and comparing the common coding sequences (CDS) of the four samples using Mafft 7.463 [96].The YN algorithm was chosen in the KaKs calculator 3.0 [97] to consider the evolutionary characters of the sequence [14,98] for the purpose of elucidating the Ka/Ks value and conducting selection pressure analysis.

Phylogenetic Analysis
For the purpose of investigating the phylogenetic positions of the three Styrax varieties, a total of 31 Styrax chloroplast genomes were utilized for phylogenetic analysis, with two additional Sinkjackia species (Sinojackia xylocarpa and Sinojackia sarcocarpa) serving as outgroups.
In order to explore the phylogenetic relationships of S. japonicus, 31 Styrax contained 3 additional complete chloroplast genomes of S. japonicus (obtained from NCBI, accession number: MZ287543; NC047429; MT187456).Two datasets (complete chloroplast genomes and CDS sequences) were generated from chloroplast genomes, and after alignment, the optimal model was determined using ModelFinder [99].Maximum-likelihood (ML) phylogenetic analysis was conducted using IQ-TREE 2.2.0 [100] with 1000 bootstrap replicates.

Conclusions
To explore the diversity of the chloroplast genome of S. japonicus and its effects on its morphology, we conducted the complete chloroplast genome sequencing of four samples from a natural population in Shandong, China, comprising three variants and one wild-type sample.Comparison with the previously published chloroplast genome of S. japonicus revealed a slight contraction in the IR region of four samples in this study, which may be indicative of adaptive evolution and indirectly suggests clues to genetic diversity.Both ycf1 and ndhD underwent positive selection, with ycf1 being identified as a pseudogene, indicating the adaptive evolution of S. japonicus to environmental stress.Phylogenetic analysis revealed that S. japonicus did not exhibit clustering or sister clades, suggesting that genetic diversity, geographic isolation, interbreeding, and other factors may have contributed to this phenomenon.In conclusion, a comparative analysis of the chloroplast genome suggests that S. japonicus may be undergoing divergent evolutionary paths, primarily driven by adaptive evolution.Furthermore, the observed variations in leaf morphology provide evidence for this phenomenon.While morphological changes are primarily regulated by the nuclear genome, the evolution of the chloroplast genome indirectly contributes to the development of these changes to some extent.Although the specific mechanisms underlying leaf deformation remain unknown, the diversity present in the chloroplast genome lays a foundation for future breeding and gene editing endeavors.
Step 1: The seed sequence of the chloroplast genome was obtained by assembling the cpDNA sequences with SPAdes v3.10.1.
Step 2: kmer iterative extend seed if the result of step 2 was a contig, the result was determined to be a pseudogenome sequence, and step 6 was performed directly.
Step 4: Use Gapfiller v2.1.1 [85] to make up the gaps for the scaffold sequences obtained in step 3.
Step 5: If a gap still exists after these operations, design primers, conduct PCR sequencing, and reassemble until the complete pseudogenome sequence is obtained.
Step 6: The sequencing reads were aligned to the pseudogenome to correct the genome.
Step 7: According to the structure of the chloroplast, the corrected pseudogenome was rearranged to obtain the complete chloroplast circular genome sequence.

Figure 1 .
Figure 1.Comparison of chloroplast genome structure in four S. japonicus samples.IR (inverted repeat), LSC (large single copy), and SSC (small single copy) regions and border genes are indicated.Note: JLA: junction between LSC and IRa; JLB: junction between LSC and IRb; JSA: junction between SSC and IRa; JSB: junction between SSC and IRb.

Figure 2 .
Figure 2. Nucleotide diversity of chloroplast genomes in S. japonicus.(A): Pi in CDS; (B): chloroplast genome Pi value.Note: window length: 600 bp; step length: 50 bp; x axis: position of the midpoint of each window; y axis: Pi of each window.The blue line represents the trajectory of the value of Pi.

Figure 3 .
Figure 3. Analysis of SSR sites and repetitive sequences in five chloroplast genomes.(A): Distribution of SSRs in the four S. japonicus samples; (B): Number of different SSRs loci types; (C): number of different repeat types.Note: In (A), symbol (+) represented the position of SSRs and the proportion of text displayed; in (C), C: complementary repeats, F: forward repeats, P: palindromic repeats, R: reverse repeats.

Figure 4 .
Figure 4. Synonymous codon usage and selective pressures in the evolution.(A): Codon content of 20 amino acids and stop codons in all protein-coding genes of four S. japonicus chloroplast genome; (B): Distribution of codon preference in S. japonicus; (C): Ka/Ks values of protein-coding genes of the five comparative combinations.Note: in (A), the top panel shows the RSCU for the corresponding amino acids, and the colored blocks, which are shown below, represent different codons.In (C), Ka: nonsynonymous; Ks: synonymous.

Figure 5 .
Figure 5. Phylogenetic tree analysis using the maximum likelihood (ML).(A) Phylogenetic topologies of 33 Styrax species; (B) phylogenetic tree constructed from complete chloroplast genome sequence (left) and CDS sequence (right).

Figure A1 .
Figure A1.Gene map of the four S. japonicus chloroplast genomes.

Figure A2 .
Figure A2.Variation level of three varieties of chloroplast genome sequences.The y-axis indicates the level of variation (between 50 and 100%), and the x-axis represents the coordinate in the chloroplast genome.

Table A1 .
Characteristics of assembled chloroplast genome basic structure.

Table A2 .
Lists of genomic genes for Jian Ren, Lian Yi, Yan Zi, and wild type.Notes: Gene*: Gene with one introns; Gene**: Gene with two introns; #Gene: Pseudogene; Gene(2): Number of copies of multi-copy genes.