Next Article in Journal
Genetic Analysis of Novel Behaviour Traits in Pigs Derived from Social Network Analysis
Next Article in Special Issue
The Complete Chloroplast Genomes of Primula obconica Provide Insight That Neither Species nor Natural Section Represent Monophyletic Taxa in Primula (Primulaceae)
Previous Article in Journal
Variable Clinical Appearance of the Kir2.1 Rare Variants in Russian Patients with Long QT Syndrome
Previous Article in Special Issue
Comparative Chloroplast Genome Analysis of Wax Gourd (Benincasa hispida) with Three Benincaseae Species, Revealing Evolutionary Dynamic Patterns and Phylogenetic Implications
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 4
10.3390/genes13040560
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species

by
Xun Wang
1,
Daru Wang
1,
Ning Gao
2,
Yuepeng Han
3,
Xiaofei Wang
1,
Xiang Shen
4,* and
Chunxiang You
1,*
1
State Key Laboratory of Crop Biology, College of Horticultural Science and Engineering, Shandong Agricultural University, Taian 271018, China
2
Institute of Advanced Agricultural Sciences, Peking University, Weifang 261325, China
3
Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture, Wuhan Botanical Garden, Chinese Academy of Sciences, Wuhan 430074, China
4
National Research Center for Apple Engineering and Technology, Shandong Collaborative Innovation Center for Fruit and Vegetable Production with High Quality and Efficiency, Shandong Agricultural University, Taian 271018, China
*
Authors to whom correspondence should be addressed.
Genes 2022, 13(4), 560; https://doi.org/10.3390/genes13040560
Submission received: 3 March 2022 / Revised: 19 March 2022 / Accepted: 19 March 2022 / Published: 22 March 2022
(This article belongs to the Special Issue Advances in Evolution of Plant Organelle Genome)
Browse Figures
Versions Notes

Abstract

:
The genus Malus is rich in species and many of its plastid genomes have been released. However, limited resources and few markers are not conducive to the comparison of differences among species and resource identification and evaluation. In this study, the complete chloroplast genome of Malus zhaojiaoensis was studied by NGS sequencing, with a total length of 159998 bp. It consists of four regions, LSC (88,070 bp), IRB (26,359 bp), SSC (19,210 bp) and IRA (26,359 bp). M. zhaojiaoensis cp genome contained a total of 111 genes made up of three classes: 78 coding sequences, 29 tRNA genes, and four rRNA genes. In addition, a total of 91 SSRs and 43 INEs were found in the M. zhaojiaoensis cp genome, which was slightly different from M. baccata and M. hupehensis in number. The analysis of codon usage and RNA editing showed that high-frequency codons tended to end at A/U bases and RNA editing mainly occurred at the second codon. Comparative genome analysis suggested that the cp genomes of eight Malus species had higher overall similarity, but there were more variation hotspots (rps16_trnK-UUU, trnG-UCC_atpA, atpH_atpF, trnT-GGU_psbD, etc.) in the LSC region. By building evolutionary trees, it can be clearly observed that M. zhaojiaoensis formed a large group with eight species of Malus, but was relatively independent in differentiation. In conclusion, this study provides high-quality chloroplast genome resources of M. zhaojiaoensis and discusses the genetic variation characteristics of Malus genus. The findings of this study will provide a good reference for plastid genome assembly and interspecific comparison in the future.

1. Introduction

The genus Malus is rich in species and widely distributed and cultivated throughout the world [1]. According to the principle of ecological geography, Malus genus can be divided into two groups: cultivated species and wild species [2]. Among them, cultivated species are greatly affected by human activities and are mainly used for fresh food and processing. Of course, some species are also used as ornamental plants [3]. Wild species exist in natural distribution areas and are often used as rootstock resources and breeding materials in production due to their good adaptability of introduction and affinity for grafting [4]. As a wild species of Malus native to China and widely distributed in northeast, northern and southwest China, M. baccata (L.) Borkh. has diverse populations and variation types [5], which is of great value in producing activity and scientific research. M. zhaojiaoensis Jiang is a new species of Malus found in Zhaojue Xian, Sichuan Province [6]. According to the report [6], M. zhaojiaoensis has oval leaves, reddish flowers, and small fruits, which are of great ornamental value. It has similar characteristics with M. baccata and M. rockii [7], but its evolutionary status and genetic relationship is still ambiguous and needs further research.
There are many methods to study plant genetic evolution, including morphological comparison, biochemical methods (isozymes, etc.) and a wide range of molecular markers (Restriction Fragment Length Polymorphism, Amplified Fragment Length Polymorphism, Random Amplified Polymorphism DNA, Simple Sequence Repeat, Single Nucleotide Polymorphism) [8]. In recent years, with the emergence of high-throughput sequencing technology and the development of bioinformatics based on multiple omics, methods such as resequencing, comparative genomics, and chloroplast genome assembly provide new and efficient solutions for species classification and variation detection [9]. Chloroplast DNA exists in plant chloroplasts in a double chain ring shape, and most plant chloroplast genomes consist of four typical parts, namely a large single copy region (LSC), a small single copy region (SSC), and two inverted repeat regions (IRA and IRB) [10]. Compared with nuclear genomes, chloroplast genomes are relatively conservative and small, so they play an irreplaceable role in species identification [11].
Up to now, there have been numerous studies on chloroplast genomes of Malus plants, such as M. toringoides [12], M. kansuensis [13], M. sieboldii [14], M. sylvestris [15], M. hupehensis [16], and M. prattii [17]. This has played a positive role in the utilization and protection of Malus resources. In this paper, the complete chloroplast genome of M. zhaojiaoensis was assembled based on second-generation sequencing data and bioinformatics methods, and its basic characteristics, sequence similarity, and evolutionary relationship were analyzed and compared. The above results can provide reference for explaining the genetic variation characteristics of M. zhaojiaoensis and Malus. At the same time, the findings of this paper will play a positive role in the conservation and development of wild species such as M. zhaojiaoensis.

2. Materials and Methods

2.1. Sample Collection and DNA Sequencing

Young, fresh and clean leaves of M. zhaojiaoensis Jiang tree were collected from Fruit Research Institute, Chinese Academy of Agricultural Sciences (Xingcheng city, Liaoning Province, China). Total DNA was extracted by improved CTAB (Cetyltrimethylammonium Bromide) method [18]. Construction of genomic DNA library and DNA sequencing were performed using Illumina HiSeq X Ten system [19]. 150 bp paired-end sequences were sequenced in Illumina platform, and insert size was 350 bp. Adapter sequences were removed by Trimmomatic v0.39 software and adapter-trimmed raw reads were used for subsequent assembly analysis [20].

2.2. Chloroplast Genome Assembly and Annotation

GetOrganelle v1.7.5.0 [21] and NOVOplasty4.2 [22] were used for chloroplast genome assembly. In GetOrganelle, the values of kmers were set to 21, 45, 65, 85, 105. The complete chloroplast genome sequence of M. zhaojiaoensis was obtained after the assembly results were confirmed. After that, the sequence was submitted to PGA (Plastid Genome Annotator) program [23] for gene annotation, and the annotation information of CPGAVAS [24] and GeSeq [25] was also used for integration and merger. Chloroplast genome sequence and annotation file were manually revised and submitted to NCBI website (https://www.ncbi.nlm.nih.gov/, accessed on 20 February 2022). Its accession number was OM793283. OrganellarGenomeDRAW [26] was used to graphically display the annotations of cp genome, and different genes were distinguished by different colors.

2.3. Description of Basic Characteristics of Chloroplast Genome

The basic characteristics of chloroplast genome include size and structure, GC content, gene classification and repeat sequence, and so on. The above analysis was accomplished by Python (https://www.python.org/, accessed on 10 February 2022), Geneious (https://www.geneious.com/, accessed on 15 February 2022), MISA (https://webblast.ipk-gatersleben.de/misa/, accessed on 8 February 2022) [27], Vmatch (http://vmatch.de/, accessed on 8 February 2022), etc. Regarding search parameters of simple sequence repeats (SSRs), the minimum number of repetitions was set to 10 (mononucleotide), 5 (dinucleotide), 4 (trinucleotide), 3 (tetranucleotide), 3 (pentanucleotide), and 3 (hexanucleotide). In the identification of interspersed nuclear elements (INEs), the parameters were as follows: Minimal Repeat Size (30) and Hamming Distance (3). The cp genomes of M. baccata (MK896774) and M. hupehensis (MK020147) were used for comparative analysis.

2.4. Analysis of Codon Usage Characteristics and RNA Editing Sites

All the coding sequences were extracted based on the chloroplast genome annotation of M. zhaojiaoensis. The calculation and comparison of RSCU (relative synonymous codon usage) value was operated in CodonW v1.4.2 (http://codonw.sourceforge.net/, accessed on 12 February 2022) programme. PREP-Cp (http://prep.unl.edu/, accessed on 16 February 2022) was used to predict potential RNA editing sites in coding sequences [28].

2.5. Chloroplast Genome Alignment, Boundary and Evolutionary Analysis

Several chloroplast genomes of Malus and other genera were searched from NCBI nucleotide database, including M. hupehensis (MK020147), M. baccata (MK896774), M. halliana (MT246302), M. micromalus (MF062434), M. prattii (MH929090), M. prunifolia (KU851961), M. sieboldii (MT593044), M. toringoides (MT483999), M. yunnanensis (MH394387) and Crataegus hupehensis (Crataegus, MW201730), etc. The homology and similarity of the cp genomes of four Malus genera were compared in Circoletto (http://tools.bat.infspire.org/circoletto/, accessed on 18 February 2022) [29]. The sequences of chloroplast genomes were aligned in VISTA (https://genome.lbl.gov/vista/index.shtml, accessed on 16 February 2022); Shuffle-LAGAN was designated as alignment program [30]. IRscope was used to visualize the junction sites of the chloroplast genomes [31]. Phylogenetic and cluster analysis was performed in HomBlocks (https://github.com/fenghen360/HomBlocks, accessed on 18 February 2022) [32] and MEGA tool (https://www.megasoftware.net/, accessed on 18 February 2022) [33], and C. hupehensis (Crataegus) was selected as outgroup. Two types of evolutionary trees (Neighbour-Joining tree and Maximum Likelihood tree) were constructed in MEGA X, both of which were examined with bootstrap method (1000 replications). In addition, ClustalW (Codons) alignment method and Jukes-Cantor substitution model were selected for NJ evolutionary trees based on single-copy genes (matK and rbcL).

3. Results

3.1. Characteristics of the Chloroplast Genome of M. zhaojiaoensis

3.1.1. Chloroplast Genome Map and Structure

After genome sequencing and assembly, a complete chloroplast genome sequence of M. zhaojiaoensis was generated, with a total length of 159,998 bp and an overall GC content of 36.6%. Like other chloroplasts, the chloroplast genome of M. zhaojiaoensis has a double chain ring structure and includes four distinct and typical regions (two single copy regions and two inverted repeat regions). The sequence lengths of the tetrads were 88,070 bp (LSC: 1–88,070), 26,359 bp (IRB: 88,071–114,429), 19,210 bp (SSC: 114,430–133,639) and 26,359 bp (IRA: 133,640–159,998), respectively (Figure 1). By statistical calculation, GC content in IR regions (42.7%) was higher than that in the whole cp genome (36.6%) of M. zhaojiaoensis, indicating that the two inverted repeat regions (IRA and IRB) were relatively stable.

3.1.2. Gene Annotation and Classification of Chloroplast Genome

Genetic annotations were completed on the assembled genome of M. zhaojiaoensis, which contains a total of 111 genes made up of three classes: CDS (coding sequence, 78), tRNA (transfer RNA, 29) and rRNA (ribosomal RNA, 4). According to function, all annotated genes could be divided into four main categories, photosynthesis, self-replication, others, and unknown (Table 1). Six atp (subunits of ATP synthase), 15 psb (subunits of photosystem II), 11 ndh (subunits of NADH-dehydrogenase), 6 pet (subunits of cytochrome b/f complex), 5 psa (subunits of photosystem I) and 1 rbc (subunit of rubisco) genes were involved in photosynthesis. Nine rpl (large subunit of ribosome), 4 rpo (DNA dependent RNA polymerase), 12 rps (small subunit of ribosome), 4 rrn (rRNA), and 29 trn (tRNA) genes were associated to self-replication. In addition to the above two parts, there were nine genes belonging to other and unknown classes, namely accD, ccsA, cemA, clpP, matK and ycf (ycf1, ycf2, ycf3, ycf4).
Among 111 genes, six CDS (ndhB, rpl2, rpl23, rps7, rps12, ycf2), eight tRNA (trnA-UGC, trnI-CAU, trnI-GAU, trnL-CAA, trnM-CAU, trnN-GUU, trnR-ACG, trnV-GAC), four rRNA (rrn4.5, rrn5, rrn16, rrn23) genes had two copies, and the rest had only one copy. Due to the existence of duplicate genes, the actual number of genes in the chloroplast genome of M. zhaojiaoensis reached 129, including 84 encoding genes, 37 tRNA and 8 rRNA genes.
Table
Table 1. Gene classification of chloroplast genome in M. zhaojiaoensis.
Table 1. Gene classification of chloroplast genome in M. zhaojiaoensis.
Gene FunctionGene GroupCodeGene NumberGene Name
PhotosynthesisSubunits of ATP synthaseatp6 (6)atpA, atpB, atpE, atpF, atpH, atpI
Subunits of photosystem IIpsb15 (15)psbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ
Subunits of NADH-dehydrogenasendh11 (12)ndhA, ndhB (×2), ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK
Subunits of cytochrome b/f complexpet6 (6)petA, petB, petD, petG, petL, petN
Subunits of photosystem Ipsa5 (5)psaA, psaB, psaC, psaI, psaJ
Subunit of rubiscorbc1 (1)rbcL
Self replicationLarge subunit of ribosomerpl9 (11)rpl14, rpl16, rpl2 (×2), rpl20, rpl22, rpl23 (×2), rpl32, rpl33, rpl36
DNA dependent RNA polymeraserpo4 (4)rpoA, rpoB, rpoC1, rpoC2
Small subunit of ribosomerps12 (14)rps11, rps12 (×2), rps14, rps15, rps16, rps18, rps19, rps2, rps3, rps4, rps7 (×2), rps8
rRNArrn4 (8)rrn4.5 (×2), rrn5 (×2), rrn16 (×2), rrn23 (×2)
tRNAtrn29 (37)trnA-UGC (×2), trnC-GCA, trnD-GUC, trnE-UUC, trnF-GAA, trnG-GCC, trnG-UCC, trnH-GUG, trnI-CAU (×2), trnI-GAU (×2), trnK-UUU, trnL-CAA (×2), trnL-UAA, trnL-UAG, trnM-CAU (×2), trnN-GUU (×2), trnP-UGG, trnQ-UUG, trnR-ACG (×2), trnR-UCU, trnS-GCU, trnS-GGA, trnS-UGA, trnT-GGU, trnT-UGU, trnV-GAC (×2), trnV-UAC, trnW-CCA, trnY-GUA
OthersSubunit of Acetyl-CoA-carboxylaseacc1 (1)accD
c-type cytochrom synthesis geneccs1 (1)ccsA
Envelop membrane proteincem1 (1)cemA
Proteaseclp1 (1)clpP
Maturasemat1 (1)matK
UnknownConserved open reading framesycf4 (5)ycf1, ycf2 (×2), ycf3, ycf4
The number in parentheses in the fourth column is the total number of duplicates contained, and the fifth column (×2) refers to two copies of the gene.
In the chloroplast genome of M. zhaojiaoensis, most genes do not have introns, but there are still some genes that contain one to two introns (Table 2). The total number of genes with introns was 23 (15 CDS and 8 tRNA genes), of which 19 had one intron, and four (ycf3, rps12, rps12, clpP) had two introns. In terms of genome location and distribution, 12 (eight CDS and four tRNA genes) of the 23 genes were located in LSC region of cp genome, accounting for the largest proportion. Four genes (rpl2, ndhB, trnI-GAU, trnA-UGC) were located in IRA and IRB, and only one gene (ndhA) was distributed in SSC region. In addition, it is worth noting that the rps12 gene encoding ribosomal subunit S12 protein consists of two distant parts of the cp genome, one in the LSC region and the other in the IR region, and it undergoes special trans-splicing during transcription.

3.2. The Type and Distribution of Repeat Sequences

3.2.1. Simple Sequence Repeats

Six types (mono-, di-, tri-, tetra-, penta-, and hexa nucleotide) of SSRs were identified in the MISA-web. In chloroplast genome of M. zhaojiaoensis, single nucleotide repeat (68) was the most frequent SSR (Figure 2A, Table S1), followed by dinucleotide repeat (16). In addition, there were five tetranucleotide SSR markers and two pentanucleotide SSR markers. No trinucleotide repeats and hexanucleotide repeats were found in chloroplast genome of M. zhaojiaoensis. For comparison, SSR types and numbers in M. baccata and M. hupehensis cp genome were also identified and counted in this study. The total number of SSRs in the three cp genomes were similar, 91 in M. zhaojiaoensis, 97 in M. baccata and 96 in M. hupehensis. Unlike M. zhaojiaoensis and M. hupehensis, there was no pentanucleotide repeat type in M. baccata, but the M. baccata cp genome had one hexanucleotide SSR (Figure 2A).
For mononucleotide SSRs, A/T repeat was most abundant with 66, 67 and 66 times in M. zhaojiaoensis, M. baccata and M. hupehensis, respectively. It is noteworthy that the composition of pentanucleotide repeats in M. zhaojiaoensis was AAGGC/CCTTG, and that in M. hupehensis was AAAAT/ATTTT (Figure 2B). The only hexanucleotide repetition in M. baccata was the ATAATT unit repeated three times, with a total base length of 18 bp.

3.2.2. Interspersed Nuclear Elements

Interspersed repeats are another important repeating element in the genome and there are abundant INEs (interspersed nuclear elements) in the chloroplast genome of plants. Four types (forward_direct, reverse, complement, palindromic) of INEs were calculated in REPuter. In the cp genome of M. zhaojiaoensis, the INEs contained 21 forward repeats, 20 palindromic repeats, 1 reverse repeat, and 1 complement repeat (Figure 2C, Table S2). The M. baccata cp genome had 27 F, 19 P, 6 R, and 1 C type repeats, while there were 29 forward, 23 palindromic, 5 reverse, and 0 complement repeats in M. hupehensis (Figure 2C). In general, the number of INEs in M. zhaojiaoensis (43) was significantly less than that in M. baccata and M. hupehensis, and the total number of INEs in M. baccata (53) and M. hupehensis (57) was close.
In addition, from the perspective of repeated fragment size, it can be found that most of the sequence lengths were concentrated in the range of 30–50 bp, especially 30–40 bp. There was one INE repeat with a sequence length of more than 50 bp in each of the chloroplast genomes of three Malus species (Figure 2D).

3.3. Relative Synonymous Codon Usage

The relative synonymous codon usage (RSCU) can represent the relative probability of using synonymous codons and reflect the usage bias of different codons [34]. The total number and relative frequency of all codons were obtained in CodonW software. The results showed that leucine was the most frequent amino acid in the chloroplast coding sequence of M. zhaojiaoensis, followed by isoleucine, serine, glycine and arginine. In leucine, codon UUA was widely used. In alanine, the codon GCU was the most frequently used. The arginine bias uses AGA as the codon, and the termination codon bias ends with UAA (Figure 3). There were 30 kinds of codons whose RSCU value was greater than 1, including GCA (Ala), CCA (Pro), UCA (Ser), AGU (Ser), ACA (Thr), UUG (Leu), CUU (Leu), CGU (Arg), UUU (Phe), GGU (Gly), etc. A total of 29 of the 30 codons mentioned above ends with an A/U base, indicating that the third bases of the high-frequency codons are biased toward A/U. In addition, the usage of codons such as AGC, CUG, GAC, UAC and CUC in the cp genome of M. zhaojiaoensis was relatively low.

3.4. RNA Editing Sites

RNA editing is a condition in which post-transcribed coding sequence undergoes base conversion that results in changes in amino acids. RNA editing also occurs in plant chloroplast encoded proteins, which can participate in the expression and regulation of genes, and affect and change gene function [35]. In this study, all potential RNA editing sites in the coding sequences of M. zhaojiaoensis cp were predicted. The results indicated that a total of 63 RNA editing sites were distributed on 25 genes. These genes mainly included ndh (ndhA, ndhB, ndhD, ndhF, ndhG), matK, rpo (rpoA, rpoB, rpoC1, rpoC2), accD, atp (atpA, atpB, atpF, atpI), pet (petB, petG), rps (rps2, rps14, rps16), etc. ndhB had the most editing loci with 12, followed by ndhD (8), matK (5), ndhF (4), rpoB (4) and rpoC2 (4) (Figure 4). In terms of the location of RNA editing sites, 48 sites were located at the second base of the codon, accounting for 76.19% of all editing sites (Table S3). In conclusion, in the chloroplast coding genes of M. zhaojiaoensis, RNA editing mainly occurred at the second codon. Interestingly, in the CDS of M. zhaojiaoensis cp, RNA editing was the most likely to lead to transitions from Serine (S) to Leucine (L), followed by Proline (P) to Leucine (L), and Serine (S) to Phenylalanine (F) (Table S3).

3.5. Alignment of Chloroplast Genomes in Malus Species

In order to compare the sequence similarities and differences of chloroplast genomes of different Malus genus, M. zhaojiaoensis and three Malus species were used for global alignment (Figure 5). The chloroplast sequences of M. zhaojiaoensis showed good collinearity with three Malus plants. In contrast, the similarity of M. zhaojiaoensis to M. baccata and M. hupehensis was higher, but the similarity to M. yunnanensis was slightly lower. To find the variation hotspot regions of cp genomes, eight species of Malus (including M. zhaojiaoensis) were used for comparative analysis where M. hupehensis (MK020147) was a reference. The results showed that most of the differences were located in the LSC region of chloroplast genome, and the other regions were highly similar. It can be clearly observed from Figure 6, that there are large variations in rps16_trnK-UUU, trnG-UCC_atpA, atpH_atpF, trnT-GGU_psbD, psbZ_trnfM-CAU, trnV-UAC_ndhC, accD_psaI, rps3_rpl16. In general, the hotspots were more distributed in non-coding regions, while gene regions were relatively stable and conserved.

3.6. Comparison of Chloroplast Genome Boundaries and Junction Sites

JLB (LSC-IRB), JSB (SSC-IRB), JSA (SSC-IRA), JLA (LSC-IRA) are the junction sites on the boundaries of the four regions of the chloroplast genome, which are of great significance in cp genome evolution [36]. In this study, cp genome boundaries of eight Malus species were compared in IRscope. As shown in Figure 7, M. prattii has the largest cp genome with 160,239 bp, while M. zhaojiaoensis cp genome has the smallest sequence length of 159,998 bp. It can be found that the length of LSC region varies from 88,070 to 88,355 bp in eight cp genomes. In contrast, the sequence length of SSC and two IR regions varies little.
For rps19 gene at the JLB boundary, its location distribution (159 bp in LSC and 120 bp in IRB) was consistent in seven Malus species, but the length on the LSC side_rps19 (210 bp) was longer in M. yunnanensis, and sequence length on the IRB side_rps19 decreased (69 bp). At the JSB site of chloroplast genome of M. prunifolia and M. baccata, the ycf1 gene expanded to the SSC region, resulting in the increased length of ycf1 gene. The ycf1 gene of IRB/SSC locus did not cross the JSB boundary in M. toringoides, M. sieboldii, M. prattii and M. hupehensis. In addition, the displacement of trnH gene at JLA locus also occurred in different species (Figure 7). The above results prove that boundary region expansion and contraction play a role in chloroplast genome size and construction.
Genes 13 00560 g007 550
Figure 7. Junction sites and boundary characteristics of eight chloroplast genomes. JLB, JSB, JSA and JLA represent four different junction sites in cp genome boundaries.
Figure 7. Junction sites and boundary characteristics of eight chloroplast genomes. JLB, JSB, JSA and JLA represent four different junction sites in cp genome boundaries.
Genes 13 00560 g007

3.7. Phylogenetic Relationships Based on Chloroplast Genomes

Comparison of chloroplast genome sequences can provide a basis for explaining evolutionary relationships of species. In order to study the evolutionary status of M. zhaojiaoensis, chloroplast genomes of 10 species were downloaded from NCBI, including one outgroup (C. hupehensis, Crataegus). Based on the comparison of 98,748 locations in the dataset, NJ and ML trees for 11 species were constructed. Because the two types of trees were detected by bootstrap method, and their topological structures were very consistent in the results, the evolutionary branches were trustworthy. As can be seen from Figure 8, 11 species were divided into two groups according to genus, and Malus plants were grouped into one group. In the Malus genus, M. yunnanensis was an independent clade due to sequence differences, while M. zhaojiaoensis and other eight species form a large clade. The figures in the evolutionary branching shown high support for this result (Figure 8).
matK and rbcL, as single-copy genes in chloroplast genome, evolved at a moderate rate and were both suitable for serving as barcodes for plant classification [37]. Hence, evolutionary relationships were characterized based on single-copy chloroplast genes matK and rbcL from 11 species (Figure 9). The results showed that M. zhaojiaoensis and eight species of Malus were grouped together, and the clades formed by the two genes were slightly different.

4. Discussion

With the development of high-throughput sequencing, more and more chloroplast genomes have been released [38]. Chloroplast genome is different from nuclear genome in that it has the characteristics of maternal inheritance. Because of its small sequence length and moderate base replacement, chloroplast genome has been widely used in the study of genetic variation and phylogeny [39].
In this study, the complete chloroplast genome of M. zhaojiaoensis was constructed by whole genome sequencing. Its cp size (159,998 bp) was similar to that of reported Malus species, including M. baccata (160,024 bp), M. hupehensis (160,065 bp), M. sieboldii (160,040 bp) and M. prunifolia (160,041 bp), etc. This indicates that the length of interspecific chloroplast sequence within Malus genus was relatively conserved and the variation was relatively moderate [40]. Further, by annotating the cp sequence of M. zhaojiaoensis, it was found that the ring-structured genome consists of four parts, namely LSC, SSC, IRA and IRB. Comparative analysis of different species showed that the LSC regions of the Malus cp genomes are highly variable, which contain many hot spots, for example, trnG-UCC_atpA and trnT-GGU_psbD. The proportion of SSR (75.95%) and INE (60.47%) repeat markers in the LSC region of M. zhaojiaoensis was larger than that of the other three regions, which also proved the above inference.
For chloroplast genes, 78 CDS, 29 tRNA, and 4 rRNA were annotated in M. zhaojiaoensis. But there were 111 genes in M. baccata cp genome, which included 76 protein coding genes, 31 tRNA genes and 4 rRNA genes. In addition, the number of tRNA in chloroplast of M. hupehensis was one more than that of M. zhaojiaoensis. One of the most important factors contributing to the difference in gene numbers between species may be the expansion and contraction of chloroplast genome boundaries [41]. For example, comparative analysis showed that the rps19 gene at the JLA junction site was missing in three Malus plants (M. zhaojiaoensis, M. prunifolia and M. yunnanensis).
Chloroplast coding sequences were evolutionarily conserved due to low variation. The codons of chloroplast CDS of M. zhaojiaoensis were calculated and the results suggested that the high frequency codons were mainly UUA (Leu), GCU (Ala), AGA (Arg) and so on. In Ocotea aciphylla, there were 30 high-frequency codons whose RSCU value was greater than 1, including AGA (Arg), GGA (Gly) and UCA (Ser), etc. [42]. In M. zhaojiaoensis cp, 29 of the 30 high-frequency codons ended with an A/U base, which reflects the biased use of codons. The similar phenomenon has been found in previous reports [43]. RNA editing was widely used in higher plants and was also an important means of regulating chloroplast gene expression. A total of 63 loci were identified in the coding sequences of M. zhaojiaoensis cp and distributed in 25 genes. ndh (Subunits of NADH-dehydrogenase) genes had the most potential RNA editing sites (27), followed by rpo (DNA dependent RNA polymerase) genes (10 sites). In addition, the predicted results suggested that RNA editing tends to lead to an increase in hydrophobic amino acids (Leucine, Phenylalanine), which is consistent with the study of Platanthera ussuriensis [44]. In Dipterygium glaucum and Cleome chrysantha cp, the amino acids Serine to Leucine were the majority of the conversion in RNA editing [45], which also supports the above speculation.
Chloroplast genome plays an irreplaceable role in species evolution and classification [46]. The information of chloroplast genome can be used to explain the influence of maternal inheritance [47]. The comparison of genetic evolution in Malus species has been a long-standing concern of researchers [48]. In this study, based on chloroplast genome and two single copy genes, the evolutionary relationship of Malus species was analyzed. The branches and topologies of evolutionary trees were consistent in different methods and datasets. Based on the above results, it can be clearly inferred that M. yunnanensis was located in the base group of Malus [49], and M. zhaojiaoensis formed a large group with 8 species of Malus, but was relatively independent in differentiation. In addition, there was further differentiation between M. toringoides, M. halliana and M. hupehensis, while the other five Malus species also closely influenced each other. The findings of this paper are consistent with the previous conjecture [6], indicating that M. zhaojiaoensis is relatively primitive in evolution.
Due to unique geographical conditions and climate distribution, China is rich in wild plant types and contributes significantly to the world’s biodiversity. M. zhaojiaoensis is a wild Malus species of China. The exploration of biological classification and genetic relationships of M. zhaojiaoensis can speed up the process of conservation and utilization, which provides a good reference for the management and development of wild resources.

5. Conclusions

In this study, the complete chloroplast genome of M. zhaojiaoensis was sequenced and assembled. Based on the annotation of the cp genome, the structure and sequence characteristics of the genome were analyzed. Prediction of repeat sequence types and RNA editing sites provides potential genetic markers and regulatory loci. Hot spots, boundary features, and evolutionary relationships of M. zhaojiaoensis and Malus were described by comparative genomics and phylogenetic analysis. The release of chloroplast genome of M. zhaojiaoensis can provide valuable support for genetic variation and germplasm identification of Malus species. The research in this paper is helpful to the conservation and development of wild Malus resources.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13040560/s1, Table S1: Type and number of SSRs in M. zhaojiaoensis cp genome; Table S2: Type and size of INEs in M. zhaojiaoensis cp genome; Table S3: RNA editing sites in M. zhaojiaoensis cp genome.

Author Contributions

Data curation, X.W. (Xun Wang), D.W., N.G., Y.H. and X.W. (Xiaofei Wang); Formal analysis, D.W. and N.G.; Funding acquisition, X.S. and C.Y.; Investigation, X.W. (Xun Wang), D.W., Y.H. and X.W. (Xiaofei Wang); Methodology, X.W. (Xun Wang) and X.W. (Xiaofei Wang); Project administration, X.S. and C.Y.; Resources, X.W. (Xiaofei Wang); Visualization, X.W. (Xun Wang); Writing—original draft, X.W. (Xun Wang); Writing—review & editing, X.W. (Xiaofei Wang), X.S. and C.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program (2018YFD1000100), General Project of National Natural Science Foundation of China (32172538), National Natural Science Foundation of China (32072520), National Modern Apple Industry Technology System Funded Program of China (CARS-27), Natural Science Foundation of Shandong Province (ZR2020MC132), Improved Seed Industrialization Project of Shandong Province (2019LZGC007), Fruit Industry System of Shandong Province (SDAIT-06-07).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated during this study are included in this published article.

Acknowledgments

We would like to thank the State Key Laboratory of Crop Biology.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, M.-R.; Chen, L.; Teixeira da Silva, J.A.; Volk, G.M.; Wang, Q.-C. Cryobiotechnology of apple (Malus spp.): Development, progress and future prospects. Plant Cell Rep. 2018, 37, 689–709. [Google Scholar] [CrossRef] [PubMed]
  2. Chen, P.; Li, Z.; Zhang, D.; Shen, W.; Xie, Y.; Zhang, J.; Jiang, L.; Li, X.; Shen, X.; Geng, D.; et al. Insights into the effect of human civilization on Malus evolution and domestication. Plant Biotechnol. J. 2021, 19, 2206–2220. [Google Scholar] [CrossRef] [PubMed]
  3. Zhou, T.; Fan, J.; Zhao, M.; Zhang, D.; Li, Q.; Wang, G.; Zhang, W.; Cao, F. Phenotypic variation of floral organs in Malus using frequency distribution functions. BMC Plant Biol. 2019, 19, 574. [Google Scholar] [CrossRef] [PubMed]
  4. He, J.; Zhou, J.; Wan, H.; Zhuang, X.; Li, H.; Qin, S.; Lyu, D. Rootstock–scion interaction affects Cadmium accumulation and tolerance of Malus. Front. Plant Sci. 2020, 11, 1264. [Google Scholar] [CrossRef]
  5. Chen, X.; Li, S.; Zhang, D.; Han, M.; Jin, X.; Zhao, C.; Wang, S.; Xing, L.; Ma, J.; Ji, J.; et al. Sequencing of a wild apple (Malus baccata) genome unravels the differences between cultivated and wild apple species regarding disease resistance and cold tolerance. G3 Genes Genomes Genet. 2019, 9, 2051–2060. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Jiang, N.G. A new species of Malus (Rosaceae) from Sichuan. J. Southwest Agric. Univ. 1991, 13, 599–600. [Google Scholar]
  7. Li, Y.N. A primarily modern systematics of genus Malus Mill. in the world. J. Fruit Sci. 1996, 13, 82–92. [Google Scholar]
  8. Cornille, A.; Giraud, T.; Smulders, M.J.M.; Roldán-Ruiz, I.; Gladieux, P. The domestication and evolutionary ecology of apples. Trends Genet. 2014, 30, 57–65. [Google Scholar] [CrossRef]
  9. Liao, L.; Zhang, W.; Zhang, B.; Fang, T.; Wang, X.-F.; Cai, Y.; Ogutu, C.; Gao, L.; Chen, G.; Nie, X.; et al. Unraveling a genetic roadmap for improved taste in the domesticated apple. Mol. Plant 2021, 14, 1454–1471. [Google Scholar] [CrossRef]
  10. Daniell, H.; Lin, C.-S.; Yu, M.; Chang, W.-J. Chloroplast genomes: Diversity, evolution, and applications in genetic engineering. Genome Biol. 2016, 17, 134. [Google Scholar] [CrossRef] [Green Version]
  11. Daniell, H.; Jin, S.; Zhu, X.-G.; Gitzendanner, M.A.; Soltis, D.E.; Soltis, P.S. Green giant—A tiny chloroplast genome with mighty power to produce high-value proteins: History and phylogeny. Plant Biotechnol. J. 2021, 19, 430–447. [Google Scholar] [CrossRef] [PubMed]
  12. Li, Y.; Liu, Y.; Wu, P.; Zhou, S.; Wang, L.; Zhou, S. The complete chloroplast genome sequence of Malus toringoides (Rosaceae). Mitochondrial DNA B Resour. 2020, 5, 2787–2789. [Google Scholar] [CrossRef] [PubMed]
  13. Li, L.; Ye, Y.; Zhao, M.; Xin, H. Characterization of the complete chloroplast genome sequence of Malus kansuensis (Rosaceae). Mitochondrial DNA B Resour. 2021, 6, 108–109. [Google Scholar] [CrossRef] [PubMed]
  14. Lou, G.; Wang, S.; Zhang, B.; Cheng, Y.; Wang, H. The complete chloroplast genome sequence of Malus sieboldii (Rosaceae) and its phylogenetic analysis. Mitochondrial DNA B Resour. 2020, 5, 2170–2171. [Google Scholar] [CrossRef] [PubMed]
  15. Xu, Y.; Zhao, Y.; Zhao, X.; Chen, X.; Yuan, Z. Characterization of complete chloroplast genome of Malus sylvestris L. Mitochondrial DNA B Resour. 2019, 4, 2357–2358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Zhang, X.; Rong, C.; Qin, L.; Mo, C.; Fan, L.; Yan, J.; Zhang, M. Complete chloroplast genome sequence of Malus hupehensis: Genome structure, comparative analysis, and phylogenetic relationships. Molecules 2018, 23, 2917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Fan, L.; Qin, L.; Yan, J.; Mo, C.; Rong, C.; Meng, Y.; Zhang, M. The complete chloroplast genome sequence of Malus prattii (Rosaceae) and its phylogenetic analysis. Mitochondrial DNA B Resour. 2019, 4, 2171–2172. [Google Scholar] [CrossRef] [Green Version]
  18. Allen, G.; Flores-Vergara, M.; Krasynanski, S.; Kumar, S.; Thompson, W. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat Protoc. 2006, 1, 2320–2325. [Google Scholar] [CrossRef]
  19. Fu, X.; Ji, X.; Wang, B.; Duan, L. The complete chloroplast genome of leguminous forage Onobrychis viciifolia. Mitochondrial DNA B Resour. 2021, 6, 898–899. [Google Scholar] [CrossRef]
  20. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  21. Jin, J.-J.; Yu, W.-B.; Yang, J.-B.; Song, Y.; dePamphilis, C.W.; Yi, T.-S.; Li, D.-Z. GetOrganelle: A fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020, 21, 241. [Google Scholar] [CrossRef]
  22. Dierckxsens, N.; Mardulyn, P.; Smits, G. NOVOPlasty: De novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017, 45, e18. [Google Scholar]
  23. Qu, X.-J.; Moore, M.J.; Li, D.-Z.; Yi, T.-S. PGA: A software package for rapid, accurate, and flexible batch annotation of plastomes. Plant Methods 2019, 15, 50. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, C.; Shi, L.; Zhu, Y.; Chen, H.; Zhang, J.; Lin, X.; Guan, X. CpGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences. BMC Genom. 2012, 13, 715. [Google Scholar] [CrossRef] [Green Version]
  25. Tillich, M.; Lehwark, P.; Pellizzer, T.; Ulbricht-Jones, E.S.; Fischer, A.; Bock, R.; Greiner, S. GeSeq—Versatile and accurate annotation of organelle genomes. Nucleic Acids Res. 2017, 45, W6–W11. [Google Scholar] [CrossRef]
  26. Lohse, M.; Drechsel, O.; Kahlau, S.; Bock, R. OrganellarGenomeDRAW—A suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets. Nucleic Acids Res. 2013, 41, W575–W581. [Google Scholar] [CrossRef]
  27. Beier, S.; Thiel, T.; Münch, T.; Scholz, U.; Mascher, M. MISA-web: A web server for microsatellite prediction. Bioinformatics 2017, 33, 2583–2585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Mower, J.P. The PREP suite: Predictive RNA editors for plant mitochondrial genes, chloroplast genes and user-defined alignments. Nucleic Acids Res. 2009, 37 (Suppl. 2), W253–W259. [Google Scholar] [CrossRef]
  29. Darzentas, N. Circoletto: Visualizing sequence similarity with Circos. Bioinformatics 2010, 26, 2620–2621. [Google Scholar] [CrossRef]
  30. Frazer, K.A.; Pachter, L.; Poliakov, A.; Rubin, E.M.; Dubchak, I. VISTA: Computational tools for comparative genomics. Nucleic Acids Res. 2004, 32 (Suppl. 2), W273–W279. [Google Scholar] [CrossRef]
  31. Amiryousefi, A.; Hyvönen, J.; Poczai, P. IRscope: An online program to visualize the junction sites of chloroplast genomes. Bioinformatics 2018, 34, 3030–3031. [Google Scholar] [CrossRef] [PubMed]
  32. Bi, G.; Mao, Y.; Xing, Q.; Cao, M. HomBlocks: A multiple-alignment construction pipeline for organelle phylogenomics based on locally collinear block searching. Genomics 2018, 110, 18–22. [Google Scholar] [CrossRef] [PubMed]
  33. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  34. Parvathy, S.T.; Udayasuriyan, V.; Bhadana, V. Codon usage bias. Mol. Biol. Rep. 2022, 49, 539–565. [Google Scholar] [CrossRef]
  35. Yan, J.; Zhang, Q.; Yin, P. RNA editing machinery in plant organelles. Sci. China Life Sci. 2018, 61, 162–169. [Google Scholar] [CrossRef]
  36. Contreras-Díaz, R.; Arias-Aburto, M.; van den Brink, L. Characterization of the complete chloroplast genome of Zephyranthes phycelloides (Amaryllidaceae, tribe Hippeastreae) from Atacama region of Chile. Saudi J. Biol. Sci. 2022, 29, 650–659. [Google Scholar] [CrossRef]
  37. Li, H.; Xiao, W.; Tong, T.; Li, Y.; Zhang, M.; Lin, X.; Zou, X.; Wu, Q.; Guo, X. The specific DNA barcodes based on chloroplast genes for species identification of Orchidaceae plants. Sci. Rep. 2021, 11, 1424. [Google Scholar] [CrossRef]
  38. Green, B.R. Chloroplast genomes of photosynthetic eukaryotes. Plant J. 2011, 66, 34–44. [Google Scholar] [CrossRef] [PubMed]
  39. Liu, Q.; Li, X.; Li, M.; Xu, W.; Schwarzacher, T.; Heslop-Harrison, J.S. Comparative chloroplast genome analyses of Avena: Insights into evolutionary dynamics and phylogeny. BMC Plant Biol. 2020, 20, 406. [Google Scholar] [CrossRef]
  40. Xue, S.; Shi, T.; Luo, W.; Ni, X.; Iqbal, S.; Ni, Z.; Huang, X.; Yao, D.; Shen, Z.; Gao, Z. Comparative analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina. Hortic. Res. 2019, 6, 89. [Google Scholar] [CrossRef] [Green Version]
  41. Guo, Y.-Y.; Yang, J.-X.; Bai, M.-Z.; Zhang, G.-Q.; Liu, Z.-J. The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions. BMC Plant Biol. 2021, 21, 248. [Google Scholar] [CrossRef] [PubMed]
  42. Trofimov, D.; Cadar, D.; Schmidt-Chanasit, J.; Rodrigues de Moraes, P.L.; Rohwer, J.G. A comparative analysis of complete chloroplast genomes of seven Ocotea species (Lauraceae) confirms low sequence divergence within the Ocotea complex. Sci. Rep. 2022, 12, 1120. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, S.; Yang, C.; Zhao, X.; Chen, S.; Qu, G.-Z. Complete chloroplast genome sequence of Betula platyphylla: Gene organization, RNA editing, and comparative and phylogenetic analyses. BMC Genom. 2018, 19, 950. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Han, C.; Ding, R.; Zong, X.; Zhang, L.; Chen, X.; Qu, B. Structural characterization of Platanthera ussuriensis chloroplast genome and comparative analyses with other species of Orchidaceae. BMC Genom. 2022, 23, 84. [Google Scholar] [CrossRef]
  45. Alzahrani, D.; Albokhari, E.; Yaradua, S.; Abba, A. Complete chloroplast genome sequences of Dipterygium glaucum and Cleome chrysantha and other Cleomaceae Species, comparative analysis and phylogenetic relationships. Saudi J. Biol. Sci. 2021, 28, 2476–2490. [Google Scholar] [CrossRef]
  46. Volk, G.M.; Henk, A.D.; Baldo, A.; Fazio, G.; Chao, C.T.; Richards, C.M. Chloroplast heterogeneity and historical admixture within the genus Malus. Am. J. Bot. 2015, 102, 1198–1208. [Google Scholar] [CrossRef]
  47. Nikiforova, S.V.; Cavalieri, D.; Velasco, R.; Goremykin, V. Phylogenetic analysis of 47 chloroplast genomes clarifies the contribution of wild species to the domesticated apple maternal line. Mol. Biol. Evol. 2013, 30, 1751–1760. [Google Scholar] [CrossRef] [Green Version]
  48. Coart, E.; Van Glabeke, S.; De Loose, M.; Larsen, A.S.; RoldÁN-Ruiz, I. Chloroplast diversity in the genus Malus: New insights into the relationship between the European wild apple (Malus sylvestris (L.) Mill.) and the domesticated apple (Malus domestica Borkh.). Mol. Ecol. 2006, 15, 2171–2182. [Google Scholar] [CrossRef]
  49. Wang, S.; Yang, J.; Zhang, B.; Shang, F.; Wang, H. The complete chloroplast genome sequence of Malus halliana (Rosaceae), an important ornamental plant. Mitochondrial DNA B Resour. 2020, 5, 2155–2156. [Google Scholar] [CrossRef]
Genes 13 00560 g001 550
Figure 1. The map of complete chloroplast genome in M. zhaojiaoensis. Different colored blocks on the outer ring represent different genes, while the shadow on the inner ring indicates GC content.
Figure 1. The map of complete chloroplast genome in M. zhaojiaoensis. Different colored blocks on the outer ring represent different genes, while the shadow on the inner ring indicates GC content.
Genes 13 00560 g001
Genes 13 00560 g002 550
Figure 2. The comparison of SSRs and INEs of chloroplast genomes in three Malus species. (A): Numbers of SSRs with different length; (B): Frequency of different types of SSRs; (C): Numbers of different types of INEs; (D): Frequency of INEs with different length.
Figure 2. The comparison of SSRs and INEs of chloroplast genomes in three Malus species. (A): Numbers of SSRs with different length; (B): Frequency of different types of SSRs; (C): Numbers of different types of INEs; (D): Frequency of INEs with different length.
Genes 13 00560 g002
Genes 13 00560 g003 550
Figure 3. RSCU characteristics of chloroplast genome in M. zhaojiaoensis.
Figure 3. RSCU characteristics of chloroplast genome in M. zhaojiaoensis.
Genes 13 00560 g003
Genes 13 00560 g004 550
Figure 4. RNA editing sites of CDS in M. zhaojiaoensis cp genome. Percentages represent the proportion of potential editing sites in different genes.
Figure 4. RNA editing sites of CDS in M. zhaojiaoensis cp genome. Percentages represent the proportion of potential editing sites in different genes.
Genes 13 00560 g004
Genes 13 00560 g005 550
Figure 5. The similarity of the cp genomes of M. zhaojiaoensis and three Malus species. Different colors represent different comparison areas with different scoring rates. Score/max ratios were colored with blue (≤0.25), green (≤0.50), orange (≤0.75), and red (>0.75).
Figure 5. The similarity of the cp genomes of M. zhaojiaoensis and three Malus species. Different colors represent different comparison areas with different scoring rates. Score/max ratios were colored with blue (≤0.25), green (≤0.50), orange (≤0.75), and red (>0.75).
Genes 13 00560 g005
Genes 13 00560 g006 550
Figure 6. Hotspots of variation in chloroplast genomes of Malus species. The abscissa represents the region scope of the cp genomes, and the ordinate shows the similarity ratio of the cp genomes. CNS in the legend refers to non-coding sequences.
Figure 6. Hotspots of variation in chloroplast genomes of Malus species. The abscissa represents the region scope of the cp genomes, and the ordinate shows the similarity ratio of the cp genomes. CNS in the legend refers to non-coding sequences.
Genes 13 00560 g006
Genes 13 00560 g008 550
Figure 8. The evolutionary trees of M. zhaojiaoensis based on chloroplast sequence. M. zhaojiaoensis is highlighted with solid blue circle and hollow blue circle respectively, and the number of the scale reflects the genetic distance. (A) Neighbor-joining tree; (B) Maximum likelihood tree.
Figure 8. The evolutionary trees of M. zhaojiaoensis based on chloroplast sequence. M. zhaojiaoensis is highlighted with solid blue circle and hollow blue circle respectively, and the number of the scale reflects the genetic distance. (A) Neighbor-joining tree; (B) Maximum likelihood tree.
Genes 13 00560 g008
Genes 13 00560 g009 550
Figure 9. Phylogenetic relationships based on single-copy cp genes matK and rbcL in 10 Malus species. M. zhaojiaoensis is highlighted with red underline in (A,B).
Figure 9. Phylogenetic relationships based on single-copy cp genes matK and rbcL in 10 Malus species. M. zhaojiaoensis is highlighted with red underline in (A,B).
Genes 13 00560 g009
Table
Table 2. Distribution and number of introns in M. zhaojiaoensis cp genome.
Table 2. Distribution and number of introns in M. zhaojiaoensis cp genome.
GeneTypeLocationStrandStartEndExon SizeIntron Number
rps16CDSLSCReverse524763782701
atpFCDSLSCReverse12,44413,7305551
rpoC1CDSLSCReverse21,86524,65120461
ycf3CDSLSCReverse45,42347,3825072
rps12CDSLSC-IRBReverse73,406102,7143722
rps12CDSLSC-IRAForward73,406146,1533722
clpPCDSLSCReverse73,68475,7465912
petBCDSLSCForward78,68480,1286481
petDCDSLSCForward80,31981,5254831
rpl16CDSLSCReverse85,09486,4894081
rpl2CDSIRBReverse88,26189,7718251
ndhBCDSIRBReverse98,873101,07415331
ndhACDSSSCReverse124,864127,08010921
ndhBCDSIRAForward146,995149,19615331
rpl2CDSIRAForward158,298159,8088251
trnK-UUUtRNALSCReverse17034271721
trnG-UCCtRNALSCForward90929860711
trnL-UAAtRNALSCForward50,58151,181851
trnV-UACtRNALSCReverse54,75855,424751
trnI-GAUtRNAIRBForward106,519107,538721
trnA-UGCtRNAIRBForward107,603108,482731
trnA-UGCtRNAIRAReverse139,587140,466731
trnI-GAUtRNAIRAReverse140,531141,550721
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, X.; Wang, D.; Gao, N.; Han, Y.; Wang, X.; Shen, X.; You, C. Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species. Genes 2022, 13, 560. https://doi.org/10.3390/genes13040560

AMA Style

Wang X, Wang D, Gao N, Han Y, Wang X, Shen X, You C. Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species. Genes. 2022; 13(4):560. https://doi.org/10.3390/genes13040560

Chicago/Turabian Style

Wang, Xun, Daru Wang, Ning Gao, Yuepeng Han, Xiaofei Wang, Xiang Shen, and Chunxiang You. 2022. "Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species" Genes 13, no. 4: 560. https://doi.org/10.3390/genes13040560

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop